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Abstract. We study a generic problem of dissipative quantum mechanics, a small 
local quantum system with discrete states coupled in an arbitrary way (i.e. not 
necessarily linear) to several infinitely large particle or heat reservoirs. For both 
bosonic or fermionic reservoirs we develop a quantum field-theoretical diagram- 
matic formulation in Liouville space by expanding systematically in the reservoir- 
system coupling and integrating out the reservoir degrees of freedom. As a result 
we obtain a kinetic equation for the reduced density matrix of the quantum sys- 
tem. Based on this formalism, we present a formally exact perturbative renormal- 
ization group (RG) method from which the kernel of this kinetic equation can be 
calculated. It is demonstrated how the nonequilibrium stationary state (induced 
by several reservoirs kept at different chemical potentials or temperatures), arbi- 
trary observables such as the transport current, and the time evolution into the 
stationary state can be calculated. Most importantly, we show how RG equations 
for the relaxation and dephasing rates can be derived and how they cut off gener- 
ically the RG fiow of the vertices. The method is based on a previously derived 
real-time RG technique |1I2I3I4] but formulated here in Laplace space and general- 
ized to arbitrary reservoir-system couplings. Furthermore, for fermionic reservoirs 
with flat density of states, we make use of a recently introduced cutoff scheme on 
the imaginary frequency axis [5] which has several technical advantages. Besides 
the formal set-up of the RG equations for generic problems of dissipative quan- 
tum mechanics, we demonstrate the method by applying it to the nonequilibrium 
isotropic Kondo model. We present a systematic way to solve the RG equations 
analytically in the weak-coupling limit and provide an outlook of the applicability 
to the strong-coupling case. 



1 Introduction 



General remarks. Dissipative quantum mechanics is a fundamental field in theoretical physics 
combining the concepts of quantum mechanics and nonequilibrium statistical mechanics [6j . The 
aim is to develop a microscopic description of how a small quantum-mechanical system in con- 
tact with large reservoirs (see FigUfor a sketch of the system) evolves into a stationary state 
and what this stationary state looks like. In equilibrium statistical mechanics for large quan- 
tum systems in contact with a single reservoir, the strength and nature of the system-reservoir 
interaction is not important, since it is a surface effect and negligible compared to the bulk 
of the quantum system. Therefore, in this case, depending on which conserved quantities are 
exchanged, the stationary state is always a canonical or grandcanonical ensemble. For quantum 
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Fig. 1. A small quantum system coupled to several infinitely large reservoirs via energy and/or particle 
exchange. The reservoirs are characterized by temperatures Ti and chemical potentials /ii. 



systems in contact with several reservoirs which are kept at different temperatures or chemical 
potentials (inhomogeneous boundary conditions), the stationary state is a nonequilibrium state 
and possibly current-carrying. However, if the quantum system is large, again only the inter- 
actions in the bulk are important which can be treated perturbatively via standard quantum 
Boltzmann equations ^ . After a short crossover time local equilibrium establishes and the fur- 
ther time evolution can be described by hydrodynamic equations [8]. These standard tools of 
nonequilibrium statistical mechanics break down for small quantum systems for several reasons. 
First the system-reservoir coupling is no longer a negligible surface effect but can change the 
states on the system considerably via quantum fluctuations. Even for weak coupling renormal- 
ization of the coupling parameters and the level positions of the quantum system can occur at 
low temperatures. Strinkingly new effects can occur for strong coupling such as a localization 
transition in the spin-boson model [6] or resonant transmission for a local spin coupled via ex- 
change to reservoir spins (the so-called Kondo effect) j9|10j . Secondly, it becomes energetically 
difficult to put several electrons on the quantum system due to the large capacitive interaction 
Ec ~ e^/L (L being the length of the quantum system), the so-called charging energy. Typi- 
cal experimental values for semiconductor quantum dots, metallic islands or carbon nanotubes 
are Ec ~ 1 — lOK and much larger than typical temperatures T ~ 10 — lOOmK (for single 
molecules coupled to leads the charging energy can be even larger reaching typical atomic val- 
ues). A simple perturbative expansion in the interaction is no longer possible and usual quantum 
Boltzmann equations can not be used. Furthermore, for low-dimensional systems, the Coulomb 
interaction can lead to completely new physical phenomena such as Luttinger liquid behaviour 
in 1-dimensional quantum wires [111 . Concepts like local equilibrium are not applicable. Phase 
coherence is maintained over the whole system size and the quantum system acts like a scat- 
tering region for electrons entering and leaving the system rather than a region where particles 
can relax, dephase and equilibrate. 

For these reasons, new theoretical tools have been developed to understand relaxation, de- 
phasing, and nonequilibrium quantum transport through small quantum systems coupled to 
external reservoirs. For noninteracting systems, the standard tool is the Landauer-Biittiker for- 
malism |12j where the particle current is expressed by the scattering matrix together with the 
occupation of the scattering waves determined by the chemical potentials of the reservoirs. In 
this case, it is possible to consider arbitrary coupling between reservoirs and quantum system 
and the coherent properties of the quantum system are fully taken into account. For interacting 
systems (or systems with spin degrees of freedom) , the situation is much more complicated and 
no unique analytical or numerical formalism is available which can cover all regimes of interest. 
Numerical methods for nonequilibrium are currently been developed, such as time-dependent 
density matrix renormalization group (TD-DMRG) |13j . time-dependent numerical renormal- 
ization group [14j , numerical renormalization group at finite bias voltage using scattering waves 
[15], quantum Monte Carlo with complex chemical potentials [16], and iterative path-integral 
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approaches [17]. Exact solutions using scattering Bethe-ansatz are available for the resonant 
level model |18j . Concerning analytical methods two perturbative approaches are commonly 
used, depending on whether one expands in the interaction parameter inside the quantum 
system or in the reservoir-system coupling. Expanding in the interaction has the advantage 
that the unperturbed part of the Hamiltonian is quadratic in the field operators and standard 
Keldysh-Green's function techniques can be applied [19120] . Furthermore, rather large quantum 
systems can be treated since the effort scales with the number of single-particle levels rather 
than with the number of many-particle states. However, this method has its limitations since 
for typical quantum dots the Coulomb interaction is the largest energy scale of the problem. In 
contrast, expanding in the reservoir-system coupling has the advantage that arbitrary interac- 
tion strength on the quantum sytem can be treated. Furthermore, the reservoir-system coupling 
is often tunable in experiments and in most cases the lowest energy scale. Therefore it seems 
reasonable to expand around the point where reservoirs and quantum system are decoupled. In 
this case the unperturbed part of the Hamiltonian contains the full interacting quantum system 
and standard Green's function techniques can not be applied (a generic problem for all strongly 
correlated systems). One way out of this problem is the use of slave particle techniques where 
the interacting system is expressed in a quadratic form using creation and annihilation operators 
of many-particle states |21|22j . Standard Keldysh-Green's function methods can then be used 
by expanding in the reservoir-system coupling [23| . However, technical complications arise due 
to an additional constraint for the slave particle number and, most importantly, diagrammatic 
approximations are often doubtful due to the unphysical nature of the slave particles. So even 
the noninteracting case is quite nontrivial [52] and vertex corrections are essential to obtain the 
relaxation and dephasing rates for the physical particles [241 . In contrast, the most standard 
method of dissipative quantum mechanics is to integrate out only the noninteracting reservoirs 
and describe the dynamics of the reduced density matrix of the quantum system via a kinetic 
equation. This can be achieved by projection operator techniques in Liouville space |25j or via 
path-integral methods ^6 26 . Most recently, a quantum field-theoretical version of this strategy 
has been developed in Refs. [28 29 30], see Ref. [27] for a review. The advantage is that Wick's 
theorem is used to integrate out the reservoirs and an exact diagrammatic representation of 
the kernel determining the kinetic equation is obtained. This allows for a direct calculation of 
this kernel in terms of irreducible diagrams whereas with projection operator techniques the 
calculation is complicated by cancellations of reducible expressions. Furthermore, the usage of 
representations similiar to Feynman diagrams simplifies the implementation of renormalization 
group ideas. 

RG for the Kondo model. The analytical methods expanding either in the interaction or 
the reservoir-system coupling can of course be used to calculate all physical quantities of inter- 
est in perturbation theory. In this way many time-dependent and nonequilibrium phenomena 
have been described in various fields, such as spin boson models, quantum optics, mesoscopic 
systems, quantum information theory, etc. However, perturbation theory is often plagued by 
various diverging terms in higher orders if some low energy scale like temperature, voltage, 
magnetic field, etc. becomes too small. In this limit perturbation theory becomes ill-defined 
and the natural question arises whether renormalization group methods can be generalized to 
the nonequilibrium situation to resum the original perturbation theory in an appropriate way 
so that it becomes well-defined again. To be specific let us introduce a simple but nontrivial 
example, the nonequilibrium Kondo model. This model is currently one of the basic unsolved 
problems of nonequilibrium condensed matter physics and serves as an important benchmark 
model for theoreticians to test the applicability of their nonequilibrium techniques. The model 
consists of the most simplest fermionic quantum system one can imagine, namely a spin-^ 
system, which interacts with two fermionic reservoirs (being kept at temperature T and two 
different chemical potentials /^l and ^r) via exchange processes, see Fig[2| The Hamiltonian 
for the reservoir-system coupling reads 
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Fig. 2. A spin-i quantum system coupled via exchange to two reservoirs. Jl = Jll and Jji 
involve exchange between the spins of the left/right reservoir with the local spin, and Jlr 
transfers a particle from one reservoir to the other during the exchange process. 
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Here, S_ denotes the spin operator of the quantum system, a are the Pauli matrices, and a, are 
the fermionic annihilation and creation operators for the particles in the reservoirs characterized 
by reservoir index a = L, R, spin cr =1, J, and state index k. The exchange couplings are denoted 
by Jaa' ■ For a symmetric system, there are only two different exchange couplings 

Jd = Jll = Jbm , Jnd — Jlr = Jrl ■ (2) 

The diagonal exchange coupling Jd characterizes spin exchange which involves only the spin of 
one reservoir, whereas the nondiagonal coupling Jnd describes a process where one particle is 
transferred between the reservoirs (see e.g. Ref. [4J for a derivation of this form of the interaction 
via a standard Schrieffer- Wolff transformation from a conventional tunneling Hamiltonian). 
The latter process leads to a particle current at finite bias. For simplicity we consider only 
the isotropic case without magnetic field. The reservoirs are characterized by a noninteracting 
Hamiltonian 

a=L,R a=L,R 

and the statistics is given by an equilibrium grandcanonical distribution 

Pres = W Pa . Pa = , (4) 

a=L,R " 

where Tl = Tfi = T denotes the temperature and p,L = —fJ.R. = V/2 the chemical potentials of 
the reservoirs {V is the bias voltage and we use units e = Ti = 1). 

Even in equilibrium V = the Kondo model is a highly nontrivial model and numerous 
many-body techniques have been used to study its properties, see e.g. [31] for a review. We 
summarize here shortly its basic properties and consider the simplest case J = Jd = Jnd- 
Higher-order perturbation theory in J for transition rates (or the linear conductance) leads 
generically to logarithmic divergencies ~ {J \n{D / Ac))" , where D denotes the bandwidth of 
the reservoirs and ylc ~ T is the low-energy scale. For the linear conductance, the logarithmic 
terms start in 0{J^) 

G = Go^ j2(l + 4Jln(DM,)) , (5) 

with Go ~ jh. To resum the most divergent terms in each order of perturbation theory (the 
so-called leading-order analysis), poor man scaling methods have been developed |32| where the 
high energy scales of the reservoirs are successively integrated out. If A denotes the effective 
bandwidth of the reservoirs (with initial value D) , an infinitesimal reduction A A — dA is 
compensated by a renormalization of the exchange coupling J while keeping the scattering t- 
matrix invariant (in leading order). As a result one finds the so-called poor man scaling equation 
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G — Gq — j— — Go 1 2/m;m \ • (^) 



(valid for J < 1) 

dJ 2J2 

dl = -— ' 

with the solution 

-^(^^-sh^ ' ' (^) 

where Tr- denotes the Kondo temperature which is an invariant of the RG equation ^ . In the 
antiferromagnetic case J > 0, the effective coupling diverges at yl = indicating a complete 
screening of the local spin by the reservoir spins. This leads to a spin-singlet ground state and 
it can be shown that the remaining potential scattering terms lead to unitary conductance 
(Kondo effect). The Kondo effect has been measured for semiconductor quantum dots, carbon 
nanotubes, and molecules [3] (for theoretical works see Refs. [10 23 30 or Ref. 47J for a review). 
This is the so-called strong-coupling regime where the perturbative RG equation ([6]) is no longer 
valid. However, the poor man scaling equation has to be cut off by the low-energy scale Ac ^ T 
and the strong coupling regime can not be reached for Ac ^ Tk- The weak coupling regime is 
defined by Jc = J {Ac) <Si 1, where perturbation theory has to be carried out in the renormalizcd 
coupling Jc with an effective bandwidth given by Ac. Replacing J Jc and D Ac in 
one obtains 

37r2 ^, _ 37r2 1 

16" In^ (T/Tk) 

Since Jc is logarithmically increased compared to the bare coupling Jq = J{D), the onset of 
the Kondo effect is indicated by a logarithmic enhancement of the conductance as function 
of temperature. Note that simple perturbation theory in the original coupling Jq can already 
break down in this regime since the two conditions 

JoHD/Ac) = ^-Joln{Ac/TK) ^ 0(1) and J, = 2ln{Ac/TK) ^ ^ 

can easily be fulfilled, provided that 1 ^ ln(ylc/TR-) ^ 1/Jo (although the typical experimental 
situation is rarely strictly in this regime, it provides a well-defined regime where perturbative 
renormalization group methods can be tested). 

In nonequilibrium V ^ and V > T, the situation is not so clear. The standard poor man 
scaling procedure suggests that the low-energy scale Ac is given by the maximum of T and 
y, i.e. the voltage serves as a cutoff parameter in the same way as temperature. Alternatively 
it was conjectured in Ref. [M] that Jd is cut off by T (since it involves only the spins of one 
reservoir) and Jnd is cut off by max{T, V}, opening up the possibility of a strong-coupling fixed 
point for V ^ Tk ^ T. Finally, in Refs. [35|36j it was proposed that also spin relaxation and 
dephasing rates can cut off the RG flow. For the isotropic Kondo model in the absence of a 
magnetic field the relaxation and dephasing rates are the same and are given by the Korringa 
rate 

r = TT J2 max{r, V"} . (10) 

This has the consequence that F ^ Tk if F ^ Tk (note that Tk is exponentially small for 
J ^ 1 whereas F scales quadratic with J) and a strong-coupling fixed point can not be reached 
for V ^ Tk S> T. In this paper we will provide a microscopic nonequilibrium RG approach 
clarifying all these issues. Indeed, it turns out that the Korringa rate cuts off the RG flow of Jd 
and Jnd- Roughly the poor man scaling equation ([6|) has to be replaced by the RG equations 
(see Eqs. (HH and (gHH) of Sec. ^ 

with Ot = 0{A — T), where we have omitted additional frequency dependencies of the vertices 
(see Sec. 15.31 for more details where also the RG equation for F is shown). As one can see. 
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temperature and Korringa rate cut off all terms on the r.h.s of the RG equation whereas the 
voltage is only present in certain terms. This shows that nonequilibrium induces new features 
into the RG equations. Neither Jd nor Jnd are completely cut off by the voltage and from these 
RG equations alone there is no justification to cut off the RG flow by the voltage \i V > T . 
Instead, the RG flow should be cut off by F leading to additional ln( Jc) contributions for the 
exchange couplings. However, as will be discussed in detail in Sec. 15. 3[ there is an additional 
RG equation for the current rate which is cut off by the voltage and not by F . As a consequence 
the cutoff parameter which enters into the conductance is indeed the voltage and we obtain 
precisely the result ([8]) with temperature replaced by bias voltage. However, it is important 
to recognize that the fact that F does not enter the final result for the conductance does not 
mean that it can be neglected. For T = and V 3> Tk: the RG equations pT|) and would 
diverge for yl —> if the cutoff parameter F were not present, leading also to a divergence of 
the conductance since the RG equation for the current rate is cut off smoothly by the voltage, 
see Sec. 15.31 Therefore, it is an important issue for nonequilibrium RG methods to incorporate 
the physics of relaxation and dephasing rates into the RG equations. 

The situation becomes even more interesting if an additional energy scale A is present, 
such as magnetic field, external driving frequency, level spacing, charge excitation energy etc. 
In this case, it turns out that additional logarithmic terms can occur which are generically 
of the form \n.{D/\n^ — A\) with n being an integer number, see Sec. 14.31 and Refs. [43l33j . 
At V ^ 2A/n these logarithmic terms diverge corresponding to certain resonance positions. 
Relaxation and dephasing rates will also cut off these divergencies. Thus, their effect is two-fold: 
relaxation and dephasing rates cut off the RG flow for the vertices and also the RG flow for 
all physical observables, such as the current rate, susceptibilities, occupation probabilities, etc. 
Furthermore, an interesting question is what happens in the strong-coupling regime, where the 
RG flow of the relaxation and dephasing rates is cut off by themselves. It is expected that 
they saturate to the Kondo temperature and prevent the RG flow from diverging. This issue is 
discussed in Sec. 15.31 

RG methods. Motivated by these considerations and due to the progress in experimental 
techniques to study quantum transport through small nanodevices, there is an increased interest 
in the development of renormalization group methods for nonequilibrium systems (which are 
either driven out of equilibrium with an external bias or are prepared in an out of equilibrium 
initial state). The flrst proposal for a nonequilibrium RG technique, commonly called real-time 
RG (RTRG), was provided in Ref. [I] for problems of dissipative quantum mechanics, see Ref. ^ 
for a tutorial introduction. The method is perturbative in the reservoir-system coupling and 
is based on the diagrammatic kinetic equation approach reviewed in Ref. [27' . Instead of the 
bandwidth, the relative time of the reservoir Green's function was used as a cutoff parameter and 
the kernel of the kinetic equation was kept invariant during the RG flow by setting up a formally 
exact hierarchy of RG equations. The method has been applied to charge fluctuations in the 
single electron box [37j . transport through metallic quantum dots [38 39 , charge fluctuations in 
semiconductor quantum dots [1140] . the polaron problem [41 , the influence of acoustic phonons 
on transport through coupled quantum dots [42j , and to the dynamics of the spin boson problem 
[2]. Within all these applications, a linear coupling between the quantum system and the particle 
or heat reservoir was considered. In this case, it was sufficient to use the RG method with a cutoff 
defined in time space. For models like the nonequilbrium Kondo model, where spin fluctuations 
are important, a quadratic coupling involving one creation and one annihilation operator of the 
reservoirs occurs, see ([T]). To treat this case in a convenient way, a version of the RTRG method 
with a cutoff defined in real frequency space has been developed recently [4]. An essential 
ingredient of these techniques is the generation of non-Hamiltonian dynamics during the RG 
flow in terms of an effective Liouville operator of the quantum system with nonzero imaginary 
parts of its eigenvalues representing the relaxation and dephasing rates. This effective Liouville 
operator appears also in the RG equations of the vertices cutting them off at these energy scales. 
However, a problem of this approach is the fact that a certain time-ordering procedure has to 
be defined for the renormalized vertices which leads to the generation of terms which are not 
present in the original perturbation theory and have to be corrected by counter terms in higher 
orders. This leads to the problem that the cutoff of the RG flow by relaxation and dephasing 
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rates can only be seen in leading order and a proof that this effect happens in all orders is 
not possible (leaving the question open whether a strong-coupling fixed point is in principle 
possible). In addition, it is not unambigiously clear what the precise scale of the rates cutting 
off the RG flow is. Finally, concerning the numerical stability for solving the RG equations, it 
turns out that cutoff functions which are defined in real time or real frequency space are not 
the best choice. In this paper, we will show how these problems can be circumvented by some 
new technical tricks which have been used recently in Ref. [43| to calculate analytically the 
nonlinear conductance, the magnetic susceptibility, the spin relaxation and dephasing rates, 
and the renormalized g-factor up to 2-loop order for the nonequilibrium anisotropic Kondo 
model at finite magnetic field in the weak-coupling regime. The first technical step is to show 
that the RG approach can be set up purely in frequency space [44j using the diagrammatic 
rules of Ref. [27j . As a consequence the occurence of correction terms can be avoided and the 
rates determining the cutoff of the RG flow obtain the right scale. Secondly, it turns out that 
it is important to integrate out the symmetric part of the Fermi or Bose distribution function 
of the reservoirs before starting the RG flow. With this choice it is possible to show generically 
for all models of dissipative quantum mechanics that relaxation and dephasing rates cut off 
the RG flow in all orders of perturbation theory and within all truncation schemes of the RG 
equations. Another technical advantage is that the dependence on the Keldysh indices can be 
avoided, minimizing the effort to solve the RG equations. Finally, it is more convenient to 
define a cutoff not in real but in imaginary frequency space by introducing the cutoff into the 
Matsubara poles of the reservoir distribution function. This idea was first proposed in Ref. [5] in 
the context of nonequilibrium functional renormalization group within the Keldysh formalism. 
This choice improves the numerical stability considerably. Furthermore, for fermionic reservoirs 
with a flat density of states, it is possible to reformulate the RG equations in pure Matsubara 
space (all integrals over the real frequencies can be performed analytically) with the difference 
to equilibrium that a whole set of Matsubara axis occurs shifted by multiples of the chemical 
potentials of the reservoirs and by the real part of the Laplace variable determining the time 
evolution of the system. With this technique it is possible to provide a generic procedure how 
to solve the RG equations analytically in a well-controlled way in the weak-coupling regime, 
see Ref. [33]. Whether it is also applicable to the strong-coupling regime is an open question 
discussed in Sec. 15.31 for the special case of the nonequilibrium isotropic Kondo model. For 
later reference, we will call this version of nonequilibrium RG real-time RG in frequency space 
(RTRG-FS). 

Especially for the Kondo model, the development of nonequilibrium RG has been pioneered 
in Refs. [36,45j (for perturbation theory see Refs. [46 24, ; a poor man scaling approach is 
described in Ref. [32] )• In these works the slave particle approach was used in connection with 
Keldysh formalism and quantum Boltzmann equations. A real-frequency cutoff was used and the 
RG was formulated purely on one part of the Keldysh-contour disregarding diagrams connecting 
the upper with the lower branch. This procedure turns out to be sufficient for the Kondo model 
to calculate logarithmic terms in leading order but the cutoff by relaxation and dephasing 
rates can not be obtained in this way. In these works it was investigated for the first time how 
the voltage and the magnetic field cuts off the RG flow and how the frequency dependence 
of the vertices influences various logarithmic contributions ^ ln(^S^^^-^^) or ~ lii(-^) for the 
susceptibility and the nonlinear conductance. 

An alternative approach to RTRG-FS for combining relaxation and dephasing rates mi- 
croscopically with nonequilibrium RG was proposed in Ref. [35 using flow-equation methods 
[49] . Speciflcally for the isotropic Kondo model without magnetic field it was shown that the 
cutoff of the RG flow by the decay rate F occurs due to a competition of 1-loop and 2-loop 
terms on the r.h.s. of the RG equation for the vertex. This is a completely different picture 
compared to RTRG-FS, where the cutoff parameter F occurs already in the 1-loop terms as 
an additional term in the denominator of the resolvents, see the RG-equations (|lip and 
Thus, RTRG-FS is closer to conventional poor man scaling equations and the physics of relax- 
ation and dephasing rates occurs naturally as a resummation of a geometric series similiar to 
self-energy insertions in Green's function techniques. Nevertheless, the flow equation method is 
well-definied and controlled, and represents a technical alternative to RTRG-FS. 
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Nonequilibrium RG methods which expand perturbatively in the Coulomb interaction of the 
quantum system have also been developed recently I 5|52|53|54 ]. In these works, the Keldysh 
formalism has been combined with functional RG techniques known from quantum field theory 
[50 51j. It turns out that a real- frequency cutoff is not useful since it violates causality and 
leads to various numerical instabilites [53152] . For fermionic models in Id a more useful cutoff 
scheme on the imaginary frequency axis has been proposed [5], which is also very useful for 
RTRG-FS (see the discussion above). An interesting result was obtained that Luttinger liquid 
exponents can depend on the nonequilibrium distribution function of the quantum system [5j. 
For zero-dimensional systems (quantum dots) with a single spin-degenerate level coupled by 
tunneling to two reservoirs (the so-called single- impurity Anderson model), the situation is 
more complicated because there is no controlled truncation scheme for a perturbative RG in 
the Coulomb interaction U, especially in the interesting regime where the Coulomb interaction 
becomes larger than the energy scale F of the reservoir-system coupling. Furthermore, the 
Green's function in the absence of the Coulomb interaction has a pole in the complex plane 
with imaginary part F (corresponding to energy broadening of the dot level due to coupling to 
the reservoirs). This pole is close to the real-axis compared to the cutoff-parameter of the RG 
flow. As a consequence the cutoff procedure defined on the Matsubara axis is not sufficient for 
this problem and F itself was proposed as the flowing cutoff parameter [S^ . With this choice it is 
possible to study the nonequilibrium Anderson impurity-model up to values oiU ^ 2 — AF even 
in the Kondo regime T, y <C Tk [54] (the truncation scheme in this work neglects the 3-particle 
vertex but takes the important parts of the frequency-dependence of the 2-particle vertex into 
account). For U ^ F charge fluctuations are suppressed and the model is equivalent to the 
Kondo model with J ^ F/U . However, it is not yet possible to approach the limit J <C 1, where 
the Kondo temperature shows the typical exponential behaviour. It remains an interesting and 
open question whether this approach can be improved and a full solution of the nonequilibrium 
Anderson-impurity model including the nonequilibrium Kondo model can be obtained. 

Finally, we mention that in this introduction we have only summarized the nonequilibrium 
RG approaches relevant for problems of dissipative quantum mechanics. There is also a rapidly 
increasing interest in the development of nonequilibrium RG methods for bulk quantum systems 
motivated by the interesting possibilities to measure the time evolution of strongly correlated 
quantum systems in cold atom gases. For completeness we mention some of the most interesting 
recent developments, e.g. the study of quantum phase transitions in nonequilibrium |55|56j . the 
far-from-equilibrium quantum fleld dynamics of ultracold atom systems [57], and the time- 
evolution of fermionic system after initial interaction quenches [58 . 




Outline. The present paper concentrates on the scheme of RTRG-FS which is especially 
useful for a generic study of problems in dissipative quantum mechanics, where the expansion 
parameter is the reservoir-system coupling. We will describe the formal technique for a generic 
model but will always accompany the formalism by an illustration for fermionic models where 
charge, spin or orbital fluctuations dominate. Especially we will apply the formalism in all 
details to the nonequilibrium isotropic Kondo model. The paper is written for students with 
some basic knowledge of many-body theory. Besides second quantization and Wick's theorem 
nothing special is required. The paper is organized as follows. In Sec. [2] we will introduce the 
generic model and some specific examples. In Sec. [3| we present a nonequilibrium quantum 
field-theoretical diagrammatic description in Liouville space by introducing creation and anni- 
hilation superoperators acting in Liouville space (i.e. with an additional Keldysh-index) . This is 
especially useful for obtaining compact diagrammatic rules. Based on this diagrammatic expan- 
sion we present a derivation of a formally exact kinetic equation for the reduced density matrix 
of the quantum system and provide the rules to calculate observables. In Sec. HI we develop the 
nonequilibrium RG approach. We discuss the general idea of invariance, various ways to define 
the cutoff function, and set up simple rules to obtain the RG equations in arbitrary order in the 
coupling. We discuss their general properties and prove the central theorem that the RG fiow 
is always cut off by relaxation and dephasing rates. For fermionic models with spin or orbital 
fiuctuations, following Ref. [43], we will describe the generic scheme how to solve the RG equa- 
tions analytically in the weak-coupling regime. Finally, in Sec. [5] we apply the technique to the 
nonequilibrium isotropic Kondo model, show the results of Ref. [43] in 1-loop order, and present 
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preliminary results for the strong-coupling case. We close with a summary and an outlook in 
Sec. El 



2 The model 

2.1 Generic case 

Generically any Hamiltonian of a problem in dissipative quantum mechanics can be decomposed 
into a part for the reservoirs, the local quantum system, and the coupling between reservoirs 
and quantum system 

H = Hres + Hs + V = Ho + V . (13) 

The unperturbed part is denoted by 

Ho = Hres + Hs . (14) 

For the local quantum system Hs we make no assumption and represent it via its eigenstates 
\s) and eigenenergies Eg 

s 

In practice the eigenstates \s) can easily be found for a quantum system with a low number 
of single particle states including arbitrary interactions. In reality each quantum system has 
an infinite number of states. However, if the system is very small (as we assume) the single 
particle level spacing is very large and high-lying states can be neglected (for sufficiently small 
temperature and bias voltage). Therefore, it is often sufficient to take only a few single-particle 
levels into account and Hs can be easily diagonalized. Since we want to include the possibility 
that particles can be exchanged between reservoirs and quantum system, the states have to be 
classified according to their particle number. We denote by Ns the particle number for state |s). 
We include both possibilities that the particles can be bosons or fermions. E.g. bosonic particles 
can be realized in cold atom gases whereas fermionic particles (electrons) occur typically in 
nanoelectronic systems. 

The reservoirs are assumed to be noninteracting and infinitely large with a continuum den- 
sity of states. Therefore, we use a continuum notation for the creation and annihilation operators 
of the particles in the reservoirs and the reservoir Hamiltonian reads in second quantization 

Hres='^Ha= ^ j duj {oj + p,a) al{uj)a^{uj) , (16) 

OL i/=aa... 

with the continuum commutation relations for the field operators 

[a,{uj),al{J)]^ ^ 5,,:5{u^J) , (17) 

where the upper (lower) sign correponds to bosons (fermions) and [A, i3]zp = AB =F BA is the 
(anti-)commutator. v is a, discrete index which labels all remaining quantum numbers of the 
reservoir particles 

v = aa ... , (18) 

where, by convention, a is the index numerating the reservoirs (e.g. a ~ L, R = ± for two 
reservoirs) and a is the spin index (e.g. a =1,1= ± for a spin-i system). Further possibilities 
are orbital indices (if orbital symmetries are present), channel indices (for transverse modes), 
etc. The chemical potential of reservoir a is denoted by Ha, and u; is the energy of the reservoir 
state measured relative to fj,a (for phonon or photon baths, the chemical potential is absent). 
The reservoirs are assumed to be described by a grandcanonical distribution function 

Pres=l[pc. , p„ = ^e-^(«°-'^=^°) , (19) 
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where is the temperature of reservoir a. 

We note that the crossover from a discrete to a continuum notation for the reservoir field 
operators can be formally achieved by the definition 

a,y(uj) = —^^=y2 6{uj - ei,k+ tJ'a)a^k , (20) 

where fc is a discrete index labelling the states of the reservoirs, e,yfe is their corresponding 
energy, and Pu{oj) is the density of states in reservoir a (which can depend on energy and 
possibly (for ferromagnetic leads) on the spin index) . It is easy to show that with this definition 
the reservoir Hamiltonian (jl6p in the continuum form is equivalent to the discrete version 



Hres = ^ ^uk O'lk^i^k ■ (21) 
i^k 

Thereby, we have assumed that the relation between k and e^k is unique, otherwise different 
branches of the dispersion relation have to be distinguished and labelled by additional indices. 

For later convenience we introduce the following more compact notation for the various 
indices characterizing the reservoir field operators 

1 — — f \ ( al{uj) for rj = + ,„r,N 

l = r,^u; , ai = a,.(c.) = | ^'^^^ ior r, ^ - " (^2) 

Here, = ± indicates either a creation or annihilation operator. If no ambiguities occur we 
use 1 = rjVLij and 1' = rj'v'oj' ^ whereas for several indices we take 1 = rjiviuoi^ 2 = r]2t^2^2, etc. 
Furthermore we define by 

l = -r],iyuj , ai = (ai)^ , (23) 

the index corresponding to the hermitian conjugate. With these notations the commutation 
relation (|17p can be written in the compact form 

[ai,ay]^ = 5n, , (24) 

where 6ir stands for 

Sii' = S,j,j'S„^'S{uj - uj') . (25) 

In all expressions we sum (integrate) implicitly over the indices, i.e. we sum over the discrete 
part Tj and v, and integrate over the continuum variable uj. The reservoir Hamiltonian can be 
written in the compact form 

Hres = + Pa)5ii> , (26) 

and the reservoir correlation function (also called reservoir contraction) reads 

aiay = (aiai')p,,^, = (^ii' /a(^) = ^ii' |^| faiv^) , (27) 
where /+(w) = /a(w), /^(w) = 1 ± /a(w), and 

/oM = ^^/t! ^ 1 = T/-(-u;) (28) 

is the Bose (Fermi) function corresponding to temperature Tq, (note that the chemical potential 
fia does not occur in this formula because uj measures the energy relative to fia)- The upper 
(lower) case in ([27| refers to bosons (fermions) , a convention which we will use throughout this 
paper. 

For the coupling between reservoirs and quantum system we take the generic form 



V = —, 912. ..n : ai a2 . . . a„ : , (29) 
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where 512. ..n is an arbitrary operator acting on the quantum system, and we sum implicitly over 
n — 1,2, . . . and all variables r]i, Vi, uji, i = 1,2, ... ,n, which are contained in the formal notation 
i = •qi,Vi,LOi. We note that we explicitly exclude the case n = because this corresponds 
to an operator of the local quantum system which can be incorporated in Hs- The symbol 
: . . . : denotes normal-ordering of the reservoir field operators with respect to the reservoir 
correlation function, i.e. in any Wick-decomposition of an average over a sequence of reservoir 
field operators with respect to Pres, no contraction is allowed which connects field operators 
occuring in the same normal-ordered block. As a consequence the average over a normal-ordered 
block is defined as zero (: 0102 . . . a„ ') p^^^ = 0. The field operators within the normal-ordering 
can be arranged in an arbitrary order (up to a sign for fermions). The prefactor ^ accounts 
for all permutations of reservoir field operators which do not change the value of the diagrams 
in perturbation theory (see later) since the coupling vertex can always be chosen such that it 
is (anti)symmetric under exchange of two indices 

9l...i...j...n = ± ffl...j...i...n • (30) 

For simplicity, we assume here that either fermionic or bosonic field operators occur for the 
reservoirs. If both are present (e.g. for combinations of electronic particle reservoirs and phonon 
heat baths), one simply writes the coupling as a sum of several terms, each one being of the form 
([29|) . In principle also mixed couplings can be treated containing fermionic and bosonic field 
operators in the same term but we will not consider this case here because it just complicates 
the notation (in this case there is a factor ^^^^^^ in front of Eq. ipS]) . where ub {np) is the 
number of bosonic (fermionic) field operators) . 

The form ([^5]) of the coupling includes all cases of charge, spin, and energy transfer between 
reservoirs and quantum system and, for n > 1, also includes nonlinear coupling. The coupling 
vertex gi2...n describes the change of the state of the quantum system and is an operator. 
We have included the possibility that it depends in an arbitrary way on the frequencies of the 
reservoir field operators. Although such a general form (especially for n > 2) is hardly necessary 
for realistic models, this general form with all possible values for n has to be considered since 
the RG procedure described in Sec. |4] will generate an effective coupling which includes all these 
terms. 

The total operator V is of bosonic nature since the parity of fermion number must be 
conserved. However, for fermions and n odd, gi2...n in (|29p is of fermionic nature and the 
sequence of gi2...n and : 0102 . . . a„ : is important (the coupling V in the form ([29]) is even not 
hermitian in this case). To be precise, in the general case V should be written as 

^ \ ■.a„an-i...ai: gi2...n . (31) 

With this choice it can easily be shown that V is hermitian, provided the coupling vertex fulfils 
the condition 

(.912. ..n)^ = gn...l ■ (32) 

Furthermore, it is shown in Appendix A that the form (|31[) is only important for fixing the 
value of the coupling vertex gi...n for a concrete model. After having defined gi...n in this 
way, one can proceed with the much simpler form (j29p and disregard all sign factors emerging 
from interchanging fermionic dot and reservoir operators. The reason is that by calculating the 
reduced density matrix of the quantum system or averages of observables, only expressions have 
to be evaluated where an average over the reservoir distribution p^es is taken, for details see 
Appendix A. Thus, in the following we work with the simpler form (j29|l and commute dot and 
reservoir operators. This simplifies the problem of sign factors considerably. 

Finally, we mention that it is sometimes more convenient to include the density of states of 
the reservoirs into the reservoir contraction (|27p . In this case we use the definition 



bi = VPi^M ai 



(33) 
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for the reservoir field operator and obtain for the contraction 

hbv = (6i6i')p,.e. = Sn'PA^)m^) ■ (34) 
Correspondingly, the coupling V is written in the form 

V ^ ^ gi2...n -hh ■■■ bn : , (35) 
where dot and reservoir operators commute, and 

^ = ^ {m-^--Vn] 5l2...n (36) 

for the determination of gi2...n- 
2.2 Specific examples 

Here, we present some specific and experimentally relevant examples for charge, spin, orbital, 
or energy fluctuations induced by the coupling between reservoirs and quantum system. 

Charge fluctuations. Electronic quantum transport through nanodevices or quantum dots 
is characterized by charge fluctuations and is conveniently described by a tunneling Hamilto- 
nian, i.e. the coupling is of the form 

^ = EEE*^r«L.c.; + (/i.e.) . (37) 

aa k I 

Here, a and a are the reservoir and spin indices, respectively, and I is an index for an arbitrary 
single-particle state basis of the quantum system, f^f are the tunneling matrix elements and Ca-i 
annihilates a particle on the quantum system in state I with spin a. In the continuum notation 
and with v = aa, we obtain 

y = ^ J du; rj a,j^{uj) g,^^{uj) = 771 ai gi , (38) 
with arip{uj) defined by ((22|l and (|20p. and the coupling vertex is given by 

where (w) = fj^i with uj = e^k — ^J'a■ 

To obtain dimensionless coupling vertices, one can define for each v and lu some reference 
tunneling matrix element f^itu) = (independent of the state index / of the quantum system), 
with a corresponding level broadening parameter r^{uj) defined by 

r,{uj) =. 27:p,{Lu)\e{co)\^ = 27r^ 6{uj - e,k + f^») ■ (40) 

fc 

With this reference energy, we define the dimensionless reservoir field operators 

i Efc Kfcl '^(^ - fi-fc + Ma) fli-fc for = - 

and the corresponding dimensionless coupling vertex 

. ^ / j:itl{^)/r{^)\caL for 77 = + 
- {J2^tri^r/\t^iu;)\cl for^ = - ' 



(42) 
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such that we obtain again the form (|38p with a — > &. The reservoir contraction with respect to 
the 6-operators reads 

Ml' ^ (&i6i'>p... = S,y ^r,{u) mco) . (43) 

Spin and orbital fluctuations. If the charging energy of a quantum dot is very large 
(compared to temperature and bias voltage), the gate voltage determining the position of 
the charge excitation energies of the quantum dot relative to the electrochemical potentials 
of the reservoirs can be adjusted such that charge fluctuations are suppressed. In this case, 
spin and orbital fluctuations dominate transport. The elementary processes are given by an 
electron tunneling on (off) the quantum system, occupying a virtual intermediate state, and 
in a second step tunneling off (on) the quantum system. As shown in detail in Rcf. [4J, the 
virtual intermediate state (which has a very high energy due to the large charging energy) can 
be integrated out via a standard Schrieffer- Wolff transformation on a Hamiltonian level, and 
the result is a coupling of the form 

V = : aiQv : ■ (44) 

In most cases (except for systems with negative charging energies [59J), one creation and anni- 
hilation operator is present, i.e. rj = —77' (note that the forms (|29p and (|3ip are equivalent in 
this case). 

Whereas (|44p describes a generic model for spin and/or orbital fluctuations on a quantum 
dot, a special case is the Kondo model ([T]) discussed in the introduction. In this case, only 
fluctuations of a spin-^ of the quantum system are considered and, comparing (|44p with ([1]), 
we obtain for the dimensionlcss coupling vertex in the continuum form 

such that the antisymmetry and hermiticity conditions 

511' = -3i'i , glr = ffi'i (46) 

are fulfilled, compare l|30p and ([5^ . 

Since, due to the Schrieffer- Wolff transformation, high-lying charge excitations have already 
been integrated out in the model (jH]), one has to consider the reservoirs with a finite bandwidth 
D (of the order of the charging energy). Thus, all frequencies lo have to be smaller than D, 
which can be achieved by various cutoff functions. In this paper, we will use a Lorentzian cutoff 
defined by the function 

Most elegantly, we introduce this function by a modification of the reservoir contraction ([27]) 

-> Sn' p{Lo) moj) (48) 

Although this is not done here by a formal redefinition of the field operators, the effect of this 
modification of the reservoir correlation function is that all frequencies are suppressed if they lie 
above the band cutoff D. This will become clear within the diagrammatic expansion described 
in Sec. [31 where the reservoir degrees of freedom are integrated out. 

Energy fluctuations. Within the traditional field of dissipative quantum mechanics, a 
heat bath or energy-exchange with the environment is considered, mostly within the so-called 
spin-boson model defined by the Hamiltonian 

Hres ^'^LOq a\aq , (49) 
9 
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Hs = -j<Jx + ^iz , (50) 
^= II + • (51) 

9 

Here, the quantum system consists of a two-level system with tunneling coupling A between the 
two levels and bias e determining the detuning. The bath consists of a set of harmonic oscillators 
(phonons) with bosonic creation and annihilation operators aj, a^. The phonon energy is > 
and the bath couples linearly (via the spatial coordinate ~ + a^) to the two-level system. 
The coupling parameters 7^ and the density of states of the phonon modes are conventionally 
parametrized by the spectral function 

J{u;) = TT ^ 72 - cjg) = 27r a e""/-" . (52) 

n — describes the important ohmic case, whereas n > and n < are called the super-ohmic 
and sub-ohmic case. D is the bandwidth of the bath, and a describes the coupling parameter 
(which is dimensionless for the ohmic case). 

For our continuum notation, we define in analogy to (|4ip the dimensionless bosonic field 
operators 

\Y.q9qS{uJ-ujg)ag for 7] = - 
and obtain with the dimensionless coupling vertex 



the form 

for the coupling, and 



9± = Y (54) 
V ^ gih (55) 
hby = S,i,-J{Lo)r{cu) (56) 

TT 

for the bath contraction, where f^{uj) — f{^), / (^) = 1 + f{^) and f{uj) — \/{e^^ — 1) is 
the Bose function. Note that only positive frequencies oj are allowed since the phonon energies 
are positive. 

Formally, one can rewrite the expressions such that both positive and negative frequencies 
are allowed by using the definiton 

Since (i,,(ti>) = (i_,,(— w), we get 

/•OO poo POD 

/ dujbr,{Lo) — / duj dri{uj) — - / duj drj{uj) , (58) 

^ Jo ^ Jo 2 ^ 

and V can be written as (note that gi does not depend on 77) 

V = gidi , 51 = ^ • (59) 
Using f^{uj) = 1 + f(j-^) = ^fi^^), we get for the contraction in the new representation 

dTdv ^ 5n'-sign{u;)J{\uj\)r{uj) , (60) 

TT 

where all frequencies are allowed now. 
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3 Quantum field theory in Liouville space 

In this section we will develop a quantum field-theoretical diagrammatic expansion in Liouville 
space. Conventional quantum field theory in nonequilibrium is based on the Keldysh-formalism, 
where all field operators are integrated out via Wick's theorem or via Gaussian integrals within 
path integral formalism. This is not possible here because the unperturbed part of the Hamilto- 
nian, consisting of the reservoirs Hres and the quantum system iJg, contains arbitrary interac- 
tion terms in Hg. Therefore, this part is not quadratic and Wick's theorem does not apply. For 
this reason, only the degrees of freedom of the reservoirs can be integrated out. However, the 
coupling vertex gi,,,n remains an operator and the dynamics between the vertices is described 
by Hs- Therefore, a diagrammatic representation is obtained where the vertices are operators 
and one has to keep track of their time-ordering on the Keldysh-contour. In contrast to previ- 
ous versions of this procedure [3J , we introduce here a slightly different way of integrating out 
the reservoirs. Instead of expanding the time evolution on the Keldysh contour and applying 
Wick's theorem for the reservoir field operators, we will introduce quantum field superoperators 
acting in Liouville space for the reservoirs. This has the advantage that the two branches of 
the Keldysh contour can be taken together from the very beginning in a compact way. Sub- 
sequently, we will apply Wick's theorem for the quantum field superoperators. This procedure 
combines in an efficient way the advantages of Liouville superoperators (the traditional for- 
malism of dissipative quantum mechanics) and quantum field theory on the Keldysh contour 
(the common way to treat nonequilibrium systems). The traditional way of deriving kinetic 
equations by using projection operators in Liouville space [25] is a very formal procedure and 
much less useful. Although one obtains a very compact and analytic notation for the kernel of 
the kinetic equation without any need of a diagrammatic representation, the reservoirs are not 
integrated out and everything is hidden in formal projectors which finally have to be evaluated 
by decomposing all expressions into many artificial reducible terms where numerous cancella- 
tions occur. Integrating out the reservoirs from the very beginning has the advantage that the 
kernel can be defined via irreducible diagrams only and the calculations simplify considerably. 
Furthermore, it is shown in Sec.[4]that a systematic renormalization group method generates an 
additional frequency dependence of the Liouvillian and the vertices which can only be described 
within a diagrammatic representation. 

The main result of this section are the diagrammatic rules how to evaluate the kernel of the 
kinetic equation and averages of observables, which is the basis for the renormalization group 
method developed in Sec. [H 



3.1 Dynamics of the reduced density matrix 

General considerations. The dynamics of the total density matrix is governed by the von 
Neumann equation 

p{t) = ^i[H,p{t)]^ = ~iLp{t) , (61) 
where L is the so-called Liouville operator acting on usual operators A via the commutator 

LA = [H,A\^ . (62) 

Thus, in matrix notation, the Liouvillian is a super-matrix with four indices 

n'm' 

with 

Lnm,n'm' = ('^l ^ ( I"') I ) I™) = Hnn' Srnm' - 5nn' Hm' m ■ (64) 

Therefore the operators acting in Liouville space are also called superoperators. 

The initial state of the density matrix at time <o is assumed to be an independent product of 
an arbitrary part ps{to) for the quantum system and an equilibrium part pres for the reservoirs 

P{to) = Ps{to)Pres , (65) 
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where pres is the grandcanonical distribution defined in p9p . This means that we assume the 
reservoirs and the quantum system to be initially decoupled, and the coupling V is switched 
on suddenly at time to- 

The von Neumann equation can be formally solved by 

p{t) = e-^-f^(*-*°V(^o)e'-^(*-*°) = e-'-^(*-*«V(io) ■ (66) 

As a result wc obtain the following formal expression for the reduced density matrix of the 
quantum system 

ps{t) = TVres P{t) = Tr,e.e-*^(*-*«V(to)e^^(*-*'') = Tires e-''^^'-''^ psih) Pres , (67) 

where Tr^es denotes the trace over the reservoir degrees of freedom. 

The average of an arbitrary observable R can be written in two equivalent ways, either in 
Heisenberg picture or with Liouville operators as 

{R){t) = TvRpit) = Tre'^(*-*°)i?e-^^(*-*«)p(to) = Tr i-zLj,)e-''^<^'~'°^ psito) Pres , 

where the Liouvillian corresponding to the observable R is defined by 

Lr = ^ [i?, •]+ (69) 

via the anticommutator. If i? = |s)(s'| denotes a Hubbard operator for the eigenstates \s) of 
the quantum system, (j68p is an expression for the matrix element {s\ps(t)\s') of the reduced 
density matrix of the quantum system. 

Diagrammatic expansion. Analogous to the Hamiltonian (fT3|) . the Liouvillian can be 
decomposed into a part for the reservoirs, the quantum system, and the coupling 

L — Lq + Ly , Lq = Lres + Ls (70) 

with 

Lres — [Hres,-]- , Ls — , Lq ~ [Hq, , Ly ~ [V, ■ (71) 

Note that we take a different definition for Ly (containing the commutator) compared to Lji 
(containing the anticommutator and a different prefactor, see (j69|) ). From the context it should 
always be clear whether the coupling V or an observable R is considered. 

Since we are aiming at a diagrammatic representation which is perturbative in the coupling, 
we insert (j70p in (j67[) or (|68p and expand in Ly using time-dependent perturbation theory. 
This leads to a compact notation in Liouville space. Within the usual Keldysh formalism, 
one would insert the Hamiltonian in ([57)1 or (|68p and expand the forward and backward 
propagators e^'^^*"*") in the coupling V. As a consequence, an additional label (the so-called 
Keldysh-index) is needed to distinguish between the forward and the backward propagation 
and the time evolution can be visualized on the Keldysh-contour, see Fig. [3l The notation in 
Liouville space is much more compact, the Keldysh-indices are hidden in the larger dimension 
of Liouville space (which is the square of the dimension of the quantum system) represented by 
the super- matrix notation ([64]) . Alternatively, one can say that the two parts of the Keldysh- 
contour have been taken together in Liouville space, see Fig. [3l However, instead of performing 
the perturbative expansion in time space, it is much easier to do it in Laplace space (at least 
for explicitly time- independent Hamiltonians as we are considering here). First, we define the 
reduced density matrix in Laplace space by 

/>oo 

ps{E) = / dte^^(*-*°)p5(i) , (72) 
Jt„ 

where E is the Laplace variable, having a positive imaginary part to guarantee convergence. 
Using ([BT]). we obtain 

Z 1 

Ps{E) = TVres-^ T p{to) = Tr^es — ■;— Ps{to)pres ■ (73) 

— L tj — Lres — — Ly 
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" '■ • -iLr -iLv -iLy -iLv -iLv P(to) 

P(to) • • m < m • 

iV iV 



t ti t2 ti t4 to t tl t2 ts t4 to 

Fig. 3. Expansion of the time evolution of an arbitrary observable R in the coupling on the Keldysh 
contour (left figure) or in Liouville space (right figure), using the two different ways to write an average 
over an observable R either in Heisenberg picture or with Liouville operators, see Eq. (|68|l . For the 
Keldysh contour, operators are ordered along the contour in the direction indicated. The propagators 
between the vertices contain the unperturbed Hamiltonian Hq. They are given by e-'-^o('"-'^' for 
the upper branch of the Keldysh contour (forward propagation) and e"'"'^'''''"'^' for the lower branch 
(backward propagation), where ti > tj in both cases. In Liouville space all superoperators are time- 
ordered along one axis, the two parts of the Keldysh contour have been taken together. Correspondingly, 
the intermediate state in Liouville space has to be characterized by two states of the quantum system, 
one for the upper and one for the lower part of the Keldysh contour. For convenience, we have chosen 
the time direction to the left since with this convention time-ordering coincides with the sequence 
of Liouville operators (operators with larger time stand to the left of operators with smaller time 
argument). 



This expression can easily be expanded in Ly by a geometric series leading to terms of the form 
i Tvres -t; ^ J— Lv — T-^ — Lv ■■■ Lv — t-^ — ps{h)Pres ■ (74) 

The next step is to integrate out the reservoirs, i.e. the trace Tr^es over the reservoir degrees 
of freedom has to be performed in (|74p . This is achieved by decomposing (|74p into a product 
of a reservoir and system part and applying Wick's theorem to evaluate the average over the 
reservoir distribution. We first exhibit explicitly all parts of the reservoir operators in ([7^ 
by finding a representation of the coupling Ly in Liouville space similiar to the form of the 
couphng V in Hilbert space, given by (|29|) . We use the form 

Lv = \aP^ -P-Gl]J- -.Jl^.-.J^-. . (75) 



n 



Here, Jf are quantum field superoperators in Liouville space for the reservoirs, defined by 
(where A is an arbitrary reservoir operator) 



ai A for p = 
A ai for p 

p is the Keldysh index indicating whether the field operator is acting on the upper or the lower 
part of the Keldysh contour. G\^ "^" is a superoperator acting in Liouville space of the quantum 
system, and is defined by {A is an arbitrary operator of the quantum system) 

^Pi...p„ A - X A f ^ for n even 1 J gi...nA for p = + , . 

^i...n ^ - ''PPI ■■ -"PP" |crP for n odd / \ -A5i...„ for p= - ' ^" ^ 

We implicitly sum overp = ± on the r.h.s. of this definition, and (jPi - P" is a sign-superoperator 
acting in Liouville space of the quantum system, accounting for fcrmionic sign factors. For 
fermions, it is defined by its matrix representation 

lrrPl■■■P^^ - A A ( P2 ■ P4 ■ ■ ■ ■ ior Ns ~ Ns' even , . 

[a )ss'.ss' - OssOs's' . p3 . . . . for - Ns' odd ' ^ '^^ 

whereas, for bosons, it is defined by the unity operator. For n — 1, (|78p has to be interpreted 
as 

1 for iV, — N^i even 



iaP)ss',ss' - SssS,,g, { ^ _ . (79) 
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c~ :■ ' 

J2 

Fig. 4. Determination of the fermionic sign occuring by moving the reservoir field operator J2 from 
the lower to the upper part of the Keldysh contour such that the sequence J2 ... is obtained. 
All Ji belong to the same coupling vertex, s and s' indicate the intermediate states for the quantum 
system on the upper and lower part of the Keldysh contour just left to the coupling vertex. 

We note some important properties of the sign operator, which will be needed later for some 
proofs 

^Pi...p.p'i...p;„ ^ ^Pi...Pn for n even , (80) 

^Pl.-.Pr. Qp^l-Jr. ^p',...p'^ ^ ^p^...p„p{...p'^ QP^l..^^ (g^) 

The proof is quite easy and follows directly from the definition ([75)1 and from the fact that, for 
n odd, the vertex G'i...„ changes the parity of the particle number difference between the states 
on the upper and lower part of the Keldysh contour, and leaves it invariant for n even, i.e. for 
the matrix element (Gi...n)ss' ss' have 

(_1)W.-A^.' ^ . (82) 

The proof of (|75p is a straightforward exercise of algebra and is provided in Appendix B. 
The definition of the vertex operators G^^"^" may at first seem unusual and it deserves some 
comments. It relates to the choice of the sign operators, which is in fact not unique, especially 
the distinction between Ng — Ng' being even or odd in ([78]) is not necessary (as can be seen 
by the proof in Appendix B, where this distinction is not used at all). However, there is a 
special reason why the sign operator (7^i---P" has been defined in such a way, thereby fixing the 
definition of the coupling vertex G\^"^" in Liouville space. The combination 

crPi-P" : Jf ... JP" : (83) 

has the property that the sign-operator exactly compensates additional signs due to interchanges 
of fermionic reservoir field operators, which arise due to the presence of field operators on the 
lower part of the Keldysh contour, i.e. for pi — —. The sign from a^^-'-P" is precisely the sign 
obtained when permuting the later operators to the corresponding position on the upper part 
of the Keldysh contour, assuming that all fermionic reservoir field operators anticommute. This 
means that the determination of the fermionic sign within Wick's theorem (see below) can 
be determined as if all fermionic reservoir field operators lie on the upper part of the contour 
precisely in the sequence . . . and as if the sign-operator were not present. Although there 
is a correction sign from the permutation of field operators belonging to the same contraction 
to be considered (see below), this simplifies the determination of fermionic signs considerably. 
To see this, consider the situation depicted in Fig. Hj where the second field operator J2 is 
intended to be moved from the lower to the upper part of the Keldysh contour along the 
reversed direction (only virtually to determine the corresponding sign from interchanges of 
fermionic operators) . The fermionic sign is determined by the parity of the number of fermionic 
reservoir field operators it passes. Up to the first field operator , the parity is identical to 
the parity of Ng — Ng', where s and s' are the intermediate states of the quantum system on 
the upper and lower part of the Keldysh contour left to the considered coupling vertex. The 
reason for this is the fact that the total parity of fermions (reservoirs plus quantum system) is 
conserved under the coupling and, consequently, also only those matrix elements of the reduced 
density matrix are unequal to zero where the parity of fermions of the quantum system are the 
same (in other words, the external operator R = \s){s'\ in Fig. [3] does not change the parity 
of fermion number). Thus, by moving J2 finally also through Ji, we see that the parity of the 
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total number of interchanges is odd (even) if Ng — Ng' is even (odd) , which is precisely the value 
the sign-operator (j78p produces. The same is obtained if any Jj. with k even is moved virtually 
to the upper part of the contour, whereas for k odd the total parity is identical to the parity 
of Ng — Ns' , again corresponding to the definition ([78|) . The same proof can be used if several 
Jk are moved to the upper part, one just has to move all Jk subsequently to the upper part, 
starting with the smallest k. 

From the definitions, the following useful properties follow directly for the Liouville operators 

Ti-sLs = , (84) 
^ TvsG'^-r=0 , (85) 

Tires Lres ~ , (86) 

where Tr^ denotes the trace with respect to the states of the quantum system. The properties 
and (155)) will turn out to be crucial for the property of conservation of probability. 
Furthermore, from the (anti-)symmetry (j30p of we get 

'~'l...i...j...n — ^'-'l...j...i...n ■ \°') 

The hermiticity of the Hamiltonian Hs — Hg and the corresponding condition ([32| for the 
coupling vertex imply the relations 

{Ls)ss',ss' = ~{Ls)*s's,s's ) (88) 

(r'Pi---P^\ - _ / 1 ioi {n even) or {Ns - N^' even) 1 ,^p„...piN* 

{^i...n )ss',ss' - I -1 for (n odd) and (iV^ - iV^' odd) / " -i ''^''^^^'^ ' ^^-^^ 

where p — —p. This can be shown by some straightforward algebra. The prefactor in the last 
equality stems from the term in the definition (j77p of the coupling vertex G. Using the 
definiton ([79]) of a^, it can also be written in the form 



1 for {n even) or {Ng — Ng' even) 1 
-1 for {n odd) and (TV,. - N^, odd) J 



((a-)"),.^.,- . (90) 



The properties (f88|) and (j89|) can be written more elegantly in operator notation by defining 
the c-transform A"^ for an arbitrary operator A in Liouville space by 

{A'^)ss',ss' = A*s's,s's ; (91) 

which, concerning the Keldysh contour, corresponds to interchanging the states on the upper 
and lower part of the contour and taking the complex conjugate (note that this definition has to 
be distinguished from taking the hermitian conjugate, defined by {A'^)ss',ss' — Agg, ^^,). Using 
this formal notation together with ([TO)) . ([55)1 and ([55]) can be written as 

(LsT = -Ls , (92) 
(Gr.;,f")^ = -(^-rG|:;f . (93) 

Reversing the sequence of all indices, the last property can also be written in the form 

(G?:.;;f")^ = --Gf-/" , (94) 

where n minus signs occur in the superscript of the sign operator. For the proof we used that, 
for n = 2r {n even) or n = 2r + 1 (n odd) , 

{a-r = {±r<7—- (95) 

and (see ([87]) ) 

G^::r = (±)'-G?-;/" . (96) 
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Furthermore, for later purpose, we note the useful relations 

{A By ^ A^B" , (97) 
[Aa^ ^ A-'a) , (98) 

where A, B are superoperators and a is an operator. 

We turn back to the task to perform the trace over the reservoir degrees of freedom in each 
term ([74| of perturbation theory in Ly- We insert the form ([75]) for all Ly and decompose 
the whole expression into a product of a part for the quantum system and the reservoirs by 
successively moving all reservoir field operators Jf' through the resolvents to the right, starting 
from the last Ly - Thereby, we use the property 

^res — (-^res ^l) J\ i (99) 

where we have used the short-hand notation 

Xi = r]^ {uJi + Mai) I (100) 

which will be used frequently in the following (all frequencies occur only in this combination). 
From ([99)1 we get 



: Jf ^ ■■■JP" ■■ = : Jr • • • ^^ : • (101) 

E — Lres — Ls E + Xi ^ . . . + Xn — ^res — Ls 

In this way, all reservoir field operators can be moved to the right. When moving the trace 
Tires in ((74|l to the right, we use the property ([86]) and can set Lres ^ in all resolvents. As 
a consequence, the term in order r + 1 of (I74p can be written symbolically in the product form 

* F ^ T (~T^ V T (~T'^ G)2 . . . ^ — {—.<7 G)r+i ^ /9s(io) 

h — Ls n\ L + Xi — Ls ni h -\- Xr — Ls n\ L — Ls 

•((: J... J:)i(: J... J:)2 ... (: J... J:)r+i)p_ , (102) 

where {^<^G)i indicates the sign and vertex operator at the i-th position and {: J . . . J :)i is 
the corresponding sequence of reservoir field operators. The energies Xi are given by the sum 
of all Xfc-variables from the field operators J^*" which occured left to the i-th resolvent. 

The reservoir part of (|102p can easily be decomposed into product of pair contractions by 
using Wick's theorem. We define the following contraction for the Liouville field operators, 
which can easily be calculated from (P7| 

7ff: = .Jfjl ^ {y} {Jl4)p..s ^ {y} T^resJl4pres = ^iP p' | J } Uvp' U^) , 

(103) 

iTv = 5n' P' { ^ } Pu{^) Urip'io) , (104) 

if the density of states Pi,{uj) is taken into the contraction according to the choice ((33)) and 
(f34| . The upper (lower) value corresponds as usual to bosons (fermions). Note that there is a 
prefactor p' for fermions in the definition of the contraction which arises as follows. As explained 
above the fermionic sign from the sign operators is compensated by moving all ' with Pi — — 
to the upper part of the Keldysh contour. This means that we can use the rule that each 
interchange of two Jf ' gives a minus sign for fermions, independent of the value of the Keldysh 
index pi. However, in doing so, we do not obtain the correct sign from Wick's theorem where it 
is not allowed to permute two field operators belonging to the same contraction. If we consider 
a contraction with p' = we see that we permute the two field operators belonging to this 
contraction when moving Jj7 to the upper part of the contour, see Fig. [5] for an illustration. 



or 
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Fig. 5. Illustration why an additional fermionic sign is obtained when two field operators are contracted 
with the right field operator being on the lower part of the Keldysh contour (left figure). Moving Jii 
from the lower to the upper part along the reversed direction of the contour, it passes through Ji, 
giving rise to an additional fermionic sign compared to the sign from Wick's theorem. In contrast, if 
the right vertex lies on the upper part (right figure), the two field operators of the contraction do not 
pass through each other when moving Ji to the upper part. 



1 2 



3 4 12 



3 4 

Fig. 6. Diagrammatic representations of the two contributions of Eq. (|105|) . A dot with index i repre- 
sents a field operator * . Dots standing close to each other belong to the same vertex G. For clarity 
they are separated to indicate their sequence in the original expression. This is important for the deter- 
mination of the correct fermionic sign. Contractions are indicated by green lines connecting the dots. 
The horizontal black lines between the vertices represent the propagation of the local quantum system 
in Liouville space. 



Therefore, in order to get the correct sign from Wick's theorem, we have to permute these two 
field operators back, leading to an additional sign for fermions. A minus sign is only obtained 

for fermions and p' = leading to the prcfactor | y j" ^-"^ P03p . As a consequence, we use the 

following rules for the Wick decomposition 

1. Contract all J-operators in (|102p such that no contractions occur within the normal-ordered 
parts. 

2. Disentangle the contractions into a product of pair contractions by leaving the sequence of 
J-operators within one contraction invariant. For each interchange of J-operators, give a 
minus sign for fermions. 

3. Translate the contractions with ()103p and sum over all possibilites to contract the field 
operators. 

As an example, we obtain for 

I 1=; — I I — F=; — I 

/. jPl tP2 . . tP3 tP4 .\ _ jPl JP2 tP3 jPi I jPl JP2 tP3 jPi 

\. Ji ^2 ■ ■ ^3 •^4 -/pres — 'Jl ''2 "^3 "^4 ^ '^l '-'2 "^3 "^4 



= ±7fr7fr + 7rr7^r ■ (105) 

The corresponding diagrammatic representation is shown in Fig. [HI 

The energies Xi in (|102p can also be determined by a simple rule. Since the contraction 
P03p is only nonzero for 77 = — ^ = w', and a = a', we get x + x' = according to the 
notation pOOp . Thus, all contractions which either stand completely to the left or to the right of 
the i-th resolvent, do not contribute to the energy Xi. Only the 2:-variables from contractions, 
which cross over the resolvent, contribute to Xi. Thereby, the cc- variable has to be chosen from 
the J-operator standing left to the resolvent. This can be visualized by a simple diagrammatic 
rule, shown in Fig. [71 At the position of the resolvent under consideration, draw an auxiliary 
vertical line and consider all a;- variables of contractions which cross this line (always taking the 
a;-variable from the vertex standing left to the resolvent). The sum of all these a:- variables is 
the energy Xi. 

Finally, we consider the determination of the prefactor, arising from the combinatorical fac- 
tors ^ in p02p . These cancel almost completely with another factor arising from the number 
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Fig. 7. Rule how to calculate the energy variables of the resolvents. The energy X corresponding to 
the resolvent standing between index 3 and 3 is given by X = a;i + a;2 + x^, where Xi, i — 1, 2, 3 are 
the variables corresponding to the contractions which cross the auxiliary vertical red line. Thereby, the 
K-variable of the field operator standing left to the resolvent has to be chosen. Note that the variable 
X4 belonging to the contraction between index 4 and 4 does not contribute because it does not cross 
the vertical cut. 



of identical diagrams when the J-operators within each normal-ordered set are permuted. The 
value of all these diagrams is exactly the same because the two fermionic signs arising from 
interchanging two reservoir field operators from the same vertex and from interchanging the 
two corresponding indices of the coupling vertex G cancel each other, see ()87|) . However, if 
m contractions are present which connect the same normal-ordered blocks, there are m! ways 
of permuting the J-operators of both groups in the same way without giving a new diagram. 
Therefore, in this case, the combinatorical factors can not be omitted completely but a sym- 
metry factor 1/m! remains. 

Summary of diagrammatic rules. We summarize the diagrammatic rules derived in this 
section to calculate the reduced density matrix of the quantum system in Laplace space. The 
value of a diagram is symbolically given by 



psiE) -> ^(±)^^ (n 



S 

— G G...G G — - — psUo) , (106) 

E-Ls E + Xi-Ls E + Xr-Ls E-Ls'^ ' ' ^ ' 

where we use the following rules to calculate the various parts 

1. S — rii<j "^jj- is a symmetry factor, where m.y is the number of contractions between vertex 
i and j. Two diagrams are considered to be different if they can not be mapped on each 
other by permuting only the field operators of each vertex (the field operators are indicated 
by dots in a diagram, where dots standing close to each other belong to the same vertex). 

2. (zb)^'' is a fermionic sign factor, where Np is the number of interchanges of fermionic field 
operators Jf in Liouville space which are needed to write the contractions in product form. 

3-117 stands for the product of all contractions. If Jf and Jf, are contracted, and Jf stands 

left to Jf , the contraction is given by 7ff, , with (see (|103p ') 

7ff: = Vp'I^} M^) , (107) 

where f{uj) = l/(e'^/-^° =F 1) is the Bose (Fermi) function, and Sii, = 5ri.-n'Suv'S{uj — uj'). 

4. To determine the energy argument Xi of resolvent i, we draw an auxiliary vertical cut at 
the position of that resolvent. Xi is the sum of all a;- variables of the contractions which 
cross the vertical cut. The x-variable of a contraction 7ff, is defined as a; = rjluj + fj,a), i.e. 
refers to the left Jf-operator of the contraction. 

5. G = G\^"^"^ are the coupling vertices acting on the quantum system, defined by ([77]) . 
Ls = [Hs, •]- is the Liouville operator of the quantum system. 

6. The formal indices 1 = rjuuj and v = aa . . . contain the index rj — ± for creation/annihilation 
operators, the energy uj of the reservoir state (relative to the chemical potential fia), the 
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>-• — 1 
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1 2 



3 4 5 



6 7 8 9 



10 11 12 



Fig. 8. Example of a diagram with four vertices. The vertical cuts are auxiliary lines between the 
vertices to determine the indices of the resolvents 11,., in (|109p . Precisely those indices occur which 
belong to the left field operator of each contraction crossing the vertical cut. 



1 2 3 4 1 2 3 4 

Fig. 9. Irreducible (left) and reducible (right) diagrams. In the left figure any vertical cut between the 
vertices hits at least one reservoir contraction. In the right figure, the vertical cut between the second 
and third vertex does not hit any contraction and corresponds to a resolvent of the form . 



reservoir index a, the spin index cr, and possible other quantum numbers characterizing the 
reservoir state. We sum (integrate) over all these indices implicitly. 

To write the resolvents in a compact way, we use the short-hand notation 

1 



_/Ti2 n — 

E + Xi + X2 + ■ ■ -Xn - Ls 

With this convention, the diagram example of Fig. [8] is given by 

2 (=•=) 7l7 723 74,12 756 78,11 79,10 

■ Gl2 G345 7Ti45 G6789 ^?489 Gio. 11,12 



(108) 



E-L; 



■PS {to) 



(109) 



where we have omitted for simplicity the obvious Keldysh indices at the couping vertices G 
and the contractions 7. The factor ^ arises from the fact that the third and fourth vertex are 
connected by two contractions, and the sign factor ± stems from taking the contraction 74.12 
out of the rest. 

Kinetic equation. The diagrammatic expansion can be formally resummed and written 
in the form of a kinetic equation by distinguishing between irreducible and reducible diagrams. 
Irreducible diagrams are those diagrams where any vertical cut hits at least one reservoir 
contraction, i.e. in each resolvent at least one x-variable occurs. In contrast, in reducible dia- 
grams there are vertical cuts crossing no contraction, corresponding to a resolvent of the form 
1/{E — Ls) without any a;- variable, see Fig. [9] for illustration. 

We denote the sum over all irreducible diagrams by the irreducible kernel S{E). Using 
(jl06l) . we obtain the following value for any diagram of the kernel 



E{E) ^ i(±)^- 



n 



G- 



1 



E + Xi-Ls 



G ...G- 



1 



E^Xr 



Ls 



G 



(110) 



where (n7)ii-r means that we only consider irreducible diagrams. Compared to l|106p . we have 
omitted in this definition the prefactor z, the first and last resolvent, and the initial density 
matrix psito). 

Each diagram for ps{E) can be written as a sequence of irreducible parts with resolvents 
\/{E ~ Ls) in between. Similiar to Dyson equations within Green's function methods, all 
irreducible diagrams can be formally resummed to define the kernel S{E), and the total sum 
of all diagrams can be written as a geometric series 
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leading to the final result 



where 

L^ffiE) = Ls + S{E) (113) 

is an effective Liouville operator of the quantum system, which depends on the Laplace variable 
E. From (I85|) and the fact that any diagram for S{E) starts with a coupling vertex G, we get 
the important property 

TrsS{E) = TYsL%f^{E) = , (114) 

which is important for conservation of probability (see below). Furthermore, in analogy to (|92p . 
we have 

{L%ff{E)r = -Llff{-E*) , {SiE)r = -E{-E*) . (115) 

The proof of this relation is provided in Appendix C. 

Eq. (|112p is the central result of this section. It shows very clearly the effect of the coupling 
to the reservoirs. In the absence of a coupling to the reservoirs, the kernel S{E) is zero, and 
L'g^ is identical to the bare Liouvillian L5, which is hermitian. As a consequence, the poles of 
ps{E) lie on the real axis, corresponding to coherent Rabi oscillations of the quantum system 
in time space. In contrast, when the coupling to the reservoirs is nonzero, a dissipative part 
E{E) has to be added to the effective Liouvillian and the analytic structure of the reduced 
density matrix changes in Laplace space as illustrated in Fig. 1101 Generically, a branch cut 
will occur on the real axis due to the continuous spectrum of the reservoirs. Analogous to the 
theory of quantum decay, we turn this branch cut into the lower half of the complex plane by 
analytic continuation. This leads to poles in the lower half plane, which originally were on the 
real axis in the absence of the coupling to the reservoirs. These poles correspond to exponential 
decay, the negative imaginary part is the decay rate F (relaxation or dephasing rate, depending 
on whether the mode corresponds to decay of diagonal or nondiagonal matrix elements of the 
reduced density matrix of the quantum system) , and the real part corresponds to the oscillation 
frequency h (e.g. an effective magnetic field). The remaining branch cuts can e.g. lead to power 
law decay, but usually their prefactor is smaller than the one of the exponential decay modes, 
and they dominate only the long-time behaviour. Generically, there will be always a single pole 
at = 0, which corresponds to the stationary state (for certain symmetries or in the case 
of symmetry breaking, it may not be unique accidentally) . It is determined by the eigenvalue 
equation 

L'r{^^^)Ps = . (116) 

This pole will play an essential role within the RG formalism presented in Sec. IH since it does 
not decay. It will be shown that it can be included perturbatively into the initial condition of 
effective vertices, whereas the decay rates Fi and the oscillation frequencies hi provide the cutoff 
of the RG flow. An important advantage of the present formalism is that the physical decay 
rates, determining the time evolution of the reduced density matrix of the quantum system, 
follow directly from the poles of Eq. (|112p . Once the kernel S{E) has been calculated within 
perturbation theory (using the diagrammatic rules ()110|1 ) or even nonperturbatively using the 
RTRG-FS scheme set up in Sec. IH the decay modes can easily be found. In contrast, within 
slave particle formalism, the physical decay rates have to be calculated in a complicated way by 
combining self-energy insertions and vertex corrections [24j , and within flow-equation methods 
the decay rates follow from the energy scale where certain 1-loop and 2-loop contributions 
become of equal order on the r.h.s. of the flow equations. 

Conservation of probability follows from (jll4p . Acting with the trace Tr^ over the quantum 
system on Eq. (|112p . we obtain 

Tvs~Ps{E) = ^Trsps(to) , (117) 
hi 
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Im(E) 



pole at zero 
(stationary state) 



(e.g. power law 
decay) 




Re(E) 



poles 
(exponential decay) 



Fig. 10. Analytic properties of ps{E). Poles occur in the lower half plane leading to exponential decay. 
r denotes the relaxation or dephasing rate, h is the oscillation frequency. The pole at _E = is always 
present and unique, it corresponds to the stationary state. Generically, branch cuts will also occur, 
leading e.g. to power law decay. 



which, after transforming back to time space, gives 

Tvspsit) - Trspsito) , (118) 

i.e. the normalization of the reduced density matrix is invariant and stays unity if it is normal- 
ized initially TTsp{to) = 1. Note that this property holds within any diagrammatic approxima- 
tion, because any diagram for the kernel S{E) starts with a coupling vertex G. We will also 
see in Sec. |4] that the RG flow within RTRG-FS preserves conservation of probability in any 
approximation scheme. 

Furthermore, by applying the property ((98)) to (|112p . we can show from (|115p and psihY = 
Ps{to) that ps{E) fulfils the condition 

psisy = Psi-E*) , (119) 

which is equivalent to the hermiticity of the reduced density matrix in time space 

ps{t)^ = ps{t) . (120) 

Eq. pi2p can also be written in time space leading to a kinetic equation. Multiplying 
Eq. (fTT2l) with -i{E - Ls - S{E)), we obtain 

l-iEpsiE) ~ ps(io)] + iLsPsiE) - -ii:{E)ps{E) . (121) 

Using ([7^ and the definition 

/•oo 

i:(E) = / dte'^'S{t) (122) 



for the kernel in time space, we see that (|12ip is equivalent to the following kinetic equation in 
time space 

-ps[t) + iLspsit) ^ -I j dt' S{t-t')ps{t') . (123) 

The second term on the l.h.s. corresponds to the von Neumann equation of the isolated quantum 
system, whereas the r.h.s. describes the non-Markovian dissipative influence of the coupling to 
the reservoirs. 

Kinetic equations of the form (jl23p are not new in dissipative quantum mechanics and can 
also be derived by other methods, e.g. with projection operators [25] or within slave particle 
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techniques [45j . However, the crucial point is not the form of the kinetic equation (which is trivial 
and obvious on physical grounds) but the way the kernel S is calculated. Using projection 
operator techniques a purely formal but quite compact expression of the kernel is obtained 
where certain projectors Q = 1 — P , with P = presTr^es, occur between the Liouville operators 
Lv projecting on the irreducible part. However, this does not help at all for the calculation of 
E because the reservoir degrees of freedom are still present in Ly- Only after the insertion of 
Q = 1 — P, an explicit calculation can be started but an artificial decomposition into reducible 
parts is created, induced by the projection operator P. All the reducible terms finally cancel in 
a complicated way, leaving only those terms of the Wick decomposition which are irreducible. In 
contrast, the diagrammatic rule (jllOp derived in this section considers directly the irreducible 
terms, and the reservoir degrees of freedom are already integrated out. Using slave-particle 
techniques and Keldysh- formalism, the derivation of a kinetic equation (or quantum Boltzmann 
equation) is quite complicated (even in lowest order perturbation theory), and the kernel is a 
complicated mixture of self-energy contributions and vertex corrections. Therefore, we believe 
that the calculation of the kernel via the diagrammatic rule (|110|) is very efficient. This has been 
demonstrated recently within perturbation theory up to fourth order in the coupling vertex for 
problems of molecular electronics [60]. In Sec.|4]we will see that it is especially useful for setting 
up nonequilibrium RG methods which incorporate the physics of relaxation and dephasing. 



3.2 Observables 

Observables. The time evolution of an arbitrary observable R can be calculated starting from 
Eqs. and ^ 

{R){t) = TTsTTresi-iLR)e-''^^'-*''^ psito)pres , Lr^^[R,-]+ ■ (124) 

The observable R is written in the same generic form ((29|) as the coupling V 

1 

i? = — ri...„ : ai . . .a„ : , (125 
n\ 

with n = 0, 1, 2, . . . (in contrast to V, where the n = case can be incorporated in Hs, this is 
not possible for R). Due to (anti-) symmetry and hermiticity, we get similiar to (|30)) and (132]) 
the properties 

ri...i...j...n = ±ri...j...i...n , (126) 
(ri2...„)^ = rft...i . (127) 

Analogous to ()75|) . we obtain a corresponding form for the operator Lr in Liouville space 
i« = ^a^-^"<;r :Jf ...J^: , (128) 
with the vertex of the observable given by 

nPi...p„ A ^ ix X / 1 for n even 1 /ri...„^ for p = + , . 

^1-" 2 ■ ■ ■ { for n odd j {An,,,,, for p = - ' 

The difference to ([77)) stems from the form Lr = ^[R, •]+ where the anticommutator and a 
different prefactor occurs in comparism to Ly = [V,-]-. The prefactor i/2 is taken into the 
definition of the vertex R to get diagrammatic rules similiar to the ones for the reduced density 
matrix and the kernel S (see below). We note that the trace over the states of the quantum 
system always occurs left to Lr in the average (I124p . i.e. we get 



Trs Lr... ^ ^ (Lr) 

s 



(130) 
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Therefore, according to the definition ((79|) . the sign operator ct^ can be dropped in (|129p for 
the calculation of {R){t), and o-pi - P" can be replaced by 

aP'-P- P2P4... . (131) 

Nevertheless, for the formal identity (|129p the sign operators have to be used in their general 
form. 

Similiar to ([M]) , one can easily show from the definition p29p and the hermiticity condition 
p27)) the relation 

fc'-r = ■ (132) 

In Laplace space we obtain from (|124p 

1 

{R){E) ^ dte'^(*-*«)(i?)(t) = TrsTr.esifl.^^ ^ Ps{ta)Pres , (133) 

Jto ^ ^ ^res — Ls - Lv 

which, after expanding in Lv^ leads to terms of the form 

Tr5 Tr^es Lr — — Ly — — Ly ■ ■ ■ Ly — — Ps{to)Pres ■ 

— I^res — — ^res — ^ — ^res — 

(134) 

The trace over the reservoirs can be evaluated in the same way as described in Sec. l3.1l for ([7^ . 
and we obtain for a certain diagram of the average of an observable in analogy to (|106p 



{R){E) ^ 1(±)^^ (Ht' 



S 

^ E + Xi — Ls E + X2 — Ls E + Xr — Ls E — Lg ^ ^ ^ 

with the essential difference that the first vertex corresponds to the vertex of the observable R 
and the trace Tis over the quantum system has to be performed. 

Decomposing (jl35p into reducible and irreducible parts and resumming formally all diagrams 
leads to 

{R){E) = Trs Sr{E) psito) = ^iTis Sr{E) ~ps{E) , (136) 
E-Lg [E) 

where we have used (|112p in the last equality. Here, the kernel Sfi{E) is defined analogous to 
X!{E) with the only difference that the first vertex is replaced by the observable vertex R, i.e. 
in analogy to (|110p the observable kernel Sii(E) is given by the diagrams 

r„(E)^i(±)".(n-l„«^^:^G,.^Gj^-i-^G ^ ,137) 

Eqs. (|136p and (|137p are the final result for the average of an observable. Once the kernels 
Eii{E) and S{E) have been calculated, they can be inserted into Eq. (|136p and the average 
of an observable can be calculated for all values of the Laplace variable E, i.e. the full time 
evolution can be obtained by transforming to time space. The stationary value of the average 
of an observable is given by 

(Ry* = lim(i?)(t) = lim ElkjiE) , (138) 

which, using ()136p . gives 

{Ry - -tTvs Sr{iO+) pf , (139) 

where p^* is the stationary value of the reduced density matrix of the quantum system, which 
can be calculated from L^^-^ {iO^)pf — 0, see (|116p . 
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Finally, we note that similiar to (jllSp . one can prove 



{UR{E)y = ^Sr{~E*) , (140) 

which, together with the hermiticity condition (|119|) for the reduced density matrix, proves that 
P36p respects the hermiticity of the observable R 



{R){Er = {R){-E*) , (141) 

implying that (R) (t) is real in time space. 

An important example of an observable is the particle current operator flowing from 
reservoir 7 to the quantum system (for electrons, the charge current is obtained by multiplying 
with — e). It is defined by 

I' = -j^N^es , (142) 

where N^^^ is the particle number operator for reservoir 7 and the derivative is calculated via 
the Heisenberg picture for the Hamiltonian 

r = -t[H,N:ij = -t[v,N:ij = 

= -7 gi...n [: ai . . . a„ :, iV^^J 
n\ 

1 " 

— iV" 77, (5q^^5i...„ : ai . . .a„ : , (143) 
1=1 



where the form ([^ for V has been inserted, and we used the identity [ai,N^^g] = —'qSa^ai 
in the last step. Note that we mean by the time derivative -^N^es of the particle number in 
reservoir 7 not the total time derivative, including the one from the coupling of the reservoir 
to the bath maintaining the temperature and the chemical potential fj,-y of the reservoir 
(leading to the grandcanonical distribution for reservoir 7). This total time derivative would 
be zero on average since the average particle number in the reservoir is a constant. Therefore, 
to get a definition of the local current operator at the position where the reservoir is coupled 
to the quantum system, we include in the time derivative only the term i[V,Nres] due to the 
coupling between reservoir and quantum system. 

In summary, the current operator can be brought into the form ()125p 



n - — P 



1 

n! 



ai . 



(144) 



with 

n 

= * n>.5an9l...n ■ (145) 

i=l 

Inserting this form of il „ into ()129p for ri...„, and using ([77]) . we obtain from p28|l for the 
current operator in Liouville space 

Li. = la''-^"^(r)f:/" :Jf ...J^: , (146) 

with 

1 " 

a^)?^:/" = -2 E V^5^n5p.P■■■^.pPG\■:l ■ (147) 

i=l 
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4 Nonequilibrium RG in Liouville space 

In this section we develop a nonequilibrium renormalization group method based on the di- 
agrammatic rules derived in the previous section. It is a formally exact RG-approach in the 
sense that an infinite hierarchy of RG equations will be set up, which, if solved completely, 
would provide the full solution of the problem. Of course , in practice, approximations in form 
of truncation schemes have to be used. Nevertheless, it is quite useful to know what terms have 
been neglected in order to be able to improve the calculations systematically. Formally exact 
RG approaches are used in many fields of physics at the moment and are based on different 
ideas. The most conventional one is to integrate out energy (or momentum) scales step by step 
(starting from high energies) and leaving the whole set of diagrams invariant by renormalizing 
certain parameters, based on the pioneering ideas of Wilson [61 . This basic idea will be also 
the guideline of our RG procedure, although the way we define the cutoff function is quite 
different. Furthermore, our diagrammatic language is more complicated since the vertices are 
still operators on the local quantum system and their time-ordering has to be considered. As a 
consequence, wc have to use the Laplace transform (instead of the usual Fourier transform) and 
renormalization group will generate an additional energy dependence Ls{E) and Gi^"^" (E) of 
the Liouvillian and the vertices, which can only be described within a diagrammatic representa- 
tion (in contrast to approaches based on projection operators). This energy will turn out to be 
the equivalent of the Laplace variable, opening the possibility to address the full time-evolution 
of the problem (however, also for the calculation of the stationary state, the additional energy 
dependence has to be taken into account in the RG equations). We will formulate the RG ap- 
proach such that all vertices stay in normal-ordered form but it can also be formulated without 
normal-ordering. In standard quantum field theory (where all degrees of freedom are integrated 
out and the vertices are c-numbers), similiar RG approaches have been developed based on 
the same idea, see the normal-ordered version by Salmhofer |51j and the non normal-ordered 
version by Polchinski [S^ (usually the RG equations are derived within path integral formalism 
but they can also be obtained on a pure diagrammatic level using the idea of invariance '63]). 
In this sense, the RG method presented in this section is a generalization to nonequilibrium 
and to the case where the vertices are operators. Furthermore, we note that using the Laplace 
transform has the consequence that no energy conservation is associated with the vertices but, 
as we will see, it has the advantage that the structure of the Keldysh indices becomes very 
simple. We will show that by integrating out the symmetric part of the reservoir distribution 
from the very beginning, the final RG equations do not contain Keldysh indices at all. This 
simplifies the calculations considerably. Furthermore, we will see that the physical relaxation 
and dephasing rates (describing the rates associated with quantum transport and not with 
quantum decay or single-particle life times as in usual Green's function methods) occur quite 
natural and can be incorporated directly into the RG equations. Therefore, our RG schemes 
seems to be an appropriate way to combine RG with the physics of relaxation and dephasing, 
an important subject in mesoscopics, quantum information processing and cold atom systems. 

We mention that there are other formally exact RG approaches in quantum field theory 
which do not use the conventional idea of integrating out energy degrees of freedom and leaving 
the sum of all diagrams invariant by renormalization. One of them is the RG developed by 
Wetterich [50j where cutoff dependent 1-particle irreducible vertex function are defined, such 
that the initial condition at cutoff yl = oo gives the bare vertices whereas for A = the full 
physical vertex functions are obtained. Differential equations are then set up directly for these 
vertex functions from which they can be calculated systematically using appropriate truncation 
schemes. This scheme is specifically adapted to the conventional diagrammatic representation 
of quantum field theories and, therefore, can not be overtaken directly to the diagrammatic 
language used in this paper. The concept of irreducibility is defined completely different here 
and the diagrams generated by our RG procedure are automatically irreducible. Within the 
RTRG-FS approach the irreducible diagrams, represented by the kernel S{E), are fed back 
into the RG equations for the vertices via an effective description of the dynamics of the local 
quantum system. This is analogous to the Wetterich RG scheme where the propagators include 
self-energy insertions. For problems where a controlled expansion in the Coulomb interaction is 
possible, the Wetterich RG method is very powerful and, recently, has also been generalized to 
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the nonequilibrium case |53|5|52j . Another RG method is the flow equation method developed 
by Glazek, Wegner, and Wilson |49l , which recently has been generalized to the nonequilibrium 
and time-dependent case by Kehrein et al. [48 58J. This scheme is based on a pure Hamiltonian 
level and the Hamiltonian is diagonalized by unitary transformations step by step. This idea is 
competely different from all other formally exact RG methods and can be viewed as a technical 
alternative to the scheme presented here. 



4.1 Basic ideas and cutoff function 

The central quantities of interest are the effective Liouvillian L'^J^ {E) (or equivalently the 
irreducible kernel S{E)) and the kernel Sji{E), given by the diagrammatic representation (see 
(fTTOll and (fml) ) 



^RiE) j \Rn=a} ' \^^' I itt\R] E + Xx-Ls "' E + Xr - Ls 

(148) 

The first term represents the n = case and is written explicitly because Gn=o is by convention 
incorporated in Ls- The second term represents all diagrams with more than one vertex or, 
equivalently, with at least one reservoir contraction. There, the difference is only the first vertex, 
which is G for Vg^ {E) and R for Sfi{E). From these quantities the reduced density matrix 
and the average of the observable R can be calculated by using pi2p and (|136p 

Ps{E) = /e//,^, /^g(M ' {R){E) = -tTTsSR{E)ps{E) , (149) 

E - Lg (E) 

and the irreducible kernel S{E) can be obtained from 

L^/^iE) ^ Ls + S{E) . (150) 

To find a mathematically well-defined formulation of renormalization group (RG), we note 
that the physical quantities of interest are a functional of the reservoir contraction 7, the vertices 
G and i?, and the Liouvillian Ls 

L%ff{E) = Ls + T{^,Ls,G) , (151) 
ERiE) = i?„=o + Tr{-i,Ls,G,R) . (152) 

The functionals J- and T-n represent the diagrammatic rules of the second term of Eq. (|148l) . 
Note that the n — case is written explicitly in the first term, so that the arguments of the 
functionals do not include the case n = for the vertices Gi...„ and Ri...n- 

The idea of RG is to integrate out the reservoir degrees of freedom and to account for these 
by a renormalization of the system parameters without changing the diagrammatic rules. This 
leads to an effective theory for the dynamics of the local quantum system. To achieve this, 
we replace the reservoir correlation function 7 in the functionals ()151|) and (I152|) by a cutoff 
dependent contraction 7^, and try to find a corresponding yl-dependence of the Liouvillian 
Ls ^ Lg and the vertices G,R ^ G'^,R'^, such that, without changing the diagrammatic 
rules (i.e. the functionals), the sum of all diagrams stays invariant, which gives 

Ll"iE) = Li{E) + Hl\Li,G^) , (153) 
SniE) = RtoiE) + TRi^^'^Li^G^^R^) . (154) 

This is the important property of invariance since the r.h.s. gives always the same independent 
of A. In this way an effective theory can be formulated where both the reservoirs and the 
quantum system are modified, but the physics of the total system stays invariant (in principle 
also the functionals and J-'-jz can be replaced by a cutoff-dependent functional but we will not 
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consider this case here). However, as we will see in Sec. 14.21 invariance can only be maintained 
if the Liouvillian and the vertices get an additional energy dependence 

Ls ^ LiiE) , (155) 
CP..^P„ - (G^)?.V./"(i^) , (156) 
^P..^P„ - {RX'::n"'iE) , (157) 

and a new diagrammatic rule has to be set up which states that the energy argument E + Xi 
occuring in some resolvent of (jl48p is the same for the Liouvillian occuring in that resolvent 
and for the vertex to the right of that resolvent (the first vertex gets the energy argument E) , 
i.e. (|148p has to be replaced by 



^ ■ G'^iE + X,) ... 1— — G^(^ + Xr 



E + Xi-mE + Xi) ' ' E + Xr-L^iE + X. 

The advantage of this formally exact scheme is that one can choose the value for the pa- 
rameter A and the way one defines the yl-dependence of the reservoir contraction 7^ in an 
arbitrary way, opening up many possibilities for RG schemes. However, only those schemes are 
of course of practical use, where the yl-dependence of the Liouvillian and the vertices can be 
found in a systematic way, and where the final evaluation of (jl58p for a certain value of A has 
some advantage compared to the original series (|148p . We discuss some of these schemes in the 
following. 

Continuous RG scheme. If the parameter yl is a continuous parameter with initial value 
Ain and final value Afi, one defines the boundary conditions of the contraction such that 

7^- = 7 , 7^^' = , (159) 

i.e. the initial contraction is the original reservoir correlation function and the final one is 
defined as zero. As a consequence the initial values of the Liouvillian and the vertices are the 
original ones 

(160) 
(161) 
(162) 

such that the original perturbative expansion (|148p is reproduced. For the final value A = Afi, 
the reservoir contraction is zero, and only the diagrams of (I158P survive where no contractions 
are present, i.e. the n = vertices, represented by Lg{E) and R^^q{E). Thus, the result for 
the physical quantities is given by 

L^^^iE) ^ Li^'iE) , Sj^iE) ^ Rtl:oiE) . (163) 

We see that Lg{E) and R^^q{E) fiow finally into the physical quantities and, therefore, we 
interpret them as effective physical quantities at scale A. If A is an energy cutting off the high- 
energy scales of the reservoirs (see (jl69p below), one can interpret the yl-dependent physical 
quantities as containing all energy scales between Am and A. 

Using this picture, one can define a yl-dependent reduced density matrix and a yl-dependent 
average of an observable R by equations similiar to (|149p via 

pME) = E_lA^E^ Psih) , {R)''{E) = -iT:TsSi{E)pi{E) , (164) 

with 

Si{E) = i?lo(^^) • (165) 







'{E) 


= Ls , 




I \Pl---Pr. 

)l...n 




— '-^l...n 




)l...n 
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Ps{E) and {R)^{E) start with the value for the isolated quantum system at yl = Ain and flow 
into the full solution Sit A — Afi. 

In summary, the central task is to find the vl-dependence of the Liouvillian and the vertices 
with the initial conditions given by (|160[) - (|162[) . To express them via differential equations, the 
so-called RG equations, one takes the derivative of the invariance properties (|153p and (|154p 
with respect to A. The l.h.s. gives zero since the physical kernels do not depend on A. On 
the r.h.s. one contribution contains the derivative of the contraction and the other ones the 
derivative of the Liouvillian and the vertices. All these contributions have to cancel each other 
to fulfil invariance and it is a technical task to find this cancellation on a diagrammatic level 
using the representation (I158p . We will discuss in Sec. 14.21 how to achieve this in an efficient 
way. As a result we obtain RG equations symbolically of the form 

.Fr(7'^,L^,G^) , (166) 

.F^«(7^L^,G^) , (167) 

T^^{^^,Li,G^,R^) . (168) 

The minus sign on the l.h.s. indicates that the sum of all terms has to vanish to fulfil invari- 
ance. Thus, the RG functionals JF^^ contain those terms where one derivative of a reservoir 
contraction occurs. To solve the RG equations with the initial conditions (|160p - (|162p . one usu- 
ally has to approximate the RG functionals in a systematic way by expanding in the coupling 
parameter, the so-called perturbative RG scheme. However, it is important to note that one 
expands not in the original coupling parameter, but in the coupling parameter of the renormal- 
ized vertices at scale A occuring on the r.h.s. of the RG equations. This can improve the result 
considerably compared to bare perturbation theory and, therefore, perturbative RG methods 
are a very important technical tool to avoid divergencies occuring often in bare higher-order 
perturbation theory. 

Choice of cutoff function. An important question within the continuous scheme described 
before is of course the most appropriate choice for the yl-dependence of the reservoir contraction, 
which, so far, has not been specified at all besides the boundary conditions (|159p . The criterion 
for the right choice is to achieve that the r.h.s. of the RG equations is a well-defined series in 
the renormalizcd coupling constant so that a perturbative RG scheme is possible. Alternatively, 
one can also stop the solution of the RG equations at a certain value of A and calculate the 
kernels from perturbation theory in the couplings at this value of A using (jlSSp . Whether this 
perturbation theory is well-defined or not depends on the size of the renormalized coupling 
constants (coming out of the solution of the RG equations and can be smaller or larger than 
the original couplings), and the low and high frequency behaviour which can induce divergencies 
for the frequency integrals on the r.h.s. of the RG equations. Divergencies at high frequencies 
are under control by choosing yl as a high-frequency cutoff in the contraction, e.g. the most 
obvious choice would be 

Y^'r"" = 7rrO{A'H) , (169) 

meaning that at cutoff scale A only reservoir energies below A are considered. Within this 
scheme, the initial cutoff is infinity Ain = oo, and the final one zero Afi = to fulfil the boundary 
conditions (|159[) for the contractions. This is the conventional bandwidth cutoff of poor man 
scaling approaches. All frequency integrals in (|158p are cut off by A for high frequencies. For low 
frequencies, there can be further divergencies induced by the sharp step of the Fermi functions or 
the strong increase of the Bose function. These divergencies are cut off by some low energy scale 
Ac provided by the smearing of the distribution functions by temperature or by the other energy 
scales appearing in the denominator of the resolvents in pSSp . like the chemical potentials of 
the reservoirs, the eigenvalues of Lg{E + Xi), or the Laplace variable E. When A approaches 
Ac, the high and low frequency cutoffs are the same and divergencies are absent. This is the 
idea of perturbative RG which works provided that the renormalized coupling constants are 
still small at scale Ac. The latter condition is called the weak-coupling regime which can be 
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solved systematically by perturbative RG. When the coupling constants are already of order 
one at scale Ac, a strong coupling problem occurs and the perturbative RG scheme becomes 
uncontrolled. Truncating the series on the r.h.s. of the RG equations is not justified and it is 
not clear which approximation scheme can be trusted. Nevertheless, it has turned out that even 
for coupling constants of the order of one, only a few terms of the series have to be taken into 
account to get reliable results. Whether this is a generic feature or just accidental for certain 
models is one of the most interesting questions in modern renormalization group theory. We 
will address this issue in Sec. 15.31 for the strong-coupling regime of the Kondo problem. 

As we have seen, an important issue is whether a low-energy scale Ac is present which cuts 
off small frequencies in (|158p (an equivalent discussion can be performed for the r.h.s. of the RG 
equations, see Sec. 14.31) . When temperature is small and when the real part of all low-energy 
scales in the denominator of the resolvents cancel each other (which happens when the sum of 
the real part of the Laplace variable, the chemical potentials, and the real part of the eigenvalue 
of Lg{E + Xi) becomes small), the only low-energy cutoff is provided by the imaginary part 
r of the eigenvalues of the Liouvillian Lg{E + Xi). Therefore, if we denote by Tk the energy 
scale where the coupling constants become of order one (the Kondo temperature for the Kondo 
model), we expect that a weak-coupling problem occurs if the minimum of all relaxation and 
dephasing rates is much larger than Tr- 



However, to justify this statement technically, it is important to prove that the zero eigenvalue 
of the Liouvillian (which generically occurs and corresponds to the stationary state, see (|116|) ) 
does not lead to an accidental situation where no low-energy scale occurs. To analyse this 
problem in more detail, we discuss first the properties of the eigenvector with eigenvalue zero. 
Except for maybe some special values of A and E, the Liouvillian can be diagonalized according 
to the eigenvalue problem for the right and left eigenvectors in Liouville space 



The right and left eigenvectors \x^{E)) and \x^{E)) will in general not coincide since the 
renormalized Liouvillian contains the dissipative influence of the reservoirs and, therefore, will 
be non-hermitian. The eigenvectors are operators with matrix representation 



where \ss') is the Dirac notation for the basis operators |s)(s'| in Liouville space. The normal- 
ization and completeness relations of the eigenvectors read 



min {Fi} 3> Tk weak coupling 



(170) 



LUE)\xt{E)) = \t{E)\xt{E)) 
{xt{E)\Li{E) = Xt{E){xUE)\ 



(171) 
(172) 



{ss'lxtiE)) = {xi{E)) 



SS' 



(173) 



{xt{E)\xf{E)) = 



(174) 
(175) 




fc 



The eigenvalues 



(E) consist of a real and an imaginary part, denoted by 



(176) 



Using (|175p . the reduced density matrix p64p can be written as 




(177) 



and the poles Zj 



of the reduced density matrix follow from the self-consistent equation 



4 = ^t{zt) = ht{zt)-ir^{zt) 



(178) 
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provided that the energy-dependence of the eigenvectors does not induce additional poles. The 
real and imaginary parts of determine the oscillation frequencies and the relaxation/dephasing 
rates at scale A. 

As shown later in Sec. 14.31 the effective Liouville operator at scale A fulfils the same prop- 
erties as the full effective Liouvillian, see (|114l) and (jllSp 

TTsLjiE)=0 , L^iEY = -Lji-E*) . (179) 

Applying the second condition to the eigenvalue equations (|171|) and (|172p . we obtain 

x^iEy = xti^E*) , x^iEy = x^i-E*) , X^iEY ^ -Xti~E*) . (180) 

If we act with the trace Tts over the quantum system on the eigenvalue equations and use the 
first property of (jl79p . we get the relation 

X^iE) J2 4{EU = , (181) 

s 

implying that either the eigenvalue or the sum over all diagonal elements of the eigenvector 
must be zero 

X^{E) = or ^ xtiE),, = . (182) 

s 

Due to the completeness relation, not all eigenvectors can fulfil the property that the sum over 
the diagonal elements is zero. Therefore there exists at least one eigenvector with zero eigen- 
value, where the sum over the diagonal elements is unequal to zero. We denote this eigenvector 
by = and assume that it is unique since it corresponds to the stationary state at scale A. 
In this case, also the left eigenvector with zero eigenvalue is unique and we get 

mE)\x^iE))^0 , (x^iE)\LjiE) =0 . (183) 

The right eigenvector \xq{E)) will change as function of A and E, and will also dependent on 
the model under consideration because each system has its own stationary state. However, the 
important point is that the left eigenvector \xq{E)) is always the same for all A and E, and 
independent of the model, due to the first property of (|179|) 

TrsL^iE) ^Y.^siE)ss,.. = (184) 

This is equivalent to the eigenvalue equation (|183p for Xq{E) if we choose 

x^{E),s' - S,,, , (185) 

up to a normalization constant, which we have chosen to be unity to fulfil the normalization 
condition (fTTi]) for A; = 

{x^iE)\x^iE)) = J2 "^oiEhs = 1 . (186) 

Thus, acting with the eigenvector Xq{E) from the left is equivalent to acting with the trace 
over the quantum system 

{x^{E)\G ^TvsG , (187) 

where G is an arbitrary superoperator. This relation is one of the most important properties 
which is preserved under the RG flow since (|179p is fulfilled. 

We now turn back to the question how to choose the yl-dependence of the contraction such 
that the zero eigenvalue of Lg{E) can not lead to a problem at low energies. Inserting the 
completeness relation ()175p between the vertices of the diagrammatic expression ()158|) we see 
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that the left eigenvector acts on the right vertex and the contribution from fc = leads precisely 
to a term of the form (I187|) 

{4{E)\{GX':J"{E + X,) = Trs{G^)l]-J-{E + X,) , (188) 

where the trace over the quantum system acts on the vertex. We now remind ourselves of the 
property ([55)1 of the vertex which is also preserved under the RG flow (see Sec. 14.31 for the 
proof) 

Trs{G^)l]J-{E)=Q . (189) 

It means that if the average over the Keldysh indices were taken in (|188p . the contribution of 
the zero eigenvalue would be exactly zero! Therefore, the task is to choose the yl-dependence 
of the contraction in such a way that for small values of A only the average over the Keldysh 
indices of each vertex occurs in the perturbation series (|158|) . 

The Keldysh indices occur in ()158|) via the vertices and via the reservoir contractions given 
initially by (fTUi)) 

(190) 



7n' ='5iI'P'|'/| Pu{^) faivp'^) , 



where we have taken the form where the density of states is incorporated. If we were able to 
avoid the dependence of the contractions on the Keldysh indices, only the average of the vertices 
over the Keldysh indices would occur in the perturbation series. The contraction depends only 
on the second Keldysh index p' which appears as a prefactor and in the argument of the 
distribution function. Therefore, only the symmetic part of the distribution function leads to 
a dependence on the Keldysh indices. We now see trivially what we have to do to avoid the 
occurence of the zero eigenvalue. We have to replace the distribution function fa{^) by a A- 
dependent distribution function f^{uj) such that for small A of the order of the relaxation and 
dephasing rates, we get an antisymmetric distribution function 

f^{-uj) = -f^icj) for A ^ min {F,} . (191) 

We conclude that not only the density of states should depend on A but also the distribution 
function, leading to a choice of the form 

pp'A _ r _ ^' J ^7 



= S,,,p'yi^ p^{u;)f^{r^p'u;) . (192) 

In summary, we have now two requirements to choose the cutoff function to avoid problems 
at high and small frequencies in the renormalized perturbation series (|158|) . First we need that 
only energy scales smaller than A should contribute to the frequency integrals, demanding a 
cutoff scheme similiar to the bandwidth cutoff proposed in (|169p . This corresponds to a A- 
dependent density of states in (|192|) given e.g. by 

p^icj) = p,iij)e{A-\ij\) , (193) 

where 9{uj) can be a sharp step function or, if needed, a certain smearing procedure can be 
introduced to get simple analytic properties or to improve numerical stability. Secondly, we de- 
mand the distribution function to be antisymmetric for small values of A to avoid the occurence 
of the zero eigenvalue leading to problems at small frequencies in the perturbative series. There 
are several possibilities to achieve the latter property which we will discuss in the following. 

Integration by discrete RG step. The easiest way to get rid of the symmetric part of 
the distribution function is to integrate it out in one step before the continuous RG starts. This 
means that we decompose / into a symmetric and an antisymmetric part 



1 

^2 



(194) 
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The symmetric part is a constant given by =pi. The corresponding decomposition of the con- 
traction ()190p reads 

iTv = + 111' , (195) 



with 



(196) 



As expected, we see that only the symmetric part depends on the Keldysh indices. 

We replace now the contraction 7 in the functionals ()151|) and ()152p by the antisymmetric 
part and renormalize at the same time the Liouvillian Lg ^ Lg and the vertices G ^ G"", 
i? — > i?° in such a way that invariance holds, i.e. we get in analogy to (|158p the new perturbation 
series 

^ .G%E + X,) \———&[E + Xr) , 



E + Xi - L%{E + Xi) ' ' "' E + Xr-L'i.iE + Xr 

where only the antisymmetric part of the distribution function occurs which is independent of 
the Keldysh indices. Therefore, only the average of the vertices over the Keldysh indices occurs 
which is denoted by G and R 

Gl...n{E) = E Gr.;„^"(^) ' = E <:n'"(^) ■ (198) 

Pl...p„ Pl---P„ 

In this way the symmetric part of the distribution function has been integrated out and has 
been shifted into the initial condition of the Liouvillian and the vertices. As we will describe 
in the following sections, the quantities Lg, and i?" can be calculated from a well-defined 
perturbation theory because the symmetric part of the distribution function is a frequency- 
independent constant, leading to no divergencies in the perturbation series. 

After having integrated out the symmetric part of the distribution function, the zero eigen- 
value of the Liouvillian can no longer occur in the resolvents and we can subsequently apply 
a well-defined continuous RG scheme to integrate out the energy scales of the antisymmetric 
part. This can be done e.g. by the bandwidth cutoff leading to a yl-dependent antisymmetric 
contraction of the form 



iti' = <5ii' { \ \ P^A^) 



(199) 



with p^(w) given by (|193p . However, as we will see in the following, also this choice of a real- 
frequency cutoff is not the most convenient one. Instead, one should again choose a yl-dependent 
antisymmetric distribution function instead of a yl-dependent density of states. 

Cutoff function on the imaginary frequency axis. The real- frequency bandwidth cutoff 
scheme to avoid divergencies for high energies has two important problems. As we will see in 
Sec. 14.21 the r.h.s. of the RG equations (|166p - (|168p has a similiar structure to the perturbation 
series (I158P but in all resolvents one of the frequencies is replaced by the cutoff A. This has 
the consequence that if A crosses the rest of the real part of the denominator of the resolvents, 
the real part of the resolvent rapidly changes sign on the scale of the relaxation/dephasing 
rate F and becomes quite large ~ -p in the vicinity of this point. Although these negative and 
positive contributions almost cancel each other (like for a principle value integral) numerical 
problems occur, leading to low accuracy and slow algorithms, see Refs. |53I64I52] for more 
details. Another problem is the fact, that in each RG step A ^ A — dA, an infinitesimal energy 
shell A > \u!\ > yl — dyl is integrated out in the reservoirs. For A large compared to all other 
physical energy scales, this means that no energy conserving processes are integrated out. As 
a consequence, the generation of relaxation and dephasing rates happens only for values of A 
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Fig. 11. Theta function smeared by temperature. 

of the order of some physical energy scale like temperature or voltage, see Ref. [4]. Therefore, 
for large A, the imaginary part in the denominators of the resolvents is very small and the 
resolvent becomes very large (or can even diverge) when A cancels the other real parts of the 
denominator (note that combinations of real frequencies occur in the denominator which can 
become of the same order as A). Of course, this is not a real problem because the different 
signs on both sides of the divergence cancel each other but a stable numerical solution is very 
problematic. 

To avoid these problems another cutoff scheme on the imaginary axis has been proposed for 
a nonequilibrium system ^ . Instead of integrating out the real frequencies of the contraction 
(|190p step by step, one analyses the analytic structure of the frequency dependence and tries 
to integrate out the poles in the complex plane step by step. For a broad density of states 
depending weakly on frequency of the form (f47|) 

the poles of p{uj) occur at ±iD and have a very large imaginary part. Therefore, we disregard 
them in the following, and analyse only the analytic properties of the distribution function (for 
more complicated problems also the analytic structure of the density of states has to be taken 
into account, see the discussion at the end of Sec. (HI). 

Using the exact representation of in terms of the Matsubara frequencies oj" = 2mrTa 
(w" — [2n + l)7rrQ) for bosons (fermions) 

Uu:) = ±tY,^^ , (201) 



with e = 0+, we see that the distribution function has poles at the Matsubara frequencies with 
residuum ztT. The convergence factor e"^" determines the symmetric part of the distribution 
function, for e = 0^ we obtain fa{uj), whereas for e = 0^ we get — /„(— tj) ~ ±(1 ± fa{uj)). 

To integrate out the Matsubara poles step by step, we introduce a yl-dcpendcnt distribution 
function of the form 

f^(,,) ^ ±tJ2^^(^tAA-K\) , (202) 

n " 

where 

. . ^ / for |.;| > 

^^('^) - I i + ^ for < (203) 

is a theta function smeared by temperature, sketched in Fig. [TT] Via (I192|) . this leads to the 
following yl-dependent contraction 



XA-K\) , (204) 



with X = r]{Lo + fia) according to ()100|1 . and 

Mq = ??Aia ■ (205) 
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If convenient, one can also add a yl-dependent density of states, but we do not consider this 
case here, assuming that the density of states has a well-defined analytic structure, as e.g. given 
by (EOQl). 

With the choice (|202p the parameter A cuts off the Matsubara poles of the distribution 
function. Therefore, as we will see in Sec. 14.21 the complex parameter iA will occur in the 
denominator of all resolvents on the r.h.s of the RG equations and all other imaginary parts 
are also positive. For this reason the resolvents will stay small and the numerical problems 
that occur if one uses the real-frequency cutoff can be avoided. This can also be seen by 
rewriting the original perturbation series (|148|) directly in Matsubara space before starting the 
RG procedure. If we assume that all frequencies are allowed (i.e. the bosonic case with particle 
number conservation is not considered here) and the frequency dependence of the density of 
states and of the vertices is weak, we can close all integrations over the cc-variables in the upper 
half of the complex plane and see from p04p (with yl = oo, i.e. leaving out the 9t^ function 
initially) that only the poles a,t x — jla + ioj" with positive Matsubara frequency > have 
to be considered. Note that all energies Xi in the resolvents of the perturbation series ()148|) are 
sums of a;-variables, so that all of them have positive imaginary parts when the integrations are 
closed in the upper half of the complex plane. Furthermore, also the Laplace variable E has 
positive imaginary part. Performing all integrations over the a;- variables in this way, we get a 
new series of the same form as (|148p but with the replacements 

X ^ fla + ILU" , IJLU iuj"^ , (206) 
7fr' -^11' P.(iry02«re'<^ , (207) 

and we sum only over positive Matsubara frequencies. Starting from this series, we can now 
introduce the same ideas of renormalization group as before and introduce a cutoff function 
0r^(yl — |w"|) into the contraction (|207p . We obtain a renormalized perturbation series of the 
form (|158p with the same replacements (|206p and (|207p . A similiar series occurs on the r.h.s. of 
the RG equations, see Sec. 14.31 Assuming that the imaginary parts of the eigenvalues of Lg{E) 
are all negative (corresponding to positive relaxation and dephasing rates, see (|176p '). we see 
that the imaginary part of all terms in the denominator of the resolvents are positive so that no 
cancellations can occur. As a consequence, this form is very stable for numerical calculations. 

We now show that also all other requirements for a suitable yl-dependent contraction are 
fulfilled by choosing the Matsubara cutoff function (|202p . This concerns the requirement l|19ip 
of antisymmetry for sufficiently small yl, the behaviour for large frequencies, and the question at 
what value of A relaxation and dephasing rates are generated. To discuss this, we first analyse 
the analytic form of the /l-dependent distribution function. 

For zero temperature, f^{oj) is depicted in Fig. [12] for the fermionic case (the bosonic case 
is obtained from changing the sign, disregarding the divergence occuring for |a;| smaller than 
temperature). Analytically we get 



f^{u;) = ±^ I doj'—— , (208) 



2n 

which can be calculated in two limiting cases 

« ^, - ^ f^{^) = T - Si(yle) ± i sign(a;) , (209) 

e TT z 

\uj\,A<$:- ^ /-^(w) = ± - arctan(yl/cj) , (210) 
e TT 

with Si(a;) = Jq dy '^'"J^-' . Using Si(cx)) — ^, we see that for \uj\ <C A, i, the distribution function 

is just moved up (down) by ^ for bosons (fermions) when yl crosses i, leaving the antisymmetric 
part 

I I I 1 

f'^(Lu)=±— / duj' ± - arctan(yl/w) for|w|,y4<- . (211) 

2tt J_y^ u) — iuj' 77 e 
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Fig. 12. Sketch of the yl-dependence of the Fermi function for zero temperature. D is the bandwidth 
setting the scale which frequencies are important. Since e is arbitrarily small, we take i ^ D and three 
limiting cases can be considered for A. For A above i, we have the original Fermi function which is 
a step function. For A between i and D the Fermi function has been moved down by i, i.e. only the 
antisymmetric part of the original Fermi function remains. For A below D the high-energy scales get a 
smaller weight and the Fermi function is given by — i arctan(^). 



This is the form for e 0+ and A finite. The regime e — > 0+ can also be calculated analytically 
for finite temperature. We obtain 

hm, f^iiu) = ±r V -i— OtAA - Kl) (212) 



for bosons 
for fermions 



where At^ is the value of the Matsubara frequency which lies closest to A, and "^(z) denotes 
the Digamma function with asymptotic values 

V(z) ^ ln(z) - ^ + 0{\) for |z| -> oo , (213) 

il){z) ^ -7 - i + 0{z) for |z| ^ , (214) 
z 

where 7 is Euler's constant. 

Thus, we see that it is important to take e finite because otherwise the symmetric part of 
the distribution function is missing at any finite A, or, in other words, the symmetric part is 
integrated out when A crosses \. Since e can be taken arbitrarily small, this means that for A 
reaching any physical scale, the distribution function is antisymmetric to any desired accuracy, 
i.e. the requirement p9ip is fulfilled. 

Instead of integrating out the symmetric part smoothly using a finite value of e, it is more 
convenient for analytical calculations to integrate it out in one single step as described above 
(see (|197p ) and, subsequently, use the choice (|212p for the vl-dependence of the antisymmetric 
part of the distribution function. In this case, the parameter e does not occur at all which 
simplifies the calculation. Thus, for the continuous RG fiow, we set e in (|204p and get 

T^o ^ r 1 _ 1 '■^ 



^-v \ P.Htt / doj' -_ — (216) 

for the yl-dependent contraction, and the perturbation series (|197p turns into the renormalized 
series 

1 I -n^ , V \ 1 riA i 



G'\E + Xi) . . . -r— — G^(£; + Xr 



E + Xi-mE + Xi) ' ' E + Xr - L-^iE + X. 
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where all vertices are now averaged over the Keldysh indices, as defined in (jl98p . and the 
zero eigenvalue of the Liouvillian can no longer occur. This 2-stage procedure of RG, i.e. first 
integrating out the symmetric part of the distribution function and then choosing the Matsubara 
cutoff for the antisymmetric part, is the procedure which is most suitable and will be used in 
the following sections. 

Concerning the behaviour at high energies we see from Fig. [12] that the yl-dependent dis- 
tribution function always contains all energy scales and falls off at high frequencies like Aju). 
This means that the scale where high frequencies become less important is set by A and, thus, 
the divergencies of the original perturbation theory at high energies are eliminated during the 
RG flow. Furthermore, in each infinitesimal step A A — dA of the RG procedure, reservoir 
energies on all scales are integrated out, leading to the effect that relaxation and dephasing 
rates are generated from the very beginning of the RG flow (even the discrete RG step at the 
beginning integrating out the symmetric part of the distribution function leads to an initial 
condition for the rates). This also helps to improve the stability of the flow. 

In summary, we conclude that the Matsubara cutoff is a very suitable choice of a cutoff func- 
tion for fermionic problems with a flat density of states. However, when complicated frequency 
dependencies of the density of states appear, one has to consider the analytic structure of the 
density of states as well and other cutoff schemes might be more appropriate. We summarize 
the criteria for a suitable yl-dependent contraction, which we have found in this section 

1. High frequencies should become less important during the RG flow. This means that the 
contraction should suppress frequencies above A. 

2. The distribution function should become antisymmetric during the RG flow to avoid the 
appearance of the zero eigenvalue of the Liouvillian corresponding to the stationary state. 
Closely connected with this property is the fact that the effective vertices can be averaged 
over the Keldysh indices. 

3. The resolvents should not become large during the RG flow leading to numerical problems. 

4. Relaxation and dephasing rates should be generated already at the beginning of the RG 
flow. 



4.2 Derivation of RG equations 

Single RG step. In this section we describe in detail of how the renormalized Liouvillian 
and the renormalized vertices have to be defined in order to fulfil the central properties of 
invariance, given by (|153p and (|154p . together with the renormalized perturbation series as 
shown in Eq. (I158|) . Instead of directly aiming at the derivation of the continuous RG equations 
(|166p - (|168p . we consider first a single discrete RG step, where a certain finite part of the 
reservoir contraction is integrated out in one step. We consider an arbitrary decompostion of 
the contraction in two parts 

7 = 7^ + 7^ , (218) 

and want to integrate out the 7'^-part, i.e. the aim is to replace 7 7^ in the diagrams and 
define a renormalized Liouvillian Lg{E) and renormalized vertices G^{E) and R^{E) such 
that (fT58l) holds 



Sn{E) / - \ R^^E) S'^' Vii ' y \ R^{E) / 

^ -G^iE + X,) ...^^ l^—-^,G^{E + Xr) 



E + Xi- L§{E + Xi) ' ' E + Xr- L^iE + Xr) 

An example of the decomposition (|218p is the splitting (|195p of the contraction into a symmetric 
and an antisymmetric part, with the task to integrate out the symmetric part in one single 
step. Furthermore, we can also view the continuous RG as a sequence of infinitesimal steps 
A A~ dA, where each step can be brought into the form (|218p by writing 



(220) 
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Fig. 13. Example of an y4-irreducible block. A- and B-contractions are distinguished by green and 
red lines. Vertices standing close to each other (the pairs (45), (26) and (78)) belong to the same G. 
There is one internal B-contraction connecting the vertices 6 and 7 of the block. Two _B-contractions 
are leaving the block, which are connected to the vertices 1 and 2, running either to the right (vertices 
1 and 2 for the left diagram and vertex 2 for the right diagram) or to the left (vertex 1 for the right 
diagram). All the B-contractions which cross the block are summarized by the upper line. E is the 
sum of the Laplace variable E plus the sum of all x-variables of the B-contractions crossing the block. 
Therefore, for the left diagram, the resolvent corresponding to the left vertical cut is given by -g^^rtjiW}^ 
and the diagram corresponds to a contribution for {G^)\^'^ {E). In contrast, for the right diagram, we 
get for the resolvent, giving a contribution for {G^)\2^^{E ~Xi). Obviously, by shifting 

E ^ E + xi, the contribution from the right diagram coincides with the one of the left diagram. 



and interpreting dA as the 7"^-part which has to be integrated out. We can then write for 
the B-quantities 

(221) 
(222) 
(223) 

from which the RG equations can be read off once a formula for the _B-quantities has been 
derived. 

To be general, we start from a perturbation series of the form (|219p . where the LiouviUian 
and the vertices have already a dependence on the energy variable E. We denote the initial 
values by Ls{E), G{E) and R{E) (corresponding to Lj{E), G^{E) and R'^iE) if we consider 
the infinitesimal step A ^ A — dA of a continuous RG scheme). Replacing each contraction 
by the sum of 7"* and 7^ , we obtain diagrams where the vertices are either connected by A- 
or B-contractions. We can group each diagram into sequences of A-irreducible blocks, where a 
vertical cut between two vertices hits at least one A-contraction. A single A-irreducible block 
can have additional ^-contractions, which either connect vertices within this block, or connect 
this block with another block, or cross over this block, see Fig. [13] for illustration. Let us label 
the indices of the vertices of those B-contractions, which connect the block to other blocks, 
by 1, 2, . . .n (in the sequence from left to right) and call them external indices. Obviously, for 
71 > 1, we interpret the diagram as a contribution to the effective vertex {G^)\^"^" {E) (if 
the first vertex is G) or to {R^)\^"^" {E) (if the first vertex is R). The energy variable E is 
chosen such that the resolvent standing left to the block is of the form , according to 

the diagrammatic rule that the energy variable of the effective vertex must be identical to the 
one of the resolvent standing left to the vertex. Thereby we get the same result, independent 
on whether the B-contractions run to blocks to the left or to the right of the considered one, 
see Fig. 1131 Therefore, we determine the diagram by the convention that all i?-contractions 
go to the right. In this case, the energy variable E is identical to the one corresponding to 
all ^-contractions crossing the block. Subsequently, we have to sum over all possibilities to 
commute the external indices 1, 2, . . . , n within the block, see Fig. [TH except when the indices 
belong to the same vertex G or R. To each permutation we associate a corresponding fermionic 
sign for two reasons. First, the effective vertex {G^)\^"^"^ [E) should be antisymmetric for 
fermions under permutation of two indices. Secondly, if the external vertex is used as input for 
the perturbation series (|219|1 . it is assumed that the B-contractions leave the effective vertex 



Ll{E) 
G^{E) 
R^{E) 



L 



A-dA 



(E) 



LiiE)-'JmdA 



qA^cIA^^^ = G^{E) - 
R^"^'^{E) = R^{E) - 



dA 
dG^{E) 

dA 
dR^{E)_ 
dA 



dA 
dA 
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Fig. 14. The same diagram as shown in Fig. [13] but the vertices 1 and 2 are interchanged. Obviously, 
since the vertices 1 and 2 do not belong to the same vertex G, this is a different diagram to the effective 
vertex {G^)i2^'^ (E), which has to be counted separately. An additional minus sign has to be considered 
for fermions since the effective vertex assumes the sequence 12 for the reservoir lines leaving the vertex. 



in the sequence 1,2, . . .rt for the determination of the correct sign when this effective vertex 
is connected to other effective vertices. Therefore, if the original sequence differs from the one 
used in the effective vertex, as is the case e.g. in Fig.[T4j the corresponding fermionic sign has to 
be corrected by appropriate rules determining the effective vertex. Furthermore, we attribute to 
the effective vertex all minus signs which occur when the A- and ^-contractions are separated 
within the A-irreducible blocks from each other and from the external indices (again this is an 
information which is no longer available when the effective vertex is used in the perturbation 
series ()219|) ). In this way, we get for the diagrams depicted in Figs. fT3l and [Ml the following 
result 

(G^)i2(i?) ^ (±) (±) 731 75I lir G^iE) 7T3 G^{E + x^) TTig (224) 

■G45iE + Xi3)ni5G26iE + Xi5)ni256 Gvs{E + Xi256) ± (1^2) , 

where we have omitted for simplicity the obvious notation of the Keldysh indices. The resolvents 
and energy variables are defined via the short-hand notations 

ni...n = TbTTT— T ' Xi...n =Xi + ...+ Xn , (225) 

where we used the full effective Liouvillian in the denominator, see below. The two sign factors 
in (|224p arise when the two A-contractions are separated. As one can see, the rule to translate 
the effective vertex is quite simple. The indices of the energy variables of all vertices is always 
identical to the indices of the left preceeding resolvent. The resolvents get the indices of all the 
left vertices of contractions crossing the vertical cut at the position of the resolvent (the indices 
of all ^-contractions crossing the diagram and not belonging to the block are contained in the 
energy variable E). 

If no ^-contractions leave the A-irreducible block, we interpret the diagram as a contribution 
to the effective Liouvillian Lg{E) (if the first vertex is G) or to R^^q{E) (if the first vertex is 
R) . An example is shown in Fig. [151 which translates to 

LsiE) ^ i 7m l23 158 lii Gi2{E) i7i2 GzAz{E + X12) i75 Gq{E + X5) TTse G7s{E + xse) , 

(226) 

where the factor ^ is a symmetry factor since two vertices are connected by the same con- 
traction within the A-irreducible block. As indicated by the definition of the resolvent (|225p . 
such diagrams can be resummed on each line connecting the vertices within the A-irreducible 
blocks, leading to the full effective Liouville operator Lg{E) occuring in the denominator of 
the resolvent. In this way a self-consistent equation is obtained. 

Summing up all ^-irreducible diagrams, we obtain the following perturbation series for the 
effective Liouvillian and the effective vertices 

G^{E) ^ G{E) n^''). ^(^) (227) 

G{E + Xi) . . . i^^— — G{E + Xr) , 



E + Xi - Ll{E + Xi) ' ' E + Xr-Ll{E + Xr) 
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Fig. 15. An example of an A- irreducible diagram contributing to the effective Liouvillian Lg(E). 

where the first term on the r.h.s. represents the original quantities, and the second term con- 
tains all A-irreducible contributions, indicated by the sub-index A — irr at the product of the 
contractions. As derived above, we summarize the additional rules 

1. All free S-contractions are directed to the right. Free ^-contractions are not allowed. 

2. The external indices of the free ^-contractions are numerated as 1,2, ...,n from left to 
right. 

3. Sum over all permutations of the external indices and assign a minus sign for fermions 
according to the parity of the permutation. Omit permutations of indices belonging to the 
same vertex. 

4. Associate a fermionic sign which arises when separating all A- and i?-contractions from each 
other and from the external indices. 

Inserting the result (|227p for the ^-quantities into the effective perturbation series (|219p for 
the physical quantities of interest, we obtain three different terms. The first term on the r.h.s. 
of ()219|1 contains via the two terms on the r.h.s. of ()227|) the original quantities Ls{E), G{E) 
and R{E), together with all A-irreducible diagrams. The second term on the r.h.s. of (|219[) 
contains all ^-reducible but i?-irreducible diagrams. In this way one can see that the derivation 
of a single RG step is nothing else than an obvious classification into different topological 
sectors, combined with a convenient resummation of A-irrcduciblc contributions into effective 
quantities. 

We note that the effective vertices and are again in normal-ordered form, i.e. in 
the effective perturbation series (|219|) no i?-contractions are allowed between field operators 
belonging to the same _B-vertex. The reason is that we have incorporated all such terms already 
into the definition of the effective vertices and the effective Liouvillian by allowing for B- 
contractions on the r.h.s. of Eq. (|227|) . In this sense, our RG scheme is quite similiar to the 
RG scheme developed by Salmhofer [51] for quantum field-theoretical problems described by 
usual Feynman diagrams. However, besides the fact that we have to deal with operator vertices 
here and have to use the Laplace instead of Fourier transform, an important difference is also 
that we use the full effective Liouville operator Lg{E) as input for the dynamics of the local 
quantum system. In this way a self-consistent equation arises for the determination of Lg [E) 
via (|227p . Transferred to the Salmhofer RG scheme this would mean in a rough sense that 
self-energy insertions should be resummed on all propagators connecting the vertices, similiar 
to the RG scheme of Wetterich [50] (disregarding the fact that both schemes are quite different 
for many other reasons). 

Finally, we mention that one can also set up a non normal-ordered version of the RG 
procedure. This can be easily achieved by just forbidding the occurence of _B-contractions on 
the r.h.s. of (|227p . This means that no S-contractions are allowed within the ^-irreducible 
blocks which connect the block with itself, see e.g. the B-contraction between vertex 6 and 7 in 
Fig. 1131 These diagrams are then considered in the non normal ordered effective perturbation 
series (I219P by allowing for ^-contractions connecting field operators from the same i?- vertex. 
Concerning usual quantum field-theoretical problems, this would correspond to the RG scheme 
developed by Polchinski [55] together with taking the full propagators between the vertices 
including self-energy insertions. However, as one can see, the normal-ordered version includes 
more diagrams into the effective quantities and, therefore, the results are expected to be better. 

Continuous RG. We can now easily derive the rules for setting up the RG equations 
()166p - (|168p for the continuous RG flow. For a single infinitesimal step A-^ A — dA, we use the 
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Fig. 16. RG diagrams for the Liouvillian and the vertex up to 2-loop order for a model with charge 
fluctuations. The slash indicates the derivative of the contraction 



decomposition (|220p with 



= ^ , ^ A-dA (228) 



Comparing (j^^ - ip^ with (with Ls ee L^, G = and R = R"^), we obtain immedi- 

ately 



R^{E) J ^ [R''{E) J 

1 ^A f IP I V \ ^ r^A I 



G'\E + Xi) . . . -r— — &\E + X,) 



E^Xx-L^{E-\-Xx) ^ "' E + Xr- LjiE + Xr) 



where, for dA — > 0, we can omit all terms ~ (dA) with fc > 1 on the r.h.s., so that only one 
contraction ^ can occur, and we can replace 

Lf = L^"''^ ^ , = G^-"^^ ^ G^ , i?^ = i?^-"*^ i?^ , 7-^-"^^ (230) 



in all other parts of the diagram. Since the total expression has to be irreducible with respect 
to ^ and only one such contraction is allowed to occur, we get the simple rule that the 
contraction must contract the first with the last vertex in the diagram, see Figs. [T51 [TTl and 
[T5]for examples. Diagrammatically this contraction is indicated by an additional slash. 

To illustrate the rules, we have shown in Fig. [16] the RG diagrams up to 2-loop order for an 
arbitrary model with charge fluctuations where the vertex depends only on one index. Omitting 
the index A and the trivial Keldysh indices, and using the notation 

711' = 71 (231) 

to exhibit explicitly the (5ij/-part of each contraction, the RG equations are given by (using 
again the elegant notation (|225p ) 

-^ = ^Gi(i^)^iGi(i; + .r) (232) 

+ P^j2GiiE)niG2iE + xi)nuG2{E + xi2)niGiiE + xi) , (233) 
dA 

- ^^¥r^ = ± ^ ^^i^) ^2 Gi(i? + X2) n,2 G^iE + X12) . (234) 
dA dA 

Note the sign factor in front of (|234p due to separating out the contraction between index 2 
and 2 from the external vertex 1. 

Figs. [T71 and [T51 show the RG diagrams up to 2-loop order for an arbitrary model with spin 
and orbital fluctuations ((44|l . the Kondo problem (j45|) being a special example. In this case we 
get the RG equations 

- ^^1^ = ^ ^2 G^2{E) n,2 G-,-AE + x,2) (235) 

+ ^ 72 73 Gi2{E) 7Ti2 G23(-B + 2:12) ilia G-^iiE + 0:13) (236) 
dA 
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Fig. 17. RG diagrams for the Liouvillian up to 2-loop order for a model with spin and/or orbital 
fluctuations. The slash indicates the derivative of the contraction 4t- 
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Fig. 18. RG diagrams for the vertex up to 2-loop order for a model with spin and/or orbital fluctuations. 



The slash indicates the derivative of the contraction 



for the Liouvillian, and 
dGn'{E) 



dA 



dA 



d'-f2 



+ -jj 73 G23iE) 1123 Gir{E + X23) -^11-23 G32{E + xii'23) 

± S^^l3Gi2{E)ni2Gv3{E + Xi2)nn'23G32{E + Xn'23) ± (1^1'] 
^l3G23iE)n23G3,iE + X23)ni2G,,2iE + Xi2) ± (l^l') 

dA 



(237) 
(238) 
(239) 
(240) 



for the vertex. The corresponding equations for R follow from replacing the first vertex G — > i?. 
Note the sign factors in front of (|239p due to separating out the contraction between index 2 
and 2, and the missing permutation of 1 and 1' in (I238|) since both indices belong to the same 
vertex. We again note that the energy arguments of the vertices are always such that the indices 
of the x-variable coincide with the indices of the left resolvent. Using these rules, it should be 
obvious how to write down RG equations in arbitrary order for generic models nearly "without 
thinking". However, the problem of course is to solve these RG equations which will be the 
subject for the rest of the paper. 



4.3 Summary of RG equations and general properties 

Summary. We first summarize the RG procedure described in the previous sections. As out- 
lined in Sec. l4.1[ we use a 2-stage RG procedure. In the first step we integrate out the symmetric 
part of the distribution function by a discrete RG step using the decomposition (|196p 



ill' = 7? , 7? = 
where we use from now on the conventions 



1 (- 



2 1 1 



UJ ^ rjLU , Ha = TlPa , X = r]{uj + Pa) ^ LU + flo 



and 



1^1. ..n — LUl + ■ ■ ■ + (jJn , Pl...n 

For fermions and T — we have 



Ail 



^l...n 



Xi 



= - ^sign(tj) 



+ Xn 



(241) 
(242) 

(243) 
(244) 

(245) 
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The effective Liouvillian Lg and the effective vertices and i?° are given by (|227p 

(246) 





f Ls 


cIe) \ - 




R^iE) J 






1 




G ... - — 1— — —G 



E + Xi- U'g{E + Xi) E + Xr- L^iE + Xr) 

where s — irr means irreducibiUty with respect to the symmetric contractions 7". The vertices 
on the r.h.s. are the initial vertices having no dependence on the energy variable E. The energies 
Xi = Xi „ contain all indices of the left vertex of all contractions crossing the i-th resolvent. 
We also use the convention 

El...n ^ E + fli,„n (247) 

such that the energies Xi can be written as 

E + X, ^ E + xi...„ = Ei,„„ + , (248) 

separating more clearly the dependence on the integration variables uji and the physical param- 
eters Ei ^n on which the final solution after integration will depend. 

For the next continuous RG procedure we need only the vertices averaged over the Keldysh 
indices. Therefore, we sum over all Keldysh indices on the r.h.s. of (|246p . Since the Keldysh 
indices of the initial vertices are all the same, see ([77|) . we obtain only the two possible combi- 
nations 

Gi...„ = J2 G'l.n' , Gi...„ = PGi'.n (249) 
p p 

on the r.h.s. of (|246p . Due to the form (|24ip of the symmetric contraction, G occurs if an odd 
number of symmetric and left-going contractions are attached to this vertex. 

In the second step we replace the remaining antisymmetric contraction 7° by the /l-dependent 
contraction (|215p with a cutoff defined on the imaginary frequency axis. By convention we write 
it in the form 

111' = ^11' li (250) 

with 



7i 



A _ ) V 



-1 



{vCo)tY^ _ \ ^ 0T„(^- Kl) (251) 



LJ — iLo: 

A 

1 r '^'-^''^^ 2^ ' 



( V \ ^^(^^) !_ f (252) 



The corresponding RG equations are given by (|229p 



dA 



G^iE)) ^ {gHE)\ (253) 



■G'^iE + Xi) ... \j— — G'^iE + Xr 



E + Xi - L^{E + Xi) ' ' E + Xr - L^{E + Xr) 

where only the vertices (jl98p averaged over the Keldysh indices occur since the antisymmet- 
ric contraction (|25ip does not depend on the Keldysh indices. The initial condition for these 
differential equations at yl = c» is given by 

Lt°^\E) = L%{E) , (254) 
{G^=ni...n{E) = GI,JE) , (255) 
{R''=ni...niE) - Rl„„{E) , (256) 



Will be inserted by the editor 



47 



where Lg, G°, and i?° are given by (|246p from the first discrete step. 

Solving the RG equations down to yl = gives the final result for the physical quantities 
via (fT63l) 

Ll"iE) = LfiE) , Sj,{E) ^ RfoiE) , (257) 

or, alternatively, one can stop the RG at any arbitrary value of A and use the effective pertur- 
bation series (|217p to calculate the physical quantities 

1 _ 1 _ 



■G'UE + Xi) ... J- -G''(E + Xr) 



E + Xi - Lj{E + Xi) ' ' E + Xr-L^iE + X^ 

Of course this makes only sense if this perturbation theory is well-defined at the scale A, i.e. if 
all large terms have already been eliminated by the previous RG flow (which happens when A 
reaches some physical low energy scale and the couplings are still small, see the discussion in 
the following sections). 

Symmetry relations. All the symmetry properties from (anti-)symmetry (see (|87|) ). con- 
servation of probability (see (|114p and (|85p). and those following from the hermiticity of the 
original Hamiltonian (see (|115p . ([93]), (|140p and (I132p ). are also fulfilled for the effective quan- 
tities. In summary we have the following relations 

G^::i^:j:l-'"iE) = ±G^:f (259) 
R^::a:'r'"iE) = ±i^::^-:'^-'"{e) , (260) 

TTsLsiE)^0 , TrsG?-:/"(i?)=0 , (261) 

Ls{EY = -Lsi-E*) , (262) 
Gll-'-iEy ^ -{a-TG^-ti-E*) = -a---G\]--ti-~E*) , (263) 
RV.:T'{Er = -{a-rit-,Ti-E*) = -a—'-Rf-t^-E*) , (264) 

where we can either use L%,G'',R°' or Lj,G'^,R'^ for Ls,G,R (for G and R, the Keldysh 
indices are of course omitted), a convention we will also use for the following discussions. These 
identities follow directly from the initial symmetries and the RG equations (|246p and (|253p , see 
Appendix D for the proof. 

We see that all symmetry properties are preserved under the RG flow. Moreover, they are 
even invariant within all truncation schemes, since they hold for each term on the r.h.s. of the 
RG equations separately, provided the complete sum over all indices is taken. 

Analytic properties. For a subsequent discussion of the frequency integrations, we first 
study the analytic properties of the Liouvillian and the vertices in the variables E and xi . . .Xn- 
For this, we need the essential ingredient that 

— - is analytic in E in the upper half plane , (265) 

E — Ls[E) 

which is equivalent to the property that the effective reduced density matrix, defined by 

Ps{E) = Psito) , (266) 

has the analytic properties depicted in Fig. 1101 i.e. contains no exponentially increasing solutions 
in time space. To prove this statement generically is difficult since it is not clear whether 
certain approximations for Ls{E) can not lead to poles in the upper half plane. The question 
is ultimately related to whether the relaxation and dephasing rates generated by the RG flow 
are positive. So far, we have not seen any mathematical proof to show this property from the 
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structure of the RG equations. Of course we know that, if all terms are taken into account on 
the r.h.s. of the RG equations, wc get the exact solution and at least at the end of the flow, 
the analytic property (|265p must be fulfilled due to physical reasons. Whether it is also fulfilled 
throughout the RG flow and within certain approximations of the RG equations is not clear 
and an important question for future research. We will assume in the following that it holds. 

From (|265|) and the analytic structure of the RG equations (|246l) and (|253p . we obtain 
immediately that 

Ls{E),G{E), R{E) are analytic in E in the upper half plane . (267) 

Furthermore, we can also study the analytic properties of the vertices in the frequency 
variables xi . . .Xn with Xi — rjiiuji + //„.) (a better choice than the frequencies uji because 
always this combination occurs in the perturbative series). Writing the vertices as function of 
the a;- variables 

) , (268) 

and assuming that the vertices are analytic in the a;i-variables initially, we can prove the fol- 
lowing analytic property which is preserved under the RG 

GiE + Xi + ...+ X„i; -Xi ...-Xm Xm+l - --Xn) 

is analytic if all variables E, xi . . . Xn lie in the upper half plane , (269) 

where 1 < m < n can take any value. The same holds if any permutation of the indices of the 
a;- variables is chosen. It means that the vertices stay analytic in the upper half plane if the sign 
of some Xi-variable is changed, provided one replaces E E + Xi. The property follows from 
the structure of the RG equations since only the frequency combinations of (|269p occur on the 
r.h.s. Consider e.g. the first term of the RG equation (|237p . which we write in the representation 
()268|) and integrate out the 5- functions of the contractions by using ()250p . This gives 

^ G,, (i?; a:ix'i) ^ • (270) 

• Gniux,ri2V2 [E] X1X2) ■ - / „ ; r G^ri2P2,ri\ y', {E -\- Xi -\- X2 \ —X2X1) , 

-tj + Xl + X2 — l^Syt^ + Xl + X2) 



or for E E + xi and xi — > —xi 
-^G^,^,^^'^^'^{E + xi;-xix[) - 



-^G^^^,^^'^^'^{E + xi;-xix[) ^ • (271) 



• G ,ri^u2 {E + Xi] —X1X2) — 77^ T G^ri2V2.ri'.v' {E + X2] —X2X1) 

-Cj + X2 — Ls(-tj + X2) 

In all cases we see that the analyticity in E, xi and x[ is preserved. This holds always because 
the sign of a a;-variable can only change due to the S12 S{xi + X2) parts of the contractions. 
Integrating out the right index of each contraction 712 leads to the argument ~xi for X2- 
However, the vertex G containing the index 2 stands to the right of the resolvent which the 
contraction 712 crosses, see Fig. [12] for illustration. Therefore, the energy arguments of this 
vertex are G{E + xi + ...;... — xi .. .). This shows that the structure of the energy arguments 
is generically of the form shown in (|269|) . 

These considerations prove that the resolvents and vertices are even analytic with respect 
to all the internal x-variables, e.g. the a;2-variable in the example above. This will turn out 
to be quite useful to discuss the frequency integrations over the internal variables. A generic 
discussion of the frequency integrations is of course not possible since the frequency dependence 
of the density of states Puiy^) and the initial vertices can differ for each model, see e.g. the spin 
boson model discussed in Sec. 12.21 We discuss here the fermionic case with a flat density of 
states, i.e. the frequency integrations are performed from minus to plus infinity and we cut off 
the frequencies of all contractions by the Lorentzian function (|47|) 



^(^) = = - ^-7^) ■ (272) 
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2 3 



Fig. 19. Illustration for the structure of the energy arguments of the vertices. According to the dia- 
grammatic rules, the resolvent and the right vertex have the structure ^^^^^^ G{E + xi; X2X3). 
Integrating out the S12 ~ S{xi + X2) part of the contraction 712, the energy arguments of the vertex 
get the form G{E + xi; —xix^). 



1 
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1 1' 



1' 1 



1 1' 



1' 1 



Fig. 20. Examples for diagrams of the effective Liouvillian when the symmetric part of the contraction 
is integrated out in one step, s (a) denote the symmetric (antisymmetric) contraction 7^ (7")- Whereas 
the first diagram is finite for D — » 00, the other two diagrams contain also terms ~ D. 



Furthermore, the initial vertices are assumed to be analytic functions of the frequencies. 

Frequency integrations for the discrete step. We start with the discussion of the 
frequency integrations for the first discrete RG step where the symmetric part of the contraction 
is integrated out. Using the analytic properties (|265p of the resolvents, the integrations can be 
performed analytically. When a symmetric contraction crosses over n resolvents in the series 
(|246p . we get an integral of the type {zi are variables with positive imaginary part) 

J '^^(^^ n - + - + z.) = n .i? + - LU^D + z.) ' ^'''^ 

i.e. uj is just replaced by iD in all resolvents. This is a very large imaginary part for each 
denominator and, therefore, this leads to a well-defined perturbation series which can be ex- 
panded systematically in the vertices and in 1/D. Therefore, it is not necessary to take the full 
Liouvillian Lg in the denominator and solve the equations self-consistently, but one can replace 
it by the initial value Ls and consider the additional diagrams from s-irreducible insertions 
perturbatively. 

Examples of diagrams for the effective Liouvillian are shown in Fig. 1201 The first diagram 
correponds to charge fluctuations and reads for fermions, using the notations ([243]) . ([2441) and 
(IMTH 

iGr|d^p(^)_^^|_^^Gi 
^1g,Gj+0{^) . (274) 

Neglecting the contributions ~ we find a result independent of D, indicating that such 
models of charge fluctuations (where the vertex Gi has only one index) have a well-defined 
limit for Z? ^ cxd if all charge states of the quantum system are taken into account. If the 
symmetric contraction crosses over more than one resolvent, the result is ~ and can be 
neglected. Thus, for such models, the determination of the initial condition is very simple. 

In contrast, the second and third diagrams of Fig. [20] are associated with a model where 
only spin and/or orbital fluctuations are considered, the Kondo model being a special example. 
For these models, certain charge excitations have already been integrated out, so that a finite 
bandwidth D has to be introduced and the limit Z) — > 00 is not allowed. In the second diagram 
the following integral occurs 

/did p(uj) f du)' piu)') = / duj niuj) = — 

' J ' uj + uj' + z J ' uj + iD + z z + 2iD 



p'llGl 
2 iD 
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El -I- uji 
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1 2 



2 1' 



Fig. 21. Examples for diagrams of the effective vertices when the symmetric part of the contraction is 
integrated out in one step, s (a) denote the symmetric (antisymmetric) contraction 7" (7°). Whereas 
the left diagram is negligible for D — > 00, the right diagram is finite. 



D 



1 tz 

^ 2D 



(275) 



which gives for the diagram the value 
1 



1 



Ls 



1 



; Gil' 

TT^ - - 

^t — DGiv Gvi 



UJl 

duj' piu)') — — 

OJ + 



1 



+ Ell' — Ls 
Ls) Gi'i 



Gi'i 



o,i, 



(276) 



As we can see also terms ~ D are generated which are unphysical since one can prove that the 
model should contain only logarithmic divergencies ~ (In(D))"' in higher-order perturbation 
theory. Such terms occur because the high temperature limit is not well-defined (only temper- 
atures T << D are allowed since otherwise charge excitations become important which are not 
considered in the model). As we will see later for the Kondo model in Sec. O the terms ^ D 
will be cancelled by corresponding terms generated in the second continuous RG flow. 

Similiarly, also the antisymmetric part of the contraction can occur in the diagrams, see e.g. 
the third diagram of Fig. [20l It contains more complicated integrals of the form (evaluated here 
at zero temperature but the same result comes also out at finite T <C D) 



duj p{uj) 



duj' p{uj') 



1 



D2 + (z + iD)^ ^''^^ 



LU + LU' + Z 

iz ^ 



duj p{lu) 



D 



ttD 



^ 2D 



ln(l 



U! - 

iz 



iD + z 



= ^lD-^^z + 0{-) . 
2 2 ^d' 

This leads to the following value for the third diagram of Fig. [20l 

1 



(277) 



LUE) 



p'lirvGZ, 
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Ell' 
duj p{uj) 



- UJl ^ 



— D Gil' Gi'i 



i-^Giv {Eiv 



Ls 

duj' p{uj') — 

UJ 

- Ls) Gi'i 



Ell' — Li 



Gi'i 



0,1, 



(278) 



A similiar structure as compared to the other diagram p76p occurs. This illustrates that the 
symmetric part can be integrated out using perturbation theory in the coupling vertices. No 
logarithmic divergencies ~ ln(D/|£'i...„ — Ls\) can occur in any order since the upper limit 
of all integrals and the cutoff provided in the denominators of all resolvents is given by the 
bandwidth D. 

Examples for diagrams for the determination of the effective vertices are shown in Fig. 1211 
The first diagram corresponds to charge fluctuations. It is of order 0{-^) and can be neglected 
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2 J UJ2 + £^2 - -l^S 0J2 + rjl2 + OJl — Ls D 

The second diagram corresponds to spin/orbital fluctuations and is finite in the limit D — > oo 
GME) - P'1'2GII \ -Gf;^' - (1 - 1') 

£/12 + Uli + UJ2 — Ls 

1 - f _ _ 1 

= i G^2 ^ ^ G,,, - (1 ^ 1') 

2 iD + Ei2 + LUi - Ls 

= (Gi2G2r - GP2G21) + 0{^) . (280) 

In conclusion, we see that for initial vertices with one index (charge fluctuations), the effec- 
tive Liouvillian is finite for D —^ 00 and the effective vertices are identical to the initial ones. 
For initial vertices with two indices (spin and/or orbital fluctuations), the effective Liouvillian 
contains also terms ~ D and the effective vertices get a finite correction to the initial ones. 
Initial vertices with three indices correspond to quite artificial models, in that case also terms 
~ can occur for the effective Liouvillian and terms ~ D for the effective vertices. 

Furthermore, we note that the effective Liouvillian can contain terms which are linear in 
the energy E, see p76p and (I278|) . This does not lead to a problem for arbitrarily large frequen- 
cies because the effective Liouvillian occurs always in the denominator of all resolvents of the 
effective perturbation series or the r.h.s. of the RG equations. In contrast, the effective vertices 
are independent of frequency, see ()280|) . which is important for a well-defined large frequency 
behaviour of all expressions. 

Frequency integrations for the continuous RG equations. We now turn to the fre- 
quency integrations of the RG equations (|253p of the continuous RG flow. From the form of 
the contraction (|251|) . we see that by closing the integrations over <D in the upper half of the 
complex plane, only the pole from the Matsubara frequency and the pole of the Lorentzian 
cutoff function (|272p at iD contributes (we again replace p,^{uj) by (|272p in all contractions). 
All other resolvents and vertices are analytic in the upper half plane, according to p65p and 
(|269p . Most importantly, the contraction which is differentiated with respect to A, crosses 

over all resolvents (to make the expression ^-irreducible), i.e. the frequency variable of that 
contraction occurs in a// resolvents! Moreover, due to the derivative, we get for this contraction 
from (P5T|) 



drfl 
dA 



ll p(^)tY^^-^StAA-\<\) (281) 



{-^l}^(-)2l 



1 1 



UJ — iAt^ oj + iAj 



(282) 



where 5t{uj) is a delta function smeared by temperature, defined as the derivative of 9t{(jj), 
given by (|203p . i.e. 

X ( \ ^ a ( \ Jo for |w| > ttT 

As a consequence, only the two Matsubara frequencies ±^Tq lying closest to ±yl contribute to 
the sum in Eq. (|28ip . and we obtain exactly Eq. (|282p . Therefore, we get the important result 
that, by closing the frequency integration of ^ in the upper half plane, the contribution from 
the Matsubara poles gives rise to the replacement 

UJ iAx^ in all resolvents of the r.h.s. of the RG equation ([253]). (284) 
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Therefore, a large positive imaginary part occurs in all resolvents so that the perturbative series 
on the r.h.s. of the RG equations is well-defined provided that the coupling constants remain 
small. In contrast, the Matsubara poles from all the other contractions occur only in those 
resolvents where the contractions cross over. 

The pole at iD from the density of states is only important if one wants to study the precise 
behaviour of the RG flow when A crosses the bandwidth D. However, nothing really interesting 
is happening in that regime, except for the fact that D sets the scale where the divergencies of 
the original perturbation theory are cut off at high energies, i.e. at yl ~ D the RG flow starts 
to renormalize the effective parameters significantly. To see this analytically, we perform the 
integration over the frequency variable cD of the ^ contraction analytically. Using (|282p with 
P72|l we obtain by closing in the upper half of the complex plane 

/ du. piu,) 1- ( _ \. + _ J.. ) {. . = 

J Ztt to — iAt^ uj + iAt^, 

- {■■•L=J . (285) 



Al^ - £12 



where {. . .}a stands for the rest of the diagram. We see that for A':s> D, there is small prefactor 
^ (-j)^ ^ 1 in front of the r.h.s. of the RG equations leading to a negligible renormalization. 
In contrast, for A the prefactor is i and the renormalization becomes important. In this 

regime, only the term {. . .}i,=iyiT„ to be considered on the r.h.s. of (I285|) and the r.h.s. of 
the RG equation becomes independent of D. 

In summary, the RG flow will roughly start a,t A ^ D but the precise value where it starts is 
not so important since the cutoff D occurs logarithmically in the divergencies of perturbation 
theory. Therefore, it is a good approximation to evaluate the r.h.s. of the RG equations for 
A <^ D, and use D just as the starting value of the flow parameter yl, where the initial values 
of the Liouvillian and the vertices are taken from the first discrete RG step. Moreover, as we 
will sec in more detail in connection with the application to the Kondo problem in Sec. [5l the 
precise ratio of the initial value of A and the bandwidth D should be chosen in such a way that 
the terms linear in D of the initial Liouvillian cancel out (as it should be since such terms do 
not exist in the original perturbation theory). 

As a consequence we consider from now on only the regime A D and omit the cutoff 
function p{u}) in all contractions. We then close the integrations over all variables uji in the 
upper half plane and replace them by the pole of the Matsubara frequencies in all resolvents 
and vertices of the diagram. As a result we find the same form of the RG equations but with 
the replacements 

(286) 



dA 



for all contractions, together with 



UJ IUJ„ 



I du! ^ 27rr„ 6It„(v1-0 (287) 

0<a;° 



for the frequency variables of all contractions 7j^, and 

u iAr^ (288) 

7 A 

for the frequency variable of the contraction -jj- . Note that only positive Matsubara frequencies 
arc allowed in the sum. 

RG on Matsubara axis. When the integrations over all real frequencies are performed, 
the energy variables E + X; of the RG equations (|253p get the form 

E + Xi ^ E + /2i..,„ + iAr^^ + iuj2...n , (289) 
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where index 1 is assumed to be the variable of -j^, and 

uJi...k =<+...+ < (290) 

is the sum over the Matsubara frequencies in analogy to the notation (|244p . We see that the 
value for the energies E appearing in the Liouvillian Ls{E) and the vertices Gi. ..„(£') are 
complex numbers with the imaginary part being some positive Matsubara frequency. Thus, we 
can reformulate the whole RG equations for frequencies on the Matsubara axis by replacing 
all variables uj by Matsubara frequencies and using for E a real part plus some Matsubara 
frequency. With this, only the following quantities occur in the RG equations 

Ls{E;uj) = Ls{E + ioj) , 

Gi,,,niE,UJ,UJi . . .UJ„) = G.,,^^^,,,r,,^^„{E,UJ,UJi . . .UJn) ^ Gi„,niE + iuj)\^^^i^. , (291) 

where uj,uji . . . cOn are positive Matsubara frequencies and E is real. By convention, if the 
frequencies are written explicitly in the arguments, the index 1 = rji/ does no longer contain the 
frequency. The RG equations in Matsubara space for Ls{E,uj) and G,,,{E,uj, . . .) are obtained 
from the original RG equations by the replacements 



dA i -1 , 

7Ti...„ ^ i7(i;i...„,c^ + wi...„) , (293) 

Gl..,n{E + Xi>,,,n') Gi,,,n{Ei',,,n>,UJ +UJi',,,n',UJi . . .UJn) , (294) 

ui ^ At^^ for $i , (295) 
^ dA 

I dio ^ 2T:Ta ^ OrM-i^n) for 71 , (296) 

0<a;° 

where all frequencies correspond to Matsubara frequencies LOi = uj"^ , and we have used the 
convention Ei,,,n = E + /2i...„ together with the definition 

^(^'^) = IT—- \ fj, , ■ (297) 

E + luj — Ls[E, Lu) 

Furthermore, the explicit frequency arguments of the resolvent and the vertices in Matsubara 
space can be avoided by the compact notation 

Gfl,,{Er,,,r,') = Gi...n{Er,,,„',uj + uji,,,,r,',uji...uj,,) , (298) 
7rf^.„ = i7(i;i...„,c^ + cc;i...„) , (299) 

which is especially useful for writting down higher-order terms of the RG equations. In this 
notation the Matsubara RG equations have nearly the same form as the original RG equations 
besides the replacements (|292p for the contractions. However, one should note that the external 
frequency variable u; is hidden in this notation and one has to specify explicitly which variable 
is set on the parameter At^ . 

As an example consider the RG equations (|232p - (l234p for a model with charge fluctuations. 
Using the above rules in the representation (|29ip . we get immediately 

- '^^^^fj^^^ = ~iGi{E,uj,AT^^)n{Euoj + AT^JGi{Ei,oj + AT^^,-AT^^) (300) 
~-Gi{E,u},AT„^)n{Ei,bj + AT^^)G2{Ei,bj'\- AT^^,uj2)n{Ei2,uj + At^^ + ljs) • 

■G2{Ei2,yJ + At^^ + uj2,-uj2) n{Ei,Ljj + At^^)Gi{Ei,ljj + At^^,-At^^) (301) 

_dGiCB^_av^ = iG2{E,Lo,AT^^)n{E2,oj + AT^^)Gi{E2,oj + AT^^,m) ■ 

■n{Ei2,oJ + oJi+ At^)G-2{E^2,^ + ^i+ At^.-At^ . (302) 
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Using the compact notation p98p and (|299p . the RG equations P35p - p40p for a model with 
arbitrary spin and/or orbital fluctuations read for fermions 

+ ^Gr2iE)nf^^GliiE^,)n^iGliiE,,)U=^,^^ , (304) 
_ dGME,^co,oj„u[) ^ ^ {G^i^E)4^Gli,iE,,) - (1 ^ l')} (305) 

- GIUE) Gf(, {E23) <,23 Gi(i?ll'23) (306) 

+ {G',iiE)4iG^f,iE,,)n^i,,,Gli{En'23) - (1 - l')} (307) 

- {G^3(^^)'r2'3'Gf((i?23)^f|Gf4(i?i2) - (1-1')} U=^.„, . (308) 
For the final physical quantities, we need according to (|257p 

Ll"{E) = Lt"{E,iu^O) , Ur{E) ^ Rt"o{E,uj^O) , (309) 

i.e. we can fix E in the RG equations and see that only E shifted by the discrete set of values 
P'i...n occurs for the first frequency variable. This means that the RG equations are local in 
Laplace space up to a shift by multiples of the chemical potentials of the reservoirs. This 
simplifies the calculation of the time evolution considerably. E.g. if only two reservoirs with 
/ix, = — /ifl = Y present, we get with n = 0, ±1,±2, . . . for i.e. only the fixed 

Laplace variable shifted by multiples of half the bias voltage can occur for the first frequency 
variable of the Liouvillian and the vertices. 

In summary, we find that the RG can be formulated on the Matsubara axis but in contrast 
to equilibrium Matsubara formalism, there is a whole set of Matsubara axis shifted by the real 
quantities 

E + nini + . . . + Tiz HZ , = 0,±1,±2, . . . (310) 

where Z is the number of reservoirs, fj,k is the chemical potential of reservoir fc, and E is the 
real part of the Laplace variable. This is the price one has to pay compared to equilibrium 
to calculate the time evolution and the effect from different chemical potentials. If different 
temperatures of the reservoirs occur (leading to nonequilibrium heat currents), the Matsubara 
frequency axis has to be defined for each reservoir individually. 

Finally, we note that within our formalism no analytic continuation to real frequencies is 
necessary to calculate the time evolution of the average over an arbitrary observable (the same 
holds for correlation functions not shown in this paper) . This is an advantage compared to the 
usual linear response formalism |65j or to recently developed nonequilibrium formalism with 
complex chemical potentials |16) . where an analytic continuation is necessary. In this way we 
can avoid the numerical problems associated with analytic continuations. 

Cutoff parameters for the RG flow. Using the form (|289p together with the repre- 
sentation (|29ip . and replacing the effective Liouvillian by its eigenvalue via (|176p . we find the 
following structure of each resolvent of the RG equation 

-A J , (311) 

E + fii..,n - + iAr^^ + + iuj2...n + 

where uj is the Matsubara frequency associated with the energy variable E. The energy argu- 
ments of the real and imaginary parts /i^, of the eigenvalues of the effective Liouvillian at 
scale A are not indicated (which is the same energy as occuring in the rest of the resolvent). 
As already mentioned in Sec. 14.11 this form of the resolvents is very useful for stable numerical 
calculations since all imaginary parts in the denominator are positive. Furthermore, it shows 
generically what the various cutoff parameters of the RG ffow are. If the flow parameter A falls 
roughly below the value {E is assumed to be real here) 

yl < max{T„, |£; + /2i...„-/i^|,w,u;2, ... , (312) 
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the resolvent (|31ip becomes approximately a constant, cutting off logarithmic divergencies 
for low energies. In this way the temperatures Tq of the reservoirs, the Laplace variable 
the chemical potentials jia of the reservoirs, the oscillation frequencies hk^ the Matsubara 
frequencies w", and the relaxation and dephasing rates Fk serve as cutoff parameters. The cutoff 
occurs in a natural way and in an intuitively expected form in the RG equations. However, in 
contrast to temperature and the imaginary parts of the resolvent, one should note that only 
the combination E + /2i...„ — is the cutoff parameter in the real part of the denominator. 
Therefore, there are interesting points when this combination is zero 

E + Aii...„ - /ifc = , (313) 

where enhanced renormalizations are expected. These points correspond physically to reso- 
nances. E.g. for the Kondo model they lead to an additional logarithmic enhancement of the 
differential conductance when the external magnetic field is identical to the bias voltage, see e.g. 
Ref. J3J for a detailed study of this problem using RTRG-FS or Refs. |45I46I47| where other 
methods have been used). A nice feature of the RTRG-FS procedure presented here is that, 
without specifying the model at all, we can generically predict that logarithmic enhancements 
are expected when the condition (|313p is fulfilled. 

In contrast, temperature is always a cutoff scale because when A falls below T^, the last 
Matsubara frequency has been integrated out and the RG flow stops. Compared to equilibrium, 
this means for nonequilibrium problems that the voltage (inducing nonequilibrium stationary 
states) and the Laplace variable E (determining the time evolution of an initially out of equi- 
librium state) induce different cutoff behaviours of the RG equations than temperature. All 
resolvents on the r.h.s. of the RG equations have a different structure concerning the combi- 
nation of energy variables E -\- /Ii...„ — h^. In equilibrium, only the eigenvalues hk can occur 
but they can not be cancelled by other energy variables. This opens up many new interesting 
physical phenomena for nonequilibrium systems to be studied in the future. 

As temperature, also the Matsubara frequencies and the relaxation and dephasing rates will 
cut off the RG flow and can not be cancelled in the denominator since they are all positive. 
Most importantly, as has already been discussed in all detail in Sec. 14.11 the zero eigenvalue 
of the effective Liouvillian can not occur in the resolvents. This means that we have found a 
generic proof that relaxation and dephasing rates will always cut off the RG flow, irrespective 
of the truncation scheme used and of the specific model under consideration. 



4.4 Weak coupling limit 

Definition of weak coupling. We start with the definition of weak coupling by defining the 
dimensionless "coupling constant" J a at scale A 

JA = A'^-^Gt.n , (314) 

where, from now on, we use the definition (|20p for the field operators such that a^{uj) has 
dimension The connection between the RG equations for Ja and is 

where I is a dimensionless flow parameter running from / = Oatyl = £'to^ = ooatyl = 0. 
Here, the first term describes the trivial renormalization of J a by the rescaling factor in (|314p . 
whereas the second one stems from the renormalization of G^. If we rescale all energies in the 
RG equation for G^ by A, the second term on the r.h.s. of ()315|) is a power series in Ja starting 
at J A for n > 1 and at J a for n = 1, see the RG diagrams of Fig. [16] for n = 1 and Fig. [18] for 
n = 2. For a term ~ J^, there are k — 1 resolvents in between, so we estimate the term in this 
order as 

'^^ 1 0(1) for.l»zi ' ^^^^> 
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where A stands symbolically for all the physical scales occuring in the denominator of the 
resolvents according to (|31ip . i.e. 

A ^ T, E + fii,„n~h, r , (317) 

except for the frequencies which are rescaled by A and will give something of order 0(1) in 
the denominator of the resolvent (note that all frequency integrations are cut off by A, they 
all give a positive contribution to the imaginary part of the denominator, and iA occurs in all 
resolvents). 

Weak coupling means that the expansion (|316p of the expression A^ in a power series 
in J A is under control for all values of the flow parameter, which is the case if the condition 

Ja < 1 or yl < Z\ (318) 

is always fulfilled. Only in this case truncation schemes on the r.h.s. of the RG equations 
are justified, leading to a well-defined renormalized perturbation series in J. We discuss this 
condition now for the possible values of n in (|315p . where n is the number of field operators of 
the considered vertex G^ . „. 

For n > 2, the weak-coupling condition is certainly not a problem since already the first 
term of (|315p leads to a reduction of J. Therefore these terms are called irrelevant in the 
renormalization group sense. 

For n ~ 2 (e.g. spin and/or orbital fluctuations) we have J a — and the first term on 
the r.h.s. of (|315|) is zero. By convention, this is called the marginal case. The weak-coupling 
condition piSp will be fulfilled if the renormalized coupling Ja stays small until A reaches the 
cutoff scale A. Under this condition, Ja will also stay small for A < A, since an additional factor 
~ (^)'"'~^ appears in this regime, see (I316p . Since the smallest cutoff scale is the minimum of 
all relaxation and dephasing rates of the problem, we obtain the condition (|170p 

min {Fi} Tk weak coupling , (319) 

i 

where Tk is defined as the energy scale where the coupling Ja=Tk ^ 0{\) if all physical cutoff 
scales A are disregarded. 

For n — 1 (e.g. charge fluctuations) the first term on the r.h.s. of (|315p leads to an increase 
of J, therefore this situation is called the relevant case. Due to the relations ([39| and (|40p. we 
expect ^ V and we obtain 




(320) 



When A reaches the physical relaxation or dephasing rate -T, we expect ^ F, and thus 

Ja^f - ^/^/T ^ 0(1) . (321) 

Therefore, already before A reaches F, it may happen that J becomes of order one and a 
truncation of the series is not strictly justified. For A < F, we expect F^ to saturate to F so 
that Ja ^ \JF j A with 7^ being independent of A. Therefore, in this regime, we estimate (|316p 
to be 

J\ {A/Ff-^ ^ (/tm)' [A/Ff-^ ^ (^/A/rf^^ , (322) 

so that, for k > 2, we expect the terms to become small again. Since, for n ~ 1, at least three 
vertices are necessary for a renormalization of the vertex, see Fig. 1161 the condition fc > 2 is 
fulfilled. As a consequence, there is only a small region around A F, where it is not clear 
whether truncation schemes are justified. Moreover, F is not the only possibility for the cutoff 
scale A. Each resolvent has its own cutoff scale and, for n = 1, the eigenvalue h alternates 
from spin/orbital-excitation energies to charge excitation energies for successive resolvents. 
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Therefore, for A ^ F, we expect further suppression factors ^ A/A ^ F/A 1 in many terms 
of the RG equations, provided that 

r < Z\ (323) 

for some other physical cutoff scale A. So also for n = 1 there will be many situations where 
the weak-coupling regime can be defined in a well-controlled way but a detailed study for all 
cases remains to be done, see e.g. Ref. [66] . 

Generic procedure. In the weak coupling limit a systematic procedure to solve the RG 
equations of RTRG-FS has been developed in Ref. f¥3] for the case n — 2, i.e. for models with 
spin or orbital fluctuations (similiar procedures for models with charge fluctuations, i.e. n = 1, 

are in progress [661). For n — 2, the idea is to deflne flrst a reference solution fo'" the vertex 
Gi2 from the lowest order term of the RG equation when A is larger than all physical cutoff 
scales. This means that we neglect the frequencies and the Liouvillian on the r.h.s. of the RG 
equation and truncate at the first term. This gives the poor man scaling equation, which reads 
according to (I305P 

^-^H^GW-d^l')} . (324) 

We use the Matsubara frequency representation (|29ip but have set all frequency arguments 
to zero (we implicitly assume this in the following when no argument is written). The initial 
condition for the leading order RG equation is the bare vertex. The order of magnitude of the 

leading order solution is denoted by the dimensionless parameter J ^ G^j'' • If observable is 
considered the first vertex has to be replaced by R 



dR[\] _ !_ 
dA ^ A 



{i?iyG^;! - (1^1')} . (325) 



This trivial replacement holds for all equations in the following, therefore we do not specify it 
in all cases. For an observable only the combination Tr5_Ri...„ . . . occurs for the calculation of 
the average. This has to be kept in mind because some properties of the vertex are only fulfilled 
if this combination is taken. 

The leading order solution is a good approximation when A is much larger than the maximal 
physical cutoff scale. We denote this scale by Ac, i.e. we define 

Ac ^ max{\E\,\fic.\,\hk\} ■ (326) 

Except for the Laplace variable E, no other frequencies are included in this definition because 
they do not enter into the final solution of the physical quantities. Of course, the leading order 
solution is only good for frequencies below A but since all frequency integrals are cut off by A 
in the RG equations this is not dangerous. Furthermore, we do not use any combinations of E, 
fia , and hk in this definition (like they occur in the resolvents) because each resolvent of the RG 
equation has its own combination. Temperature is also not included in the definition of Ac since 
temperature is a trivial cutoff scale. The effect of some finite temperature Tq, of reservoir a is 
just that the Matsubara frequencies of this reservoir do no longer occur in the RG equations 
for A < Ta- Since the scale of the relaxation and dephasing rates is given by Fk ^ AcJc, with 
Jc being the scale of the vertex at A — Ac, we get Fk < Ac due to the weak coupling condition 
Jc <C 1. Therefore, also the relaxation and dephasing rates do not enter the definition p26p . 

The next step is to expand the full solution for the Liouvillian and the vertices systematically 
around the leading-order solution for A > Ac using the exact RG equations. This leads to an 
expansion of the form 



Ls{E,u) = 4°) + 4^(i?,c.) + Lf{E 
Gi2{E,u,LUi,iU2) = G^a^ + G'f^iE,co,coi,LU2) + ... , 

— — — ("3") 

Gi234{E,uj,uji . . .UJ4) = Gy^4^{E,uj,uJi . . .UJ4) + Gy^^{E,uj,uJi . . .UJ4) + ••• , 

etc. (327) 
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where L^f \ G^^') - and = [Hg, ■] is the bare LiouviUian. The expansion for the vertex 
Gi...n starts at A: = n/2. As shown in Ref. [45] it turns out that this expansion is not a pure 
power series in J but can also contain terms ~ J^{\n.J)^^^. However, for J ^ 1, the series is 
also well-defined in this case. We stop this procedure at scale A = Ac and obtain a series of the 
LiouviUian and the vertices in powers of the leading-order solution at the same scale A^. The 
order of magnitude of the leading-order solution at Ac is denoted by Jc <C 1 which is assumed 
to be a small quantity for the weak-coupling regime. 

If all temperatures Tq are larger than Ac we can stop the RG at the minimum of all tem- 
peratures. In all other cases we solve the exact RG equations for A < Ac perturbatively in Jc 
using the LiouviUian and the vertices from (|327p at scale Ac as initial condition. This pertur- 
bation theory is well-defined because, looking at the exact RG equation (|315p for n = 2, we 
see that only terms ~ {Jc)'' {ln{Ac/ A))''~^ can be produced, with A being some other physical 
low-energy cutoff scale of the form (13171) . The lowest possible value for A is the minimum F of 
all physical relaxation and dephasing rates. If we assume that the order of magnitude of F is 
identical to its value Fc at scale A — Ac, we can estimate the maximal value of the logarithm 
as 

In(^) ^ ln(J,) , (328) 

where we have used Fc ^ Ac{Jc)'^ due to a dimensional analysis (except for A, there is no 
other physical energy scale available for A > Ac). Therefore, similiar to the series p27|) . we get 
correction terms ~ (Jc)'' (In Jc)'"~^ which are not dangerous for Jc ^ 1. If the cutoff scale is not 
F, some deviation from a resonance position occurs for Z\, leading to interesting logarithmic 
enhancements at resonance positions. This has been evaluated in Ref. [43j within a consistent 
2-loop formalism for the anisotropic Kondo model in a finite magnetic field. 

1-loop treatment. Here, we summarize the procedure of Ref. [43] in 1-loop, i.e. we will 
consider all terms which are proportional to either Jc, Jc or ln( Jc). Thereby, it will turn out 
that the real part of the eigenvalue of the effective LiouviUian contains terms of order Jc and 
Jc In(Jc), so that corrections of the order J^ can be neglected. For the imaginary part of the 
eigenvalue and for the current kernel, the series starts at J^. For the rates, terms proportional 
to Jc ln( Jc) can be calculated consistently only within a 2-loop analysis, see Ref. [43^ for more 
details. 

First, we note that the leading-order solutions G^^y and r']^^, do not generate any new 
structure in Liouville space but leave the form ((77)) of the original vertex and the form (|129[) 

of the original observable invariant (the latter holds only for Trs^^^^,). As shown in Appendix 
E we get 

[511',-]- , (329) 

'^[rii',-]+ , (330) 
with 

{512 521' - (1^ 1')} , (331) 

{'"12 521' + .912 ?'2i' - (1 ^ 1')} ■ (332) 

These are the so-called poor man scaling equations. As a consequence (see Appendix E), the 
adjoint of the leading-order vertices is given by 

(511')^ = 5i'i , iG[\l)^ - G« , (333) 
(niO^ = rvi , (i?!!!)^ = ■ (334) 

Note the difference to the property (|55]) . where the c-transformation has been considered. 



R 



(1) _ 
11' ~ 



p 

E(i?(i))??, = 



dA ~ A 
driv _ j_ 
dA ~ A 
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Furthermore, using (|329p . we define for arbitrary A 

g[\1 =y2piG^'^)7r = [511',-]+ • (335) 



p 



As is shown in Appendix E the RG equation for G^'^y reads 

'^'^11' _ 1 (Ml) Ml) Ml) M 



dA A 



1 (Ml) Ml) _ Ml) Mi)\ 

^ \'^12 '-'21' ^1'2 '^21 J 
= ^ {^12^ '^Sl' ~ ^l'2'5k''} ' (336) 

and the leading order result for the current operator can be written as 

^iv^=(^ii'G[r , cjy ^ -■^{mSai'r + V25a2i) ■ (337) 

Next, we consider the 1-loop RG equations ([303]) and pOSp for zero temperature (finite T 
just leads to a trivial cutoff of the RG flow) 

"'^"^'''"^ - Gu td.,————-^—- —^^G,, ,(338) 



dA Jq iuj2 + + iA + Ei2 — Ls{Ei2,LU2 + lu + A) 

• ^G,2^G2,, (l^l')l . (339) 



dA [ 'lA 

We have neglected the frequency dependence of the vertices and the influence of the Liouvillian 
on the RG equation for the vertex. These corrections are of the same order as the other 2-loop 
terms (|236p . (|238l) - (|240p on the r.h.s. of the exact RG equations and have been analysed in 
detail in Ref. [43] . It can be shown that they do not influence the final result up to second order 
in Jc (see a more detailed discussion of this point at the end of this section). To perform the 
frequency integration in (|338p . we use the form (coming out of the RG flow described below, 
see (1365]) ) 

Ls{E,uj) = Lf + 4'^ - {E + iuj)Z^^^ + Lf{E) , (340) 

where L^g^ and Z^^^ are terms linear in J which are defined via the differential equations (|366p . 

The dependence of Lg {E, uj) on the frequency ui has been neglected, which is a term of higher 
order. With this representation the frequency integral in p38p can easily be calculated with 
the result 

= ~^ Gu Ca{Eu Ls{E,2)) G.i , (341) 



dGi 



dA A 

where we have defined the important function 



- = i {G12G21, - (!<- 1')} , (342) 



9 A — 77 

Ca{z) = ln(-r— ) , (343) 



A - iz 
and 

The initial conditions a,i A = Aq ^ D from the integration over the symmetric part of the 
distribution function have already been derived in (|276p . (|278l) and (|280p 

Ls{E,u:)U^A. = - ^^i?G« G« + g G« {E,, + - 4))G« 

-Id G« G« - . I G« [E,, +7.^ Lf)GU , (345) 
G^r\A=Ao = G« - ^ \ (g^^,) G« - G« G^) . (346) 
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Analogous RG equations and initial conditions hold for Sb,{E) and Rw by replacing the first 
vertex on the r.h.s. of all equations by R. 

The RG equation (|342[) for the vertex Giy has the same form as the RG equation (|324p for 

the leading-order vertex G^^^, . However, the two vertices are not identical because the initial 

condition ()346p is different. Whereas G^^'l is the bare vertex initially, the vertex Gw has an 
additional second order correction ~ J^. Therefore, in the spirit of the series p27p . we write 



Gil' = G^il + g1i^, , (347) 



with 



^ = 1 {Gf ) G« -f G^^ G^, - (1 ^ 1')} , (348) 
and initial condition 

Gf,U-^|(G«G« -G^G^^)) . (349) 

Interestingly, it can be shown that this form holds not only initially but for all values of A, 
i.e. (|349p is the solution to the differential equation (|348p . This holds also for the case of the 
observable if Trg is applied from the left. For the proof we refer to Appendix E. As one can 
see, the correction has an additional factor i, so that it leads to the lowest order result for the 
complex part of the vertex. 

In contrast to all the other second order corrections, which have already been neglected in 
the approximate form p42p of the RG equation for the vertex, the correction G\y is important 
to calculate even the lowest order result in J. To see this we insert the expansion (|347p into the 
RG equation (|34ip for the Liouvillian, expand w 1 — Z'^^\ and neglect all terms ^ J'^ 

on the r.h.s. leading to 

dLs{E,uj) ^ 
dA 

= -^Giy Ca{Ei2 Ls{E,2)) g'^^ +iG[\^ Ca{Ei2 Ls{E^2)) ^''^ G^§ 

- I Gf^ Ca{E^2 Ls{Ei2)) g'^^ - I G^'a' Ca{Ei2 + - Ls{Ei2)) G^^^ . (350) 

We have to be careful not to count the powers in J incorrectly. As we can already see from the 
RG equations (|324p and (|348p an addtional factor ^ on the r.h.s. can lead to one order less in 
J when compared to the l.h.s. Therefore, in the regime A > Ac, we have to check if additional 
factors ^ can occur on the r.h.s. of the RG equation for the Liouvillian. To see this we subtract 
from the function Ca{z) the asymptotic part ~ -j and define 

;; / N ^ / N , ,2/1 — iz, iz , , 

so that Ca{z) and Ca{z) are integrated by the functions 

CAiz)^-FAiz) , FAiz)^AH^)--H ^ ^ ) , (352) 

Uz) = -^^Fa{z) , F^(z) = F^(z) - |(ln(^) + 1) , (353) 
dA I —liz 

with the following asymptotic behaviours 

M(ln(2) + 0(i)2) + f(ln(^) + l) for^>>|z| ^^^^^ 



2 iz 
■ z \2\ 



for A < Izl 



yl(ln(2) + 0(f)2) foryl»|z| 



Fa{z) - \ ^ ^^^^^A_^ + 1-4) for .1 « |z| • 
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We now insert Ca{z) = -£i^A{z) + ^ into the RG equation (|350p for the Liouvilhan and 
start with the contribution from the part -£^^a{z). For the first term on the r.h.s. of pSOp we 
perform a partial integration and get 



+ *{^Giy|^4^)G^'i^+*G«^4^){^G'«| , (356) 

with z = Ei2 + iuj — Ls{Ei2). Here, we have neglected terms which differentiate with respect 
to Ls{Ei2) occuring in z, which are terms of higher order when inserting the RG equation for 
the Liouvillian. The first term on the r.h.s. of (|356p leads to the following contribution for the 
Liouvilhan 

L^^''\e,lo) ^ ~iG^^^Ta{Ei2+ioo-Ls{Ei2))G^^I . (357) 

Its initial value cancels with the second term on the r.h.s. of the initial condition for the 
Liouvillian p45p by choosing approximately 

for the initial value of the parameter A (note that Fa{z) can be approximated by ln(2)yl in 
this regime, according to ()355p ). In this way the terms ^ D vanish. We note that there is no 
term arising from the RG which cancels the fourth term on the r.h.s. of p45p which is also 
^ D. This term is real and leads to a renormalization of the real part of the eigenvalues of the 
Liouvillian. If it is nonzero the model is not well defined and the frequency dependence of the 
original vertices is important. As we will see in Sec. [5] the term vanishes for the Kondo model. 

(|357p is a contribution ^ \Ei2 + ioj — Ls{Ei2)\J'^ to the Liouvillian for A \Ei2 + iut — 
Ls{Ei2)\- Therefore, it is in fact a second order contribution to Ls{E,lo). The last two terms 
on the r.h.s. of (|356p can be estimated by inserting the RG equation p24p for the leading-order 
vertex, which gives -jaG^^^ ~ ■ Using Fa{z) '--^ yl for yl ^ z, we get terms ~ which are 
of third order and are neglected. If we perform the same analysis for the other terms on the 
r.h.s. of p50p . one of the vertices in (|357p is replaced by G*-^-* or we get an additional factor 
Z^^\ Thus, also these terms are of third order in J and are neglected. 

The remaining terms in the RG equation stem from the par t ^ of Ca{z), giving rise to 
contributions to the Liouvillian in first and second order. Using p44p . we can systematically 
expand up to third order in J on the r.h.s. of ()350p . and can group the various terms by the 
following differential equations 

dL'^s\E,Oj) 1 Ml) IP , ■ , r(0)N^(l) /orq\ 

- ^ ^12 (^1-2 + - ^5 ) Ggj , (359) 

'^^^f'"^ = i { + - ^^s' ) + iE,2 + - ) g(^) } , (360) 

'^^"^f'"^ = ^ G« {-4^ - iE,2 + ^^)Z(l) + + 4°)Z(1)} G^ , (361) 

where the initial condition for l'^^ {E, lu) is zero, and for L'^g''\E,uj) and L^g''\E,u}) it IS given 
by the fifth and third term on the r.h.s. of (|345D . respectiviely, 

Lf{E,Lo)\A=A, = , (362) 
Lf\E, uj)\a=a, = " * J Gfi^ {E,y hf) G^^\ . (363) 

L^i^\E,c.)\A=Ao = (Eiv +^u;- 4''^)gW , (364) 
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There is an additional factor ^ on the r.h.s. of these RG equations and, therefore, one gets one 

order less in J when integrating them. L^^-* (S, w) is a contribution to the renormalization of 
the real part of the eigenvalues of the Liouvillian in first order in J, e.g. for the Kondo model it 
leads to a renormalization of the magnetic field. It has a linear frequency dependence and can 
be decomposed as 

L'^g\E,Lj) ^ L^s^ - {E + iu})Z^^'> (365) 
with L^g^ = L^s\e = 0, w = 0) and 

L^g'^\E^Lo) is a contribution to the renormalization of the real part of the eigenvalues in 

second order in J and can be neglected compared to ■'(£', w). It is not calculated consistently 
here because the higher order terms of the RG equations and the frequency dependence of the 
vertices give rise to terms of the same order, see the complete 2-loop analysis of Ref. [43] . 
Nevertheless, we have included it here for completeness. 

L^g^\E^Lij) is a contribution to the rate since G^^^ is complex according to the solution 
P49p . As we will see for the special example of the Kondo problem in Sec. [5] it is the only term 

giving rise to a finite current in lowest order in J. The important difference between L^^""* and 

is that the generation of is finished at yl = Ac. As we explained above, for A < Ac, 
we solve the RG equations perturbatively in Jc which is the scale of the leading order vertex at 
A — Ac- Therefore, in this regime, all quantities are systematically expanded in Jc and usual 
power counting applies. Therefore, the terms involving G^^-* on the r.h.s. of the RG equation 
pSOl) for the Liouvillian lead to corrections ^ for A < Ac and do not influence the result up 

to J2. 

In contrast, the first term on the r.h.s. of (|350p gives rise to a contribution ~ in the 

regime A < Ac. Replacing there the leading order vertices by their value G^^^,'^ at yl = Ac, and 
using (|352p . we can integrate this term from Ac to A with the result 



Ls{E,uj)a = Ls{E,uj] 



Ac 



{^a{Ei2 +tuj- Ls{Ei2)a) - TaAEi2 + - Ls{Ei2)a:)} G^i^' . (367) 



Thereby we have neglected the derivative of Ls{Ei2)a with respect to A (see the discussion 
below). The term involving Ta^ can be taken together with L^g'^\E , uj) a^ by using the result 
P57p at yl = ylc- Together with the relation ()353p we get 

LsiE,Lo)A = + Lf'{E,oj) + Lf^'{E,oj) + Lf^'{E,oj) + 

+ Lf"''^''{E,uj)- iG^I^'Ta{Ei2+iuj-Ls{E^2)a)g':^^" , (368) 

where the index c indicates always the value ai A — Ac, and L''g°'^^{E,uj) is defined by 

L^P^{E,Lu) = (369) 

= -\G[T{E,2+^u.-Ls{E,2)AJ (in ^ , , + A 

^ V -2i[Ei2 + iLU - Ls{Ei2)aJ / 

in the regime A < Ac- Setting A = and u; = 0, we get the effective Liouvillian L'g^{E). Using 
Fa=o{z) = 0, we get 

L'^J^iE) = Lf + Lf\E) + Lf^^^iE) + Lf^"{E) + L^g'^^iE) . (370) 
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The problem is the precise value of Ls{Ei2) in (|369p . The integration can not be performed 
analytically if the yl-dependence of is (£'12) a is taken into account in the regime A < Ac. Since 
no terms are generated linear in Jc in the regime A < Ac, we get according to (|344p . 

Ls{E)a = (4"^ + 4'^^ + LfHEU) , (371) 

with L):'{E)a = L):'{E,uj = 0)a which are given by the second order terms of (|368p in first 
approximation. Using (|352p we see that the last term on the r.h.s. of (|368p has a complicated 
logarithmic dependence on A which can only be integrated numerically. However, even from a 
more precise numerical analysis we do not expect that the variation of Ls{E)a will change the 
result for the effective Liouvillian significantly. The important part of Ls{E)a is the hermitian 
part which determines the position of the resonances where the logarithmic function in (j369p 
becomes maximal. According to (j37ip . this part can be cut off in good approximation by the 
terms up to linear order in Jc which are independent of A. Expanding (j37ip up to this order 
gives 

hc = Lf+Lf;^Z^P^Lf , (372) 

where the additional index d means that we take only the diagonal part with respect to the 
Liouvillian L^"-* (otherwise the various terms do not commute and we consider higher-order 
corrections which are not calculated consistently). Terms of order are neglected in he- Note 
that he can depend implicitly on E via the cutoff Ac when E is the maximal low-energy scale. 
Only in a small region of order F around the resonances the antihermitian part of Lg cuts 
of the logarithmic divergencies. The prefactor of the imaginary part occurs then under the 
logarithm and gives only a very small perturbative correction. Therefore, in performing the 

integral, the error is quite small when neglecting the /l-dependence of Lg {E)a- We have taken 
the value a.t A — Ac in (|369p . but, alternatively, one can also take the value at yl = 0. We denote 

the antihermitian part of L^g\E)A^ by —fc{E), which according to (|367p . is given by 

^ciE)^^L^^^^^iE)+^Lf/{E) , (373) 

where we take again only the diagonal part with respect to L^g'^ . In summary, we replace 
LsiEi2) in jSMl) by he - irc{Ei2). Using 

Ac 



{Ei2-hc + irc{Ei2)) In- 



-2i{Ei2-hc + trc{Ei2)) 



^Ue,,) (^12 -hc)+i\\Ei2~ hc\ , (374) 



with 



Hr{E) = E [ In , = + 1 ] , (375) 

we can decompose the final solution (j370p in hermitian and antihermitian parts (which holds 
when E is real according to the property (|333p ) 

Llff{E) = h'^^iE) ~ iV-ffiE) , (376) 
h^fi{E) = Lf + 4^' ' EZ^^^^ + Lf^^'^^E) + lP^E) , (377) 
r'ff{E) = *4''''"(£;) + ilf^^E) , (378) 



with 



4'^'«^(i?) = -^G«^ff^^(^^^)(i?i2-MG^f , (379) 

TT 

4 



^Lf^''-■{E) = ^G^^^\E,2~hc\G^i' ■ (380) 
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In these formulas L^g^'^, Z^^'>'^, L^g''^'^{E) and L^g''^'^{E) are determined by solving the RG equa- 
tions p66p . p60p and p6ip up to the scale A^.. The quantities ft-c and are defined by (|372p 
and ()373|) . ft-c is identical to h'^f^ up to the terms linear in Jc, with E — > ' and taking the 
diagonal part with respect to i^"*. As discussed above, without significant error one can also 
calculate FciE) from E^f^iE) at = 0, i.e. 



We note that a precise justification to take the values at yl = is given in Ref. [43]. Similiar 
formulas hold for the kernel Sfi{E) of the observable R if the first vertex is replaced by R in 
all equations. Note that all formulas contain an additional implicit dependence on E via the 
cutoff Ac defined by (l326l) . 

The final results ()377p and (j378p are generic formulas including all terms ^ Jc t Jc i Jc Jc) 
for a fermionic model with spin and/or orbital fluctuations. Terms ^ ln( Jc) occur at resonance 
\Ei2 — hc)\ — 0, where the logarithmic term in Hp becomes of order \n{Ac/Ec) ~ ln( Jc) since 

Fc ~ Jc^c according to (|373|) . Although the prefactor of Hp [E) becomes very small for i? ^ 0, 
derivatives of the effective Liouvillian with respect to some energy scale like voltage or magnetic 
field will show the logarithmic increase at resonance. However, we note that the precise position 
and the broadening of the resonance is not exactly given by the poles of the reduced density 
matrix, but can only be calculated numerically. Nevertheless, the difference is not very large 
and hardly visible. We note that the logarithmic part ~ J^ ln( Jc) of (j377p appears only in the 
renormalization of the hermitian part of the Liouvillian, i.e. it does not contribute to the rates, 
e.g. to the stationary transport current. However, for the time evolution of physical quantities, 
it will play a crucial role due to its logarithmic dependence on the Laplace variable E, leading 
to branch cuts in the complex plane. As we have shown, the logarithmic terms arise from the 
solution of the RG equations between Ac and some other cutoff scale A = \Ei2 — hc\. This 
is generically the case because the solution of the RG equations at scale Ac involves only the 
maximal cutoff scale Ac and no logarithmic terms ~ ln(ylc/Z\) with some smaller cutoff scale A 
can occur. Therefore, the solution at scale Ac is a power series in Jc with possible logarithmic 
contributions ~ J^ (In Jc)*^^^ from higher orders (but independent of A). 

What we have not analysed here are the 2-loop terms (|236p . (|238p - (|240p of the RG equations 
and the influence of the frequency dependence of the vertices. They provide further terms of 
order J^ on the r.h.s. of the RG equations and can, as we have seen above, influence the result 
for the Liouvillian in order J^ at scale Ac- This has been analysed in detail in Ref. [13] with the 
result that all these terms nearly cancel, the only effect is a small correction for the definition 
of the function H, given by (|375p . which has to be replaced by 



i.e. only the factor 2 in the argument of the logarithm is absent. This leads only to a small 
correction of order ~ J^ for the hermitian part of the Liouvillian. Neither the terms ^ Jc or 
~ Jc ln( Jc) for the hermitian part nor the terms ^ Jc for the antihermitian part are changed. 

However, the 2-loop analysis of Ref. [43] is very important if one is interested in corrections 
to the rates of order Jc \n[Ac/A) which become Jc ln( Jc) at resonance. Such terms lead e.g. 
to experimentally accessible resonances for the differential conductance when the bias voltage 
matches with certain energy excitations like magnetic fields. It turns out that especially the 
frequency dependence of the vertices is important for the generation of such terms below ylc. 
Moreover, for A > A^ even the frequency independent part of the vertices is corrected by various 
terms in second order in J. For the Kondo model, it can be shown that these corrections lead 
to a redefinition of the Kondo temperature. However, for more generic models with orbital 
fluctuations, it might be the case that the vertices in second order become a new structure in 
Liouville space and will also influence the prefactor of some Jg \n.{Ac/A) terms. 




(381) 




(382) 
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4.5 Strong coupling limit 

When the conditions of weak coupling are not fulfilled no systematic truncation scheme can 
be applied to solve the RG equations. Nevertheless, since the RG equations are a set of fully 
self-consistent equations with the renormalized Liouvillian appearing in the denominator of all 
resolvents, there is some hope that the equations can also lead to reliable results when the 
coupling constants become of order one. A divergence of the coupling constant like in poor man 
scaling equations is not expected since the relaxation and dephasing rates will also increase for 
increasing coupling and will cut of any divergence. However, whether quantitatively reliable 
results can be expected is not at all clear and certainly an interesting field for future research. 

Using RTRG-FS, the problem of solving the full RG equations numerically in any truncation 
scheme lies in the frequency dependence of the vertices and the Liouvillian. Therefore, to start 
with a minimal ansatz, one possibility is to retain only the dependence on the Laplace variable 
E since this variable takes discrete values shifted by the chemical potentials of the reservoirs, 
see piOp . So e.g. for the fermionic problem of spin/orbital fiuctuations, where the RG equations 
are given by (|303p - ()308p . a minimal approach would consist in omitting the dependence on all 
other frequencies and take only the first term on the r.h.s. of the RG equations into account. 
At zero temperature, this leads to the equations 

- ^Mg) = ^G,2{E) tduj, ^ ^ ^^^Gr,{E,2) , (383) 

dA Jo Ei2+iA + iuj2- Ls[Ei2) 

_dGn_CB) ^ \ , G2y{E^2) - (l-n . (384) 

dA Ei2 + lA- Ls[Ei2) 

Performing the integral leads to 

dLs{E) 



dA 



iGi2{E)CA{Ei2-Ls{Ei2))G2-i{Ei2) , (385) 



= , -J, ]-r .r, , G2AEl2) " (1 - l') , (386) 

dA A~ iEi2 + iLs[Ei2) 

where Ca{z) ~ ln( '^'^ZH ) ti^s been defined in (|343p . The initial conditions are given by (|345p 
and (jMS]) . 

At finite temperature a similiar set of equations can be set up, the only difference is that 
integrals have to be replaced by sums over Matsubara frequencies, and we have to use the 
smeared theta function Qt{!-^) defined in (|203p . Instead of (|383p and (|384p . we obtain the 
equations (for simplicity, we assume that the temperatures of all reservoirs are the same) 

dA \ Ei2+lAT + lUJn- Ls{Ei2) j ' 

_dGn_(^ = -»gl2(i^) r, ^-^ r ir, , G2,,{E^2) (l-l') , (388) 

dA Ei2 + iAt — Ls(Ei2) 

where uj„ are the discrete Matsubara frequencies and At is the Matsubara frequency lying 
closest to A. The sum in (|387p can be performed analytically with the result 

"^^ff^ = -tGi2(E)£A{Ei2~Ls{Ei2))G2i{Ei2) , (389) 

= Gi2{E)- \ . G.AEu) - (1-1') , (390) 

dA At — iEi2 + iLs[Ei2) 

i.e. the same result as for T — but A is replaced by At in the RG equation for the vertex, 
and the function Ca{z) has to be replaced by the more general expression 



(391) 
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where '^'(z) is the Digamma function with asymptotic properties given by pi3p and (|214p . 
Furthermore, it can be shown that the initial conditions for the LiouviUian and the vertices are 
the same as those for zero temperature. 

These sets of RG equations will be analysed in Sec. 15.31 for the isotropic Kondo model in 
the absence of a magnetic field. As expected the coupling constants stay finite and a quite 
promising agreement with results from NRG calculations in equilibrium is obtained. However, 
it turns out that the equations are unstable against exponentially small changes for the inital 
condition for the relaxation rate F. Therefore, it has to be analysed in future how the equations 
can be improved by higher order terms in order to stabilize them. One possibility is to consider 
the frequency dependence of the vertices in leading order by taking the frequency independent 
vertices on the r.h.s. of the RG equations, together with the consideration of the next order 
terms of the RG equations. A numerical analysis of the full frequency dependence is certainly 
very time consuming but it has to be expected that the numerical solution will be very stable 
since all the imaginary parts of the denominators are positive. All these investigations will be 
the subject of future research. 

5 Application: The nonequilibrium Kondo model 

In this section we will apply the renormalization group formalism developed in Sec. [3] to the 
nonequilibrium isotropic Kondo model in the absence of a magnetic field. For weak coupling, 
we evaluate the general equations of Sec. l4.4l to get all physical quantities up to second order in 
the coupling. The results presented are a special case of the more general treatment in Ref. [43] 
where the anisotropic Kondo model in the presence of a magnetic field has been considered, 
and all quantities were calculated up to third order in the coupling. In Sec. 15.31 we will discuss 
some preliminary results in the strong coupling regime. 

5.1 Model and algebra of basis operators in Liouville space 

Model. We discuss here the Kondo model introduced in Sec. 12.21 The quantum system is a 
spin i which is coupled to the spins of several reservoirs by the exchange coupling term 

V = ^911' ■ aiav ■ , (392) 

with 

_ 1 / Jc,a'S-a^^, forr? = -V = + 
^" 2 {-Jc'cS-a,,, forry = -V = - ' ^^^^^ 

Jaa' = Ja'a are the exchange coupling constants, S_ is the spin i operator of the quantum 
system, and a are the Pauli matrices. We have taken here the isotropic case where the exchange 
couplings are the same for all spatial directions, a is the reservoir index. In principle the 
formalism presented here is applicable to an arbitrary number of reservoirs kept at chemical 
potentials /ia, but sometimes we refer to the case of two reservoirs (a = ± = L, R) with 
chemical potentials given by 

Ma = « y , (394) 

where V is the applied voltage and we use units e = h = 1. 

We assume here the case of zero magnetic field, i.e. the Hamiltonian and the LiouviUian of 
the quantum system are initially set to zero 

Hs ^ , Lf^^O . (395) 
Using dZZl) we get for the vertex in Liouville space 

Gff, = i ( -^r' -TP '^ = -'^[ = + , (396) 

" 2 \ -Ja'aL!' ■ CT^v for ry = -ry' - ' ^ ' 



Will be inserted by the editor 



67 



where LP are two fundamental spin vector operators in Liouvillc space defined by {A is an 
arbitrary operator in Hilbert space) 

L+A^SA , L-A^~AS . (397) 

For convenience, we will consider in the following always the case rj = —rj = + when 
discussing some vertex Gn/, since the antisymmetry G\\i = — Gi'i gives trivially the other 
case 77 = —rj' = — . 

For J] — —T]' = + and using p47p . the vertex of the current operator in Liouville space can 
be written as 

with 



{i-^f^' = cl^'PGfr , (398) 



= -7t('^«7-'5«'7) - . (399) 



1 
2 

The kernel (I137[) for the current is denoted by S-, (E) . 

Since the Kondo model emerges after integrating out the charge degrees of freedom, a high- 
frequency cutoff is needed, which we describe by the cutoff function ()47|) 

P(^) ^ 2^'n2 ' (400) 

where D is the physical bandwidth (which is fixed here and is not changed during the RG 
procedure). 

Basis operators in Liouville space. Due to spin rotational symmetry not all 16 matrix 
elements of each operator in Liouville space are needed. To get a minimal set for the isotropic 
case, we first define two scalar operators by 

L"" ^ - + L+ -L- , = - ~ L+ -L- , (401) 

4 — — ' 4 — — ' 

and three vector operators by 

1/ ^ ^{L+ - LT - 2iL+ , (402) 

- -liL+ + L-) , (403) 

= ^{L+ - LT + 2iL+AL") . (404) 

Using spin rotational invariance together with spin conservation (i.e. the parity of s— s' and s— s' 
must be the same for any matrix element {Ls)ss' ,ss') and the property Tr^ Ls{E) — Trg L" — 0, 
the Liouvillian Ls{E) and the kernel S^{E) can be represented as 

Ls{E) = {h{E) - ir{E))L'' , (405) 
TvsS-~f{E) = ir^{E)Trs L'' , (406) 

where the trace over the quantum system has been written in the second equation since only 
this combination occurs for the observables. Many equations for the observable are only valid 
in this sense. Therefore, we will use in the following the convention that each time the current 
vertex or the current kernel occurs, we act implicitly with the trace over the quantum system 
from the left. The kernel for the current is written such that r'-y(0+) gives the stationary current. 
Using Trs L'> A = Trs A (see dHH) below) and Trg ps{E) = i/E (see (fTTT])). we get from (fTSe]) 

and nng) 

{r)iE) = ^r,{E) , {Py' = A(0+) • (407) 

The quantities h{E) and r{E) correspond to the real and negative imaginary parts of the 
eigenvalues of Ls{E), respectively. 
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Table 1. Algebra of Liouville basis operators. The table shows the product . The same comes 

out for the combination ^(L^)"^ (L^ but the sign of the operators , , and has to be 

changed in the table. 



To express the vertices Gu' = G+aa,-a'a' (note that we use implicitly rj — ~r( ~ +) in 
terms of basis operators, we need an appropriate representation for the reservoir spin labels u 
and cr'. We introduce the tensor operators 

{V^\„, = b,,: for X = a, , (408) 
(V<)..' ^ ■ o^„, forx = 1,2,3 . (409) 

The operators act simultaneously in Liouville space of the quantum system and in the space 
of the reservoir spin labels. We note that the notation for L°- and i'' is ambigious, however, it 
should be always clear from the context whether and also act in reservoir spin space or 
not (e.g. in (j405p and (|406p this is not the case). The vertices are then represented as (implicitly 

for r\ = —rf = +) 

Gn'{E)= Gl^'iE){L^U' , (410) 

X=a,6,l,2,3 

I'li'iE) - Ili{E){L^U' • (411) 

X=a,64,2,3 

In this way we only need to solve RG equations for the c-numbers G^^,{E) and l2a'iE) (or 
for the matrices G^{E) and I'^^{E) in reservoir space). The algebra of these representations is 
closed as we will see in the following. 

We note some important transformations and properties of the basis operators which will be 
frequently needed in the following. We define the transpose by only interchanging the reservoir 
spin indices 

{{L^f)^^, - {L^).'. . (412) 

As a consequence, using Gn' — — Gi'i, we get for the representation of the vertex for all cases 
of rj and t]' 

Gn'{E) = J2 {- &(J^t\lx)T fo' ^ = = t (413) 

X— a, 1,2,3 

where G^{E) is considered as a matrix in the reservoir indices. 

In contrast, the c-transform {L^Y, defined in (pT|) for operators in Liouville space, is defined 
only with respect to the degrees of freedom of the local quantum system. We get 

{L^Y = for x = a, 6, 1,3 , (L^f ^ . (414) 

Applying this transformation to the representation of the Liouvillian, the kernel, and the vertices 
in terms of the basis operators, and using the properties (|262l) - p64p . we obtain the following 
helpful relations (note that h{E) and r{E) are defined as being real for arbitrary E) 

h{E) = -h{-E*) , r{E) = r{-E*) , r^{Ey = r^{~E*) , 

GIAET = -Gl,^{-E*) , il^Ey ^ -PJ^{-E*) forx = a,6,l,3 , 
GlAET - Gl,^{-E*) , iAAEY = il?^{~E*) . (415) 
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Denoting by Tro- the trace with respect to reservoir spin indices, we get 

TraL^^2L^ torx = a,b , Tr^ i'^ = forx = 1,2,3 . (416) 
If the trace over the quantum system acts left to the basis operators, we obtain 

TrgLX = for x = a, 2, 3 , Trg ^ Trg A . (417) 

Table [1] shows the closed algebra L^L^ of the basis operators L^. The same algebra holds for 
the combination 

[{L^fiL^'fY , (418) 

but a different sign occurs for the operators L^, L^, and L^. Note that one is not allowed to 
write for (j418[) because the transpose is only acting in reservoir spin space but not in 

Liouville space of the quantum system. 

To evaluate the RG equations in lowest order, we will either encounter terms of the form 
A-12B211 (where we sum over the index 2) or K{Eiii)Aiii Byi (where we sum over the indices 1 
and 1'). Here, An' and Bn' are arbitrary vertices represented in the form (|410p . and K{Eiii) 
is an arbitrary function of the variable Eiy = E + fiiy = E + -qjia + r/'/ia' introduced in (|247p . 
Using the antisymmetry Aiy = —Ayi and Biy — —Byi, and inserting the representation 
(|413p . one obtains after summing over the two possibilities rj — —r]' = ± the helpful identities 



A^B^L^L^ for 77 = -77' = + 

" l(Ax)^(Bx')T(ix)T(ix')T ioTV^-v'^- ' ^ ' 

K{E^y)AivBi,i = 2K{Eo,^,)Al^,Bl,^Tr^L^L^' , (420) 

where we use the definition 

Eaa- ^ E + ^la ~ Ha> . (421) 

A^ and B^ in (|419p arc matrices in reservoir space and we sum implicitly over all x, a, a' . 
Using (|416p and the algebra of the basis operators, we see that (|420p is only nonzero for the 
combinations 

ixx') = (aa),(fefe),(13),(22),(31) . (422) 
For the evaluation of the RG equations in strong coupling we will also need the identity 

K[Ei2)Ai2B2v - \ (Tx\T(Tx'\T f , ■ <.'^23) 

i K{Ea^a) Al^^ B^,^^ {L^Y {L^ y for r/ = -rf = - 



5.2 RG in weak coupling 

In this section we apply the general scheme described in Sec. 14.41 to the Kondo problem. Since 
no magnetic field is assumed, the parameter Ac defined in ()326ll is given by 

Ac = max{|£;|,^„} , (424) 

i.e. Ac is just the maximum of the Laplace variable E and the voltage V. To stay in the weak 
coupling regime, we assume that this energy scale is much larger than the Kondo temperature 
Tk which is the energy scale where the RG equations for the leading order vertices diverge (see 
(p2l) below) 

max{\E\,V} Tk = Aoe'^ . (425) 

Here, Aq is the initial value for A and Jq is the overall scale of the initial exchange couplings. 

Leading order RG. We start with the RG for the leading order vertices G^^!, G^^ and 
I^^P , defined by the RG equations p24p and (|336p , and the result p37p for the current vertex. 
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Using (|396p and (|402p - (|404p . we make the same ansatz as for the initial form {rj = —rj' = + is 
always assumed) 

g[\]^~JL^ , GW=ij(Li + L3) , /7« = ij7ii , (426) 

where, in the equation for the current vertex, we assume implicitly that Tr^ is taken from the 
left which cancels the contribution of according to (|417p . J = is the symmetric matrix 
of the exchange couplings in reservoir space and the antisymmetric matrix = —{J^Y' is 
defined by 

Jlo.' = ol^' • (427) 

To prove (I426P we use this ansatz together with (I419P and the algebra of the basis operators 
to evaluate the r.h.s. of the RG equations as 

= -^J' {{L' + L^)L' {{LYiL' + LYf)} = -\j\L' + L^) , (429) 

We see that (|426p is fulfilled provided that the matrix J is defined by the well known poor man 
scaling equation of the Kondo problem 

The current follows directly from (j337p . The poor man scaling equation can be easily solved 
with the ansatz (if this is also fulfilled initially) 

a 

with J being the overall scale of all exchange couplings given by 

^ - J ^— (432) 

dA~ A ' ^ - 21n^ ' ^ 

where Tk is the Kondo temperature given by (|425p . As a consequence, under the condition 
(|425p . the leading order vertices will stay small in the regime A > Ac- 
Vertices in second order. The vertices G^^, and I^^f^ follow directly from (|349p by 
inserting (|426p and using (|419p 

Gfi^' - ^^^{-l)J'{L\L^+L')-{{LYiL'+Lrf)} = ^^J'L^ , (433) 

I'A'' - '^^l{J^JL\L^+L^)-JJ-{{LY{L'+Lrr)} 

= -i^iJ^J - JJ"')L'' . (434) 
8 

As we will see below, O'f-^, is essential for the leading order result for the current, whereas /^r'^'' 
is not important in this order. 

Liouvillian and current kernel. To evaluate the Liouvillian L'^J^ {E) and the current 
kernel S.y{E) from (i57Sl) - ([5751) . we need the identity (ji^ for B = G^^! = -JL^ or B = 
G^p = i^j'^L^. Due to (|422p . the corresponding possibilities for A are A — G^^' = —JL'^ or 
A = I^^P ^ (l/2)JTii with the resuh 

K{En')G^^^,G^^,{ = 2K{Ec,a'){Jaa'fL'' , (435) 
KiEn')I^iPG^^{ = 37rzif(£;„„0cL"^-'(-^')-'^' • (436) 
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As a consequence, together with Lg = 0, we find directly from p59p . p60p . p79p and 
that 

■* and Z^^^ follow from p66p . For we need K ^ fj.^ — ^la' in (I435P which gives zero 
due to a ^ a'. As a consequence we get also he = ^ from (|372p . For Z*^^^ we need if 1 in 
(|435ll . This gives 

^Z(i) - -(Tr„j2)L'' = ^(Tr.J)i'' , (438) 

where Tr^ is the trace with respect to the reservoir indices, and we have used the poor man 
scaling equation ()430p in the last step. As a result we find 

L^s^ = , = (Tr„J)L° . (439) 

To be precise one should subtract from Z^^^ a contribution where J is replaced by its initial 
value Jq. However, this contribution is very small and vanishes in the scaling limit Jq — > 0, 
_D — > cx), such that the Kondo temperature remains constant. 

L^s''\e) and S^^^''\e) follow from We can replace Z^^^ Tr^J in this equation 

since L°'L^ — L/^L"" ~ L^. Thus, we get u!f''\E) ^ since no combination applies according 

to (|422)) . For L^g^^E), we get the form p35)) with K E^a' ■ The part ~ f^a' does not 
contribute due to a ^ a'. This gives 

^L^s'\e) = -jE{Tr^J){Tr^j')L^ ^ }.E-^{TT^jf , (440) 
where we again used the poor man scaling equation (|430p in the last step. So we get 

= ^EiTi^jfL'^ , 42c) ^ _ ^44^^ 

L^g'^\E) is a contribution of second order in J with the same form as the linear order term 

-£;Z(i) = -E{TiaJ)L'' for L^s\e) from Therefore, it is an unimportant contribution. 

The most important terms arise from E':y^\E), "■'^'^(i?) and l'^^"'^^^ [E), which can be 

calculated from (|5SD)) . (15751) and ([5511)) . For S^^''\e), the combination (|i5S)) with X ^ E^a' 

applies. We see that only the vertex c'^^y contributes but not /^i^''. The part from E gives zero 
due to a ^ a'. It remains the part 

-^S^^'\E) = ^37r*(^„-M.')cL'-/-'(J')a.'i'' = -*Y('^"-/^-')cL'^(-^-')'i' ■ 

(442) 

Inserting c^^, — — {l/2){5a'y — 5a' -y) and interchanging a ^ q' in the second term, we get with 
s'^^^\e) = ir^{E) (see dM])) the result 

r^{E) = ^5a^{tMa~ ^J^c.'){Jc.a'f . (443) 

For A — Aq, this result agrees with the initial condition (1363]) . L^s^^^^E) and l'^^^^^^E) 
follow directly from the combination (|435p with K Hp {Eaa') and K — ^ Eaa', respectively, 
with the result 

Lf^^^E) ^ ^Hp^^^^^,^iEaa')iraa')'L'^ , (444) 



72 



Will be inserted by the editor 



where Hr{E) is defined in (|382p (differing from (|375p by an unimportant factor 2 inside tlie 
logaritlim arising from a full 2-loop analysis, see the discussion at the end of Sec. 14. 4|) . Using 
([5gT|l . fc can be calculated from L^f°^'>^"' {E) = -irc{E)L'' (due to L^L" ^ L" the operator L" 
is not needed for /^(i?)). 

Summary. Collecting all terms and inserting in (|376[) - p78p and (|407p . we obtain the fol- 
lowing result for the effective Liouvillian and the current kernel 

L^ffiE) = {h''ff{E) - ir'ff{E))L'' , (446) 
TysS-,{E) = ir^if{E)TYsL'' , (447) 



with 



h^ff{E) = -EiTl-^r) + ^E{TT^rf - i?f^(£;^,„,,)(i?oc.')(^aa')' , (448) 



TT 



r^ff^E) = - \E^o.'\{JL'y « r,{E) , (449) 
r:;ffiE) = iniE) = ^ <5„, (^„ - (J^„,)' , (450) 
where Eaa' = E + — ^^a' and 

ME) - E L'^^^^m + l\ . (451) 

is the solution of the poor man scaling equation (I430p at A = Ac = max{|i?|, V}, given by 
(PT|) and (11221) as 

•^aa' = 2 y/x^^lF^ J , = 1 , J = — jnax{|_E|,y} ' (452) 

a 2 in 5^^^ 

Note that J'^ depends imphcitly on E for \E\ > V. 

If only two reservoirs with fia — aV/2 are present, the result can be written as 



h^ff{E) = -E{Jl + J%) + ii?(J£ + j; 



cn2 

-H^^^^^{E){{Jl)^ + {r^f) - H^^^^^y^{E±V){J^,f , (453) 
r^ff^E) ^ I (( J£)2 + (J]^)2) + I |i? ± y| (j;;,)2 « f,(iJ) , (454) 

qffiE) = (/'')(£;) = 7-(J,';.)'^ , (455) 

where Ja = Jaa are the diagonal and Jnd — Jlr. — Jrl the nondiagonal exchange couplings, 
and we sum implicitly over the two possibilities ± on the r.h.s. of the equations. 

For i? = 0, (|454p and (|455p are the well-known golden rule results for the Korringa spin 
relaxation rate and the current, with the exchange couplings replaced by the renormalized 
ones from the poor man scaling equation cut off at the voltage. However, for E 0, interesting 
logarithmic contributions - E\n{V/\E\) (for £; < T/) or ~ {E±V) \n{V/\E±V\) (for E ~ tV) 
appear in second order in Jc for h'^-^^ (E) , which are present in both regimes E > V and E < V. 
This will generically lead to branch cuts in the complex plane for the reduced density matrix. 
In contrast, up to second order in Jc, the current depends on the Laplace space E only for 
E > V via J'^ — 1/(2 ln(i?/TR-)). However, as has been discussed in detail in Ref. 43J, this 
changes in third order in Jc, where also the current gets similiar logarithmic contributions as 
h'^ff{E). An interesting field of future research will be to discuss in detail the consequences for 
the time evolution of the spin population and the nonequilibrium current. 
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5.3 RG in strong coupling 

In this section we briefly discuss some preliminary results in the strong coupling regime. The 
strong coupling regime is defined by the condition that all low energy scales fall below the 
Kondo temperature 

\E\,V,T<Tk . (456) 

In this regime no systematic truncation scheme can be applied because the vertices become 
of order one. Nevertheless, it can be studied what result the RG equations in the minimal 
approximation scheme proposed in Sec. 14.51 give. We will study here the stationary case (i.e. 
E = Q) for the differential conductance 



GiT,V) = -^{rr = -^r^{E = 0) (457) 



for the special case of two reservoirs with Ha = aV/2. Note that we work in units e = ft = 1, so 
that finally the result has to be multiplied by tt to get the conductance in units of the universal 
conductance Gq = 2e^//i. From exact solutions (see e.g. Refs. |31I47I67) ) it is known that at 
T = V = 0, the conductance should be universal and the following Fermi liquid result should 
apply for low temperatures and voltages 

G/Go = l-^~^ , (458) 

where the Kondo temperature Tr- in this result is not universal. However, Tk drops out when 
considering the universal ratio of the coefficients of the terms ^ and ^ . Furthermore, for 
= 0, one can compare the linear conductance with results from numerical renormalization 
group, where the Kondo temperature is defined such that the conductance at T = Tk is given by 
G = Go/2 [68]. With this numerical method one can especially study the crossover at T ~ Tk 
very reliably. 

We restrict ourselves to the leading order form of the Liouvillian, the kernel, and the vertices 
as they have been used for the weak coupling regime 

Ls{E) = {h{E) - ir{E))L- , (459) 

Tts S-y{E) ^ ir-y{E)TTs L'' , (460) 

Gii'(£;) = G^{E)L^ + G^{E)L^ , (461) 

TTsIJ,,iE) = f\E)TTsL' . (462) 

However, we include the full dependence on the energy variable E since this is not very time 
consuming for a numerical solution. The vertex G^{E) has been included because, as we have 
seen within the weak coupling analysis in Sec. 15. 2[ the vertex G^^, L'^ is essential to obtain 

the current in leading order, whereas I^j?^ ~ does not contribute in this order, see the 
discussion before (|442p . For convenience we introduce the matrices J and K defined by 

&{E)^-J{E) , &{E)^iTTK{E) , (463) 

and take only the real part of Jaa'{E) and Kaa'{E) into account. For simplicity we calculate 
only the stationary current so that we can set i? = for ^ and F^. With these approximations 
we obtain the following relations from (|415p for E being real 

h{E) = -h{-E) , r{E) = r{-E) , r; = , 

Jaa'(^^) = Jc.'c,{-E) , K^^,{E) = K^,^{~E) , III = -IlL ■ (464) 
Due to the property — ' '^^'^ parametrize the current vertex by 

ill = \clo..Ji ■ (465) 
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Finally, we use the abbreviations 

C{E) ^ CA{E~h{E)+ir{E)) , (467) 

where the function Ca{E) is defined in (|39f P for the general case of finite temperature. At is 
as usual the Matsubara frequency lying closest to A. From (|464p we obtain 

n{E) = n{-E) , C{E) = C{-E)* . (468) 

The vertices are parametrized in such a form that the initial conditions become very simple. 
From (|426| and (|435)) we obtain 

i^aa' = ^(^^)aa' , J I = Jnd , (469) 

and Jaa' are identical to the initial exchange couplings. From the general expression for the 
initial condition (|345p for the Liouvillian and the current kernel, and using the helpful identity 
(I420p . we get with the original form (|426p for the vertices the following result for the initial 
conditions of h{E), r{E), and F^iE) 

h{E) = £;-(Tr„j2) , r{E) = D^(Tr„j2) , r,{E) = jV^{J,,d? • (470) 
16 8 4 

The initial conditions for h{E) and r^{E) are unimportant since they vanish in the scaling 
limit J — > and D oo such that the Kondo temperature remains constant. In contrast, as 
already discussed previously, the terms ~ D for r{E) are artificial and the initial value Aq for 
the parameter A has to be chosen such that these terms vanish. In Sec. 14.41 we found for the 
weak coupling case that ylo is approximately given by 

^°-T6?^^ ' (4^1) 

leading to 

r(i;) 21n(2)ylo(Tr„j2) . (472) 

Whether this is also the correct value for the strong coupling case where the RG equations are 
different is not at all clear. In fact we will see below that the strong coupling solution depends 
crucially on the inital value of r{E). 

We study now the RG equations (|389p and (|390p proposed in Sec. l4.5l as a minimal set to go 
beyond a weak coupling analysis. Using the relations (|420p and (|423p together with the algebra 
of the basis operators in Liouville space, we obtain after some straightforward manipulations 
(similiar to those already explained in all detail in Sec. 15. 2p the following RG equations 

-^h{E) = 2Im{£(i;„„0} Jaa'{E)Ja'aiEc.c.') , (473) 

■^r{E) ^ 2Rc{C{E^^,)} J^^,iE)J^,^iE^o.') , (474) 

+ E[{Ea2a') Ja2a'{E) J 0,02 {Ea2a')} (475) 

for the determination of the values for h{E) and r{E), and the RG equations 

-^F-^ = -l2T:iL{fla- ^Jia')iZa' Ka'a{^J^a- ^^.a') , (476) 
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■77(^^2 - tJ'a')ll^a' Jaa2itJ'a2 ^ A^a')} (477) 



-n{Ea2a') J a2a' {E) K aazi^aza')} (478) 



to obtain the stationary current rate F^. Note that h{E) and /^(i?) enter into the resolvent 
n{E) and the function L{E) via the definitions (|466p and (|467p . In this way, especially r{E) 
provides an important cutoff parameter for the RG flow in the strong coupling limit so that 
the coupling constants do not diverge and a finite conductance comes out. 
From (|475p one can prove the relation 

Jcc'iE) = Ja'a{Eaa') ■ (479) 

We now turn to the case of two reservoirs with fia ~ aV/2. Inserting the parametrization 
(|465p and using the symmetry relations (I464p , one obtains after some algebra the simplified RG 
equations 

-^r^ ^ -&Ti^lTnC{V)JiKRL{V) , (480) 

47J1 = -niV)JiiJL + JB)iV) . (481) 
aA 

Using the symmetry relation (|479p . the RG equations for h{E), r{E), Ja{E) — Jaa{E), and 
Jaa{E), with a = —a, can be simplified to 

47HE) = 2ImC(E) JJEf + 2lmC(E + aV) Ja,a(Ef , (482) 
dA 

-^riE) = 2ReL{E) JJEf + 2ReL{E + aV) Ja&{Ef , (483) 
dA 

4tJc,{E) = -n{E)J^{Ef ~ n{E + aV)Jo,a{Ef , (484) 
dA 

■^Jo.&{E) = -n{E)^{JL + JB.)iE)Jo.dE) - 

^n{E + aV)]^{JL + JR){E + aV)Jo,o.{E) . (485) 
The RG equations for Ka{E) — Kaa{E) and Kaa{E) read for two reservoirs 

4jKa.{E) = -n{E)J^{E)K^{E) - n{E + aV)J^o.{E)Ka^{E + aV) , (486) 
dA 

■^K^dE) = -n{E)^{JL + JR){E)K^a{E) - 

-niE + aV)Jo,c.{E)^{KL + KR)iE + aV) . (487) 

These RG equations can not be solved analytically but some general features can be studied. 
First, in the weak coupling regime V > Tk, the RG equations (|484p and (|485p reveal how 
the exchange couplings are cut off by the voltage. Independent on whether the diagonal or the 
nondiagonal exchange coupling is considered, there are terms n(E) which are independent of 
the voltage, and other terms ~ n{E±V) which contain the voltage. This means that for E — 0, 
the diagonal as well as the nondiagonal coupling are finally cut off by the spin relaxation rate F 
and not by the voltage. However, this does not lead to a problem or an increased conductance 
because the current rate F^ as well as the current vertex J/ are cut off by the voltage, as can 
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be seen from ()480p and (|48ip . Therefore, the logarithmic increase of the exchange couplings 
between V and F does not influence the conductance considerably. However, it is important 
to notice that this does not mean that one can omit the spin relaxation rate F from the RG 
equation for the couplings. In this case, the couplings would diverge at the Kondo temperature 
and, as a consequence, also the current vertex and the conductance since the cutoff functions 
are smooth functions. Therefore, although the rate F does not appear in the final result for the 
conductance for V > Tk, it is still important to have it in the RG equations for the couplings. 

In the strong coupling regime T, y < Tk, the couplings do not diverge because the rate 
F cuts off the RG flow roughly at the Kondo temperature. Since no other energy scale remains, 
F is expected to approach Tk when A is below the Kondo temperature. To analyse this 
question let us consider the simplest case E = V = T = and all exchange couplings to be the 
same Jaa' — J- Since h{0) = 0, we get the set of differential equations 



Analysing this nonlinar set of differential equations numerically, one finds that the solution 
is unstable against exponentially small changes in the initial condition for F. Thereby, the 
crucial problem is not the elimination of the linear terms in Aq of the initial condition (|472p 
for F, but to determine the initial condition such that F saturates at the Kondo temperature. 
Furthermore, if one studies the influence on the differential conductance, one finds that the 
conductance at T ^ V ~ can be tuned to any value depending on the initial condition 
for F. Due to this ambiguity, we have fitted the initial condition ioT F a.t T = V = E = 
such that the conductance becomes universal G = Gq — 2e^//i. Using this single parameter, 
we have then solved the RG equations for arbitrary values of voltage and temperature. The 
result for the linear conductance G{T)/Go is shown in Fig. [22] together with a comparism 
to numerical renormalization group results performed by Costi |68j . The coincidence is quite 
impressive considering the fact that we have made many approximations by neglecting higher 
order terms and the frequency dependence of the vertices and the Liouvillian. Especially the 
fact that the broadening of the crossover regime comes out correctly shows that the coincidence 
is nontrivial and is not done by hand with the fitting parameter. It means that the approximate 
RG equations contain the right correlation between the broadening of the crossover behaviour 
and the final height of the conductance at T = 0. It will be an interesting question for future 
research whether the inclusion of higher order terms will avoid the instability and provide the 
saturation of F at the right scale. The conductance can also be calculated as function of voltage 
but the result falls nearly on top of the curve shown in Fig. [51] by using a simple rescaling. One 
can also determine the coefficients of the Fermi liquid relation (|458p but it turns out that the 
ratio of the coefficients comes out incorrectly by roughly a factor of 2. 



6 Summary and outlook 

In this paper we have presented a formally exact renormalization group scheme to analyse the 
time evolution and the stationary state for a generic problem of dissipative quantum mechanics: 
An arbitrary local quantum system with a small number of states coupled to several grand- 
canonical reservoirs at different temperatures and chemical potentials. The RG scheme is set 
up by an expansion in the system-reservoir coupling using a compact diagrammatic represen- 
tation which is based on a quantum field-theoretical formulation in Liouville space. It allows 
for a direct determination of the irreducible kernel which is the dissipative part of the kinetic 
equation for the reduced density matrix of the local system. This kernel determines at the same 
time the effective Liouville operator of the local system which is fed back into the RG equations 
to determine the effective propagator between the vertices. The RG scheme is set up in Laplace 
space allowing for an analysis of the time evolution in the presence of inhomogeneous boundary 
conditions leading to a nonequilibrium stationary state. Generically, we have shown how to 
avoid the occurence of the zero eigenvalue of the effective Liouville operator in the propagator. 




dJ 
dA 



A + F 



(488) 
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Fig. 22. Comparism of the linear conductance as function of temperature from RTRG-FS (red curve) 
with NRG results from Costi [68]. The Kondo temperature Tk is defined such that G{T = Tk) — \Go- 



This provides a generic proof that the RG flow is ultimatively cut off by relaxation and dephas- 
ing rates. We have shown how to calculate the average of arbitrary observables. Correlation 
functions can also be calculated with similiar schemes which have not been shown here. 

The RG equations represent an infinite hierarchy of coupled differential equations for effec- 
tive vertices with an arbitrary number of reservoir field operators. The important input of this 
scheme is a cutoff-dependent correlation function for the reservoirs which can be chosen arbi- 
trarily. However, to be able to truncate the infinite hierarchy in a systematic way, we discussed 
several criteria how an appropriate cutoff function has to be defined. To avoid the problem of 
the occurence of the zero eigenvalue in the effective propagators we developed a generic two 
stage procedure. First one integrates out the symmetric part of the reservoir distribution func- 
tion by using perturbation theory. The result is used as an input for the following continuous 
RG fiow with a cutoff dependent antisymmetric reservoir distribution function. The antisym- 
metry guarantees that the zero eigenvalue can not occur in the effective propagators between 
the vertices. Connected with this property is the fact that only the vertices averaged over the 
Keldysh indices occur in the RG scheme which simplifies the calculation. For the choice of the 
cutoff dependence of the reservoir correlation function within the continuous RG scheme, we 
proposed to integrate out the poles and branch cuts in the complex plane by defining the cutoff 
on the imaginary frequency axis. This has three advantages. First, high energy scales are inte- 
grated out systematically providing a high-energy cutoff for the frequency integrations in the 
effective perturbation theory. Secondly, with respect to the real frequencies, all energy scales 
are considered in each step of the RG procedure so that relaxation and dephasing rates are 
generated from the very beginning. Third, the numerical stability of the RG flow is improved 
considerably because the resolvents occuring on the r.h.s. of the RG equation contain large 
imaginary parts in the denominator which are of the order of the cutoff parameter. 

The choice of the cutoff function is particularly simple for the special case of fermionic 
reservoirs with a flat density of states. In this case, the cutoff is deflned by integrating out the 
Matsubara poles of the Fermi functions step by step. All integrations over the real frequencies 
can then be calculated analytically and the RG scheme can be set up on the Matsubara axis 
similiar to exact RG schemes for equilibrium problems. However, no analytic continuation to 
real frequencies is necessary to calculate transport properties. Furthermore, one has to consider 
a whole set of Matsubara axis shifted by the real part of the Laplace variable (to account for 
the time evolution) and by multiples of the chemical potentials of the reservoirs (to account 
for finite bias voltages). For different temperatures, each reservoir defines its own Matsubara 
axis. We have shown that the stability of the RG fiow is expected to be very good because all 
imaginary parts in the denominators of the resolvents are strictly positive so that no cancella- 
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tions can occur. From the real part of the denominators one can generically see that enhanced 
renormahzations are expected if the resoiicLiic6 condition E — f^i...n — 

ft, = is fulfilled, where E 

is the real part of the Laplace variable, /Ii...„ consists of the sum of arbitrary multiples of the 
chemical potentials of the reservoirs, and h are the renormalized oscillation frequencies of the 
time evolution of the reduced density matrix of the local quantum system (e.g. magnetic field, 
level spacing, charging energy, etc.). 

For the case of fermionic reservoirs with a flat density of states we have shown how to solve 
the RG equations analytically in the weak coupling regime. Following Ref. [13], we considered a 
system with generic spin and/or orbital fluctuations and derived the final result for the effective 
Liouvillian and the average of an arbitray observable up to order J^. Here, Jc is the order of 
magnitude of the effective vertex in leading order at scale Ac = Max{|£'|, |/2i...„|, \h\]. For the 
effective oscillation frequencies we found a logarithmic term ~ ln(ylc/|-E — — ft+jr'l), where 
r is the relaxation or dephasing rate corresponding to the oscillation frequency h. At resonance 
E — p,i...n — h = 0, the logarithmic term is cut off by F and an enhancement ~ ln( Jc) is 
obtained. In arbitrary order one expects enhancements ~ Jci^n Jc)^^^ with fc > 2. As shown in 
Ref. [13], similiar features are generically expected for the relaxation and dephasing rates and 
for averages of observables like the current, but the order of magnitude is one power less in Jc, 
i.e. one expects enhancements ~ Jc (In Jc)*^^^ with fc > 3 at resonance. For the special case of 
the nonequilibrium Kondo model at flnite voltage V , this leads to the well-known logarithmic 
enhancements ~ \n{V/\V — h + ir\) for the differential conductance a,t V ^ h [47i36j45j. 

For strong coupling no systematic truncation scheme can be developed but we have shown 
that the presence of effective relaxation and dephasing rates cuts off the RG flow such that the 
coupling vertices will not diverge and stay of order one when the cutoff reaches the scale of the 
Kondo temperature (where usual poor man scaling equations diverge). We proposed a minimal 
set of RG equations for the strong coupling regime by truncating the RG equations and using 
a leading-order parametrization of the vertices. Preliminary results have been presented for 
the isotropic Kondo model, where a comparism to numerical renormalization group methods 
showed a promising coincidence for the linear conductance as function of temperature. However, 
within the lowest order approximation scheme, the RG equations showed an instability against 
exponentially small changes in the initial condition for the relaxation rate F. Therefore, a fltting 
procedure was necessary to fix the conductance at T = y = to the universal value. It remains 
an interesting subject for future research to include higher-order terms and study the stability 
of the RG equations against the choice of the initial conditions in the strong coupling regime 
in more detail. 

Besides the strong coupling regime interesting open questions for the future concern es- 
pecially the study of charge fluctuations, the time evolution, frequency-dependent density of 
states and bosonic systems. Especially for systems where the density of states has branch cuts 
in the upper half of the complex plane which do not coincide with the Matsubara axis, the RG 
equations get a different structure compared to fermionic systems with a flat density of states. 
Although the cutoff can be always defined by cutting off the imaginary part of all branch cuts, 
the presence of several branch cuts changes the result of the integration over the real frequen- 
cies. As a consequence, not all resolvents of the RG equations will contain the cutoff and the 
numerical stability of the equations has to be studied in more detail. However, the general two 
stage procedure and the systematic solution in the weak coupling regime can be used as well. 
Interesting systems in this connection are bosonic systems with particle number conservation 
(e.g. cold atom gases), where the condition w > > /j-q has to be fulfilled for the reservoir 
energies, or superconducting reservoirs where the density of states has a gap. Systems, where 
the frequency dependence of the density of states is analytic (besides possible poles lying at 
very high energies) or when the branch cuts can be turned into the lower half of the complex 
plane, can be described by RG equations with a similiar structure as the one for a flat density 
of states. The only difference is that the first step of the RG flow, where the symmetric part 
of the distribution function is integrated out, has to be considered again because the frequency 
integrals are different. 
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A Appendix: Commutation of dot and reservoir operators 

In this appendix we prove that all signs occuring from interchanges of fermionic dot and reser- 
voir operators cancel exactly the signs emerging from the prefactor rji . . .r]„ and the reversed 
sequence : a„ . . . ai : (compared to : ai . . . a„ :) in the form ()31|1 of the coupling 

y = ^77i...77„ : a„a„_i...ai : gi2...„ • (489) 
n! 

This means that after having defined the coupling vertex 51... „ via this equation, one can use the 
simpler form (I29p and consider dot and reservoir operators as commuting objects. This property 
can be proven for any average of a sequence of t^-operators over the reservoir distribution Pres 

{V{h)V{t2)...V{U))p^^, , (490) 

where V{t) is the interaction picture with respect to Hj-^g+Hg. These are the expressions which 
occur for the dynamics of the reduced density matrix of the dot. For observables one of the 
V^-operators is replaced by the observable but the proof remains the same since the observable 
is of the same form as V . 

We start the proof for n = \. The coupling is of the form V ~ 7710191 and (|490p becomes 
(time-dependence not indicated) 

771 ?72 •■•7?r (01510252 •••ar.9r)p,e= • (491) 

This average is only nonzero if r is even and if rji — 0, so that an equal number of annihilation 
and creation operators is present. Thus, taking all reservoir field operators to the right (starting 
with Or-), we get a total sign (—1)'"''^ from commuting them through the dot operators, and 
(|49T1) reads 

Vim-'-Vri-^Y^^ 9i92---9r{aia2...ar)p^^^ • (492) 

Using rji = and r even, we get 771 ... 77^ (—1)''/^ — (— 1)''/^(— l)*"/^ = (—1)'" = 1 for the 
prefactor, and all sign factors have cancelled leaving an independent product of dot and reservoir 
operators. 

For 71 > 1, we formally imagine 51. ..n to consist of a product 5152 ■ ■ ■ 9n of fermionic dot 
operators gi (which is allowed to determine the fermionic sign factors). From the proof for 
71 = 1, we know that if V is of the form 

^ = -7 111112 ■■■Vn ■ (0l5l)(0252) . . . (an5n) ■ , (493) 
m 

all sign factors cancel in the end. Moving all reservoir field operators in (|493p to the left (starting 
with o„), each time involving an even number of permutations of fermionic operators, we obtain 

1 , , 

y = — 7717/2 ■ ■ - Vn ■ CLnan-l . . . Ol : 5152 ■ ■ ■ 9n , (494) 



which is precisely of the form (j489p after replacing 51 ... 5„ by 51. .. 
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B Appendix: Coupling operator in Liouville space 

In this appendix we prove Eq. ([75]) for the coupling vertex, originally defined by (|7ip with V 
given by (|29p . i.e. we have to prove 

Lv ^ [gi...n : ai ... a„ :,•]_ = ^ a^-f" : ■ • ■ : ■ (495) 

Since, due to the definition (j77p of the coupling vertex, all pi in this expression are identical 
to a common index p, we distinguish the two cases p = ± and add them up finally. Thereby, 
p = + produces the first term of the commutator in (|495p . and p = — the second one. 
For p — +, we get by acting on an arbitrary operator A 

G+;;„+ : J+ ...J+ : A = 51.. .„ : ai . . . a„ : A , (496) 

where we have used the definitions ((78|l . ([77|. and (|76|l of the sign-operator, the coupling vertex 
and the reservoir field operators in Liouville space. As a consequence, the first term of the 
commutator in (j495p is produced for pi = p = +. 

For p — we first consider bosons, where cr^ ■■^ = 1. We obtain 

cr"- -" G]^;;„" : J]" . . . J,7 : A = -A : a„ . . . oi : 5i...„ 

= -Agi,,,n : ai . . . a„ : , (497) 

coinciding with the second part of the commutator of (|495p . 

For p — ~ and fermions, we distinguish between n even and n odd. For n even, we have 
cr^ ■■^ = ( — 1)"/^ and the same result is obtained for 

(^^"^ Gi.Wa ■■ Ji ---Jn ■■ A = -(-l)"/2 A : a„ . . . ai : gi,,,^ 

= -Agi...„ (-1)"/^ : a„...ai := -A.gi...„ : ai . . . a„ : . 

Here we have used in the second equality the fact that operators from the reservoirs and the 
quantum system are considered to commute (a property which holds finally if all expressions are 
averaged with respect to the reservoir degrees of freedom, as explained in detail in Sec. 12. II and 
appendix A). For n odd, we use a~ ^^~ = (— 1)("~^'/^(t~, and, since a~ cancels with another 
a~ occuring in the definition ()77p of the coupling vertex, we obtain again 

Gi:::n ■.J^...J-:A= _(_l)(»-i)/2 ^ : a„ . . . ai : gi...„ 

= -Agi„,n (-l)(""i)/2 : a„ . . . ai : = -Agi,,,^ : ai . . . a„ : , 

and the proof of Eq. (|495p is complete. 



C Appendix: Kernel properties 

Here we prove the property (jllSp of the kernel S{E) 

{S{E)y = -E{-E*) . 
We use the diagrammatic representation (jllOp for the kernel 

m) ^ (n7).„, G^^^^G G^^^^G , (499) 

and calculate the transformation {S{E)Y using the relation {ABY = A'^B'^ (see (HZ])). With 
the help of the properties (|92p and (|94p , and renaming all summation variables by p ^ —p and 
•q —77, we obtain the following replacements in (|499p 

l7v -7f^' > (500) 

G?^.;;f" - -<y---G^:::'' , (501) 
1 11 ^^^^^ 

E — Xi — Lq E* + Xi + Ls —E* — Xi — Ls 
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where we have used the form (|103|) for the contraction and the upper (lower) sign corresponds 
to bosons (fermions). Note the fact that rj — )■ —77 implies also x — rj[uj + — > —a; and, 
consequently, Xi —Xi. Besides a factor ± for each contraction and a sign operator a ■ •^ 
for each vertex (which can give rise to a minus sign only for fermions) , we obtain the replacement 
E — > —E* and a total minus sign since the number of vertices is by one larger than the number 
of resolvents in (|499p . Thus, to complete the proof of (|498p . we have to show for fermions that 
the occurence of the sign operators cancels the minus sign for each contraction. To show this, 
we shift all sign operators through the vertices G to the left by using the property ([5T|) . As 
a consequence, all sign operators can be taken together to a single sign operator of the form 
(J where the number of minus signs in this expression corresponds to the total number Z 

of reservoir field operators of the diagram. Since Z = 2r must be an even number, the total sign 
operator is identical to a ■ •^ = (—1)'' according to the definition (j78|) . Since r is precisely 
the number of contractions, this sign factor cancels exacly one sign for each contraction. 

D Appendix: Invariance of symmetry relations 

Here we show that the symmetry properties (|259p - (|264p are invariant under the RG flow. 
Obviously, they are fulfilled initially. So we assume that we can use them on the r.h.s. of the 
RG equations (|246p and (|253p . and prove them for the l.h.s. 

The (anti-)symmetry relations (|259p and (|260p follow from the diagrammatic rules stated 
after (|227p . Since we sum in the RG equations over all permutations of the external indices 
and assign a corresponding fermionic sign, the effective vertices become automatically (anti-) 
symmetric. 

The property (|26ip of conservation of probability is also trivially fulfilled, since by acting 
with Tts on the RG equation and summing over the external Keldysh indices, we get zero since 
the first vertex already fulfils this property. Note that by summing over all Keldysh indices, the 
first vertex is just averaged over the Keldysh indices since the contractions do not depend on 
the left Keldysh index. 

The properties (I262p - (l264p . originally related to the hermiticity of the Hamiltonian, follow 
immediately by taking the c-transform (defined in Eq. (I9ip ) of the RG equations and applying 
a proof in analogy to the one presented in Appendix C. The only difference is that, for the 
vertices, the total sign operator a ■ obtained from shifting all sign operators to the left 
using (|8ip . contains also the minus signs from the external indices. However, by using ()80p . 
this sign operator can be split into two sign operators, one containing the minus signs from the 
external vertices (leading to the sign operator on the r.h.s. of the properties (|262p - (|264p '). the 
other containing the minus signs from the internal vertices which cancels the signs from all the 
contractions, see Appendix C. 

E Appendix: Leading order RG equations 

Here we prove various properties for the vertices G^iy, G^^',, and G^^, together with corre- 
sponding properties for the vertices R^-^y and Rfy , which are needed in Sec. 14.41 For notational 
simplicity, wc omit in the following the index (1) for all vertices. 

First, we prove that the RG equations (|324p and (|336p for the vertices G^^! and gJ\! are 
consistent with the RG equation (|33ip for the vertex giy using the relations (|329p and (|335p 
between G and g. From p29p and (I33ip we get 

^PP _ 1 (r'PPr'PP nPP r<pp\ fKr\'i\ 

by simply acting on an operator and using G^^, — — G^fj^. Summing this equation over p gives 
the correct RG equation for Gw 

A.G^r = i ^ (G?^Gf^, - Gff.G^O (504) 
p 
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_ 1 \ ^ Ir-iPP^PP ^PP r-iPP I i^PPr^PP rtPP r~iPP\ (r,c\r,\ 

— — l'-^12'-^21' ~ '-^1'2'-^21 + '-'12'-'21' ~ *-^l'2'-^21'' (.OUOj 

P 

- \ Y.^GllG-,y-Gr,G^l) (506) 

= \ {Gl2G2y - Gv2G-2l) , (507) 

where p — —p, and we have used GffG^^, = G^yGi2 — G^fjG^^ to arrive at (|505p (note that 
2 is a summation index which can be changed to 2). 

Multiplying (|503p with p and then summing over p gives the RG equation for Gn' 

ipG'mi - PG\\GZ + pGZG% " vG^G^) (508) 

p 

= \ Y.^GZG2v-vGv2G^-^ (509) 
p 

= j{Gl2G2y -GV2G2,) , (510) 

and with a similiar proof we can show the same with G <-!■ G on the r.h.s. 

Next we show that the RG equation p25p for R^-^y is consistent with the RG equation (|332p 
for the vertex riy using the relation (|330p between R and r. As usual this holds only if one acts 
with the trace over the quantum system from the left. Doing this and acting with the r.h.s. of 
P25p on an arbitrary operator A, we get the correct RG equation for 

jTTs{Rl2G2y - RV2G2M (511) 

1 

= ^Tr5(ri252i' - fi'252i - 52i'?'i2 + g2irv2)A (512) 

= ^Trs(ri252i' - ^'2521 + -5'i2?'2i' - gi'2r2i)A (513) 

= i-^TrsrivA ^ -^TvsitivA . (514) 

For the special case of the current vertex, we prove that the ansatz p37p for the current 
vertex fulfils the correct RG equation. For this we need the property 

TrsGii'G22' = - Trs G22'Gii' , (515) 

which is based on applying several times the property TrgGn' — (see (|85|)) or 

TrsGf?, = -TrsGff, , (516) 

with p = —p. As a consequence we also get 

TrsGT,,GZ = ^TrsG^G^, = -TrsG5f,Gff, = TtsG%G% , (517) 

and (|515p can be proven in the following way 

Tr5Gn'G22' = Trs E^'^n'^'s'' (518) 

pp' 

= Trs Y.PGfvGll, + Trs E^^^n'^i' (519) 

p p 

= Trs Y.PG%G% + Trs Y.PG%G\l> (520) 

p p 

= Trs Y.PGf2'Giv = -Trs Y.PG^22'Giv - -TrsG22'Giv , (521) 
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where we have appUed (|517p in ()519p and (|52ip . 

With the help of (I515P , we can now easily proof the correct RG equation for the current 
vertex using the ansatz (|337p 

^Trs /7i, = cI,,Tts ^Gw = jcJyTrs {G^^G-^y - Gy^G^,) (522) 
= ^c7i,Tr5Gi2G2i, = ^Trs [ci^G,2G-2y + cl^,G^2G-^y] (523) 



= ixrs [c1^Gi2G-^y - cl^,G-^yGi2] (524) 
= ixrs {c^^GisG.i, - cl^Gy^G-^^] (525) 
= ixrs (/72G21, - /7,2G2i) , (526) 

where we have used c^/ = + c^^/ hi (|523p and applied (|515p to arrive at (|524p . 

We study now the RG equation for the vertices G\y and R\y which are defined by p49p . 
For both cases we arrive with the same proof at the correct RG equation (|348p in the following 
way (we have chosen R for the proof) 

= ^*2yl'^'"'^ |(^13G32 - i?23G3i)G2i' + 

+Ri2{GnG-^y - GrsGga) - (1 ^ l')} (528) 

= -«^Trs {i?12(G23G3i, -GP3G32) + 

+ i?13G32G2i. - i?23G3iG2r - (1 ^ l')} (529) 

= ^Trs {^i2Gfi! + i?f2^G2i, -(1^1')}- 

-«^Trsi?23 {G31G21, - G31G2P - (1 ^ 1')} , (530) 

where we have interchanged 2 ^ 3 in the first and third term of the r.h.s. of (|528p and have 
added and subtracted the term ^23G3iG2i/ to (|529p . The first term of (|530p gives the correct 
RG equation for R\y. The second term is zero as can be seen from 

^23 IG31G21' — G31G21' — G31/G21 + G31/G21I (531) 
= ^23 IG31G21/ — G31G21' + G2i'G3i — G2i'G3i| (532) 

= i?23 E [pG'lGi': - P'Gf.Gi': + pG^v G{f - p'Gf^. Gf } (533) 

pp' 

= 2 i?23 E {pG^Gfr + PG% Gf^ } (534) 

p 

= 2 i?23 E {G'IG% - G% Gf^ } = , (535) 
p 

where we have interchanged 2 <-> 3 in the third and fourth term of (j53ip and used i?23 — —Rzi- 
Finally, we prove the properties (|333p and (|334p . (511')^ = _and (rii/)J' = r\,\ follow 
trivially from the RG equations (|33ip and p32p . For the vertices Gw and Rw we use the 
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matrix representation 

(Gii')sis;,s2S2 = (5ii')siS2<5s;s^ - <5siS2(5ll')s^s'l J (536) 

— % 

{Rlv)s^s'^,S2s'^ = l^{{riv)s,S2Ss[s'^_+ Ss,s2{'riv)s'^s'^) , (537) 

and get for gJ^, 

{G\r)siS[,S2s'^ = {Gll'Ts^s'^^Sis'^ = (5ll')s2s/s'iSi ~ '^SiS2(5ll')s;4 

~ {9ll')siS2^s[s'^ ~ ^SiS2{9ll')s'^s[ 

= {9l'l)siS2^s{s'^ ~ '5siS2 (ffl'l)si,s'j 

— (Gl'l)sis'j,S2S2 J (538) 

which proves ([333]). In the same way one proves ([333]) for the vertex Rw. 
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